Next Article in Journal / Special Issue
Strangeness Production in Nucleus-Nucleus Collisions at SIS Energies
Previous Article in Journal
Effects of Scattering of Radiation on Wormholes
Previous Article in Special Issue
A Phenomenological Equation of State of Strongly Interacting Matter with First-Order Phase Transitions and Critical Points
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prospects of Constraining the Dense Matter Equation of State from Timing Analysis of Pulsars in Double Neutron Star Binaries: The Cases of PSR J0737 ‒ 3039A and PSR J1757 ‒ 1854

The Institute of Mathematical Sciences (IMSc-HBNI), 4th Cross Road, CIT Campus, Chennai 600 113, India
Universe 2018, 4(2), 36; https://doi.org/10.3390/universe4020036
Submission received: 4 December 2017 / Revised: 8 February 2018 / Accepted: 8 February 2018 / Published: 12 February 2018
(This article belongs to the Special Issue Compact Stars in the QCD Phase Diagram)

Abstract

:
The Lense-Thirring effect from spinning neutron stars in double neutron star binaries contributes to the periastron advance of the orbit. This extra term involves the moment of inertia of the neutron stars. The moment of inertia, on the other hand, depends on the mass and spin of the neutron star, as well as the equation of state of the matter. If at least one member of the double neutron star binary (better the faster one) is a radio pulsar, then accurate timing analysis might lead to the estimation of the contribution of the Lense-Thirring effect to the periastron advance, which will lead to the measurement of the moment of inertia of the pulsar. The combination of the knowledge on the values of the moment of inertia, the mass and the spin of the pulsar will give a new constraint on the equation of state. Pulsars in double neutron star binaries are the best for this purpose as short orbits and moderately high eccentricities make the Lense-Thirring effect substantial, whereas tidal effects are negligible (unlike pulsars with main sequence or white-dwarf binaries). The most promising pulsars are PSR J0737 − 3039A and PSR J1757 − 1854. The spin-precession of pulsars due to the misalignment between the spin and the orbital angular momentum vectors affect the contribution of the Lense-Thirring effect to the periastron advance. This effect has been explored for both PSR J0737 − 3039A and PSR J1757 − 1854, and as the misalignment angles for both of these pulsars are small, the variation in the Lense-Thirring term is not much. However, to extract the Lense-Thirring effect from the observed rate of the periastron advance, more accurate timing solutions including precise proper motion and distance measurements are essential.

1. Introduction

Timing analysis of binary pulsars leads to the measurement of pulsar’s spin, Keplerian orbital parameters (the orbital period P b , the orbital eccentricity e, the longitude of periastron ω , the projected semi-major axis of the orbit x p = a p sin i , where a p is the semi-major axis of the pulsar orbit and i is the angle between the orbit and the sky-plane, and the epoch of the periastron passage1), as well as post-Keplerian (PK) parameters like the Einstein parameter ( γ ), Shapiro range (r) and shape (s) parameters, the rate of the periastron advance ω ˙ , the rate of change of the orbital period P ˙ b , the relativistic deformation of the orbit δ θ , etc. [4]. Sometimes, the Shapiro delay is parametrized differently, with parameters h 3 and ς [5], or with r and z s [6], and these parameters can easily be expressed in terms of conventional parameters s and r.
Measurement of PK parameters leads to estimation of masses of the pulsars and their companions, as well as tests of various theories of gravity [7,8]. Note that, in principle, measurements of only two PK parameters are enough to extract two unknowns, i.e., the masses of the pulsar and the companion, while measurements of more than two PK parameters lead to tests of gravity theory through consistency. However, the uncertainty in measurement is not equal for every PK parameter, and usually the two most accurate parameters are used to obtain the best mass estimates (e.g., Figure 1 of [9]). Additionally, relativistic binary pulsars have the potential to constrain the Equation of State (EoS) of matter at extreme densities. Measurements of masses of two pulsars being around 2 M have already ruled out many soft EsoS [10,11]. On the other hand, the recent analysis of the gravitational wave event GW170817 from the merger of two neutron stars ruled out some extremely stiff EsoS [12,13]. Still, a large number of EsoS is allowed, including several hybrid [14] and strange quark matter [15] EsoS in addition to standard hadronic ones. Therefore, further progress on this issue is essential, and that is expected in the near future as, on the one hand more, neutron star-neutron star mergers are expected to be detected in future runs of the advanced LIGO, and on the other hand, many binary pulsars are being timed regularly and accurately. Besides, upcoming radio telescopes like MeerKAT and the Square Kilometre Array (SKA) will lead to significant improvement in pulsar timing.
Among all of the PK parameters, in the present article, I concentrate mainly on the periastron advance, which has the potential to constrain the dense matter equation of state due to the Lense– Thirring effect. Neutron star-neutron star binaries or ‘Double Neutron Star’ (DNS) systems are the best for this purpose as these systems are compact enough to display effects of strong field gravity, yet wide enough to have negligible tidal effects. Moreover, due to the high compactness of neutron stars, spin-induced quadrupole moments are also negligible. On the other hand, as these neutron stars are rapidly spinning (spin periods of most of the pulsars in DNSs are less than 100 milliseconds), the Lense-Thirring effect is significant. That is why I concentrate on DNSs in this article, although the mathematical formulations are valid for any kind of general relativistic binaries with the possibility of additional effects depending on the nature of the objects in the binaries. There are about seventeen DNSs known at the present time, including the first discovered binary pulsar: the Hulse–Taylor binary, PSR B1913+16. For one DNS, PSR J0737 − 3039A/B, both members are radio pulsars, and the system is known as a double pulsar 2.
For a long period of time, the double pulsar was the most relativistic binary. However, the recently- discovered DNS, PSR J1757 − 1854 shows even stronger effects of general relativity in some aspects, i.e., P ˙ and γ [17], due to its high eccentricity combined with a small orbital period. In fact, its eccentricity is 6.9-times and orbital period 1.8-times larger than those of PSR J0737 − 3039A/B. It is noteworthy that ω ˙ is larger for PSR J0737 − 3039A/B than PSR J1757 − 1854 where the contribution of eccentricity is less dominant; see Lorimer and Kramer [4] for expressions for PK parameters.
PSR J1906 + 0746 is another DNS having the second smallest value of P b , just after PSR J0737 − 3039A/B. However, its P b is slightly (1.1-times) and e is significantly (7.1-times) smaller than those of PSR J1757 − 1854. That is why although PSR J1906 + 0746 has the third largest value of ω ˙ , other PK parameters are small, even smaller than most of the other DNSs. Finally, the latest discovered DNS, PSR J1946 + 2052, has broken all the records by having the smallest value of P b and the largest value of ω ˙ [18].

2. Precession in Double Neutron Star Binaries

Precessions (of both the spin and the orbit) of neutron stars in DNSs are very important. Following Barker and O’Connell [19], the rate of the change of the unit spin vector ( s a ) of a spinning neutron star (a) in a binary can be written as:
s ˙ a = Ω s a × s a ,
where the angular spin-precession frequency Ω s a can be written as:
Ω s a = Ω sPN a + Ω sLT a .
Ω sPN a is the angular spin-precession frequency of the neutron star due to the space-time curvature around its companion and can be calculated within the ‘Post-Newtonian’ (PN) formalism. Ω sLT a is the angular spin-precession frequency of that neutron star due to the Lense-Thirring effect of its spinning companion. We can further write:
Ω sPN a = A PN a k ,
where k = L / | L | is the unit vector along the orbital angular momentum, L is the orbital angular momentum, and:
Ω sLT a = A LT a s a + 1 3 ( k . s a + 1 ) k .
Here, if a is the neutron star under consideration, a + 1 is its companion. The amplitudes in Equations (3) and (4) are given by:
A PN a = G c 3 2 / 3 n 5 / 3 ( 1 e 2 ) M a + 1 ( 4 M a + 3 M a + 1 ) 2 ( M a + M a + 1 ) 4 / 3 ,
A LT a = G c 3 β s a + 1 n 2 ( 1 e 2 ) 3 / 2 M a + 1 2 2 ( M a + M a + 1 ) ,
where G is the gravitational constant, c is the speed of light in a vacuum and n = 2 π / P b is the angular orbital frequency. M a is the mass; I a is the moment of inertia; P s a is the spin period of the a-th neutron star; and
β s a + 1 = c I a + 1 G M a + 1 2 · 2 π P s a + 1 , β s a = c I a G M a 2 · 2 π P s a .
For pulsars in DNSs, the Lense-Thirring term ( Ω sLT a ) can be ignored. Even for the case of the slow pulsar (B) of the double pulsar, where the Lense-Thirring effect due to the spin of the fast pulsar (A) contributes to Ω s B , I find A PN B = 5.07481 deg y 1 and A LT B = 3.48865 × 10 5 deg y 1 (using the values of the parameters given in Table 1). This leads to Ω sLT B 6.97731 × 10 5 deg y 1 using the fact that s A is almost parallel to k [20,21], and hence, Ω s B = 5.07474 deg y 1 , which is close to the observed median value of Ω s B = 4.77 deg y 1 [22]. It is unlikely that companions of other pulsars would be much faster than PSR A, and hence, Ω sLT a can be neglected. Therefore, Equation (1) becomes:
s ˙ a = A PN a k × s a = A PN a sin χ a u ,
where χ a is the angle between k and s a and u is a unit vector perpendicular to the plane containing k and s a in the direction given by the right-hand rule. Therefore, the measurement of the spin-precession of a pulsar helps estimate χ a when other parameters involved in Equation (5) are already known.
Similarly, neglecting the tidal and the spin-quadrupole effects, as well as the spin-spin interaction, the orbital angular precession frequency can be written as:
Ω b = Ω bPN + Ω bLT ,
where Ω bPN and Ω b LT are contributions from the space-time curvature and the Lense-Thirring effect, respectively. Note that both members of the binary contribute to each term. Ω b leads to the precession of both the Laplace–Runge–Lenz vector A , as well as L , but none is directly observable. Damour and Schafer [23] first studied the manifestation of Ω b in terms of observable parameters by decomposing Ω b as:
Ω b = d φ a d t h a + d ω a d t k + d i d t Υ a ,
where h a is the unit vector along the line-of-sight, i.e., from the Earth to the a-th neutron star (the pulsar), Υ a = h a × k | h a × k | is the unit vector along the line of the ascending node, φ a is the longitude of the ascending node of the a-th neutron star and i is the angle between the orbit and the sky-plane. Damour and Schafer [23] has also given:
φ ˙ a = 1 sin 2 i [ Ω b · h a cos i ( Ω b · k ) ] ,
ω ˙ a = 1 sin 2 i [ Ω b · k cos i ( Ω b · h a ) ] ,
d i d t = Ω b · Υ a .
As already mentioned, ω ˙ a is the parameter of interest as it has the potential to put constraints on the dense matter EoS. Therefore, I explore properties of ω ˙ a in the next and subsequent sections.

Periastron Advance

The observed rate of the periastron advance ( ω ˙ a , obs ) of the a-th member of a relativistic binary can be written as:
ω ˙ a , obs = ω ˙ a + ω ˙ Kop a = ω ˙ PN a + ω ˙ LT a + ω ˙ LT a + 1 + ω ˙ Kop a ,
where ω ˙ PN a is due to the space-time curvature caused by both members of the binary, ω ˙ LT a is due to the Lense-Thirring effect of the star a and ω ˙ LT a + 1 is due to the Lense-Thirring effect of the star a + 1 . These three terms together come into the expressions of ω ˙ a (Equation (12)) caused by the precession of the orbit. ω ˙ Kop a is a secular variation due to the gradual change of the apparent orientation of the orbit with respect to the line-of-sight due to the proper motion of the barycenter of the binary [24]. The expressions for different terms are:
ω ˙ PN a = 3 β 0 2 n 1 e 2 [ 1 + β 0 2 f 0 a ] ,
ω ˙ LT a + ω ˙ LT a + 1 = 3 β 0 3 n 1 e 2 ( g s a β s a + g s a + 1 β s a + 1 ) ,
ω ˙ Kop a = K cosec i μ α cos φ a + μ δ sin φ a , K = 502.65661 ,
where only the first and the second order terms are retained in the expression of ω ˙ PN a . μ α and μ δ are the proper motion of the barycenter of the binary (which is measured as the proper motion of the visible object) in the right ascension and the declination respectively, both expressed in units of milliarcseconds per year. All other parameters are in SI units. The parameters introduced in Equations (15)–(17) are as follow:
β 0 = ( G M n ) 1 / 3 c ,
f 0 a = 1 1 e 2 39 4 X a 2 + 27 4 X a + 1 2 + 15 X a X a + 1 13 4 X a 2 + 1 4 X a + 1 2 + 13 3 X a X a + 1 ,
where X a = M a / M , M = M a + M a + 1 is the total mass of the system. β s a is defined in Equation (7), and
g s a = X a 4 X a + 3 X a + 1 6 ( 1 e 2 ) 1 / 2 sin 2 i × ( 3 sin 2 i 1 ) k · s a + cos i h a · s a = X a 4 X a + 3 X a + 1 6 ( 1 e 2 ) 1 / 2 sin 2 i × ( 3 sin 2 i 1 ) cos χ a + cos i cos λ a
where λ a is the angle between h a and s a (see Figure 1). It is obvious from Equation (20) that the maximum value of g s a occurs when s a is parallel to the vector ( 3 sin 2 i 1 ) k + cos i h a , giving:
g s a , m a x = 3 + 1 sin 2 i 1 / 2 X a 4 X a + 3 X a + 1 6 ( 1 e 2 ) 1 / 2 .
If s a k , i.e., χ a = 0 and λ a = i :
g s a , = X a 4 X a + 3 X a + 1 3 ( 1 e 2 ) 1 / 2 .
Using the expression of g s a , in Equation (16), I get ω ˙ LT a , , and using the expression of g s a , m a x , I get ω ˙ LT a , m a x . The negative sign in the Lense-Thirring term (Equation (16)) implies the fact that actually ω ˙ LT a , m a x is the minimum value. However, depending on the values of i, χ p , and λ p , g s a (Equation (20)) can be negative making ω ˙ LT a (Equation (16)) positive. Note that the Lense-Thirring effects from both the pulsar and the companion come in the total periastron advance rate (Equations (14) and (16)). However, if the companion is much slower than the pulsar, as in the case of the double pulsar, then the contribution from the companion can be neglected. Due to the non-detection of any pulsation, we do not know the values of the spin periods of the companions for other DNSs and cannot rule out the possibility of significant Lense-Thirring effect from those companions. However, most of the pulsars in DNSs are recycled (at least mildly), suggesting the fact that pulsars were born as neutron stars earlier than their companions. The companions, i.e., the second born neutron stars in DNSs, are expected to be slow; because even if neutron stars are born with spin periods in the range of a few tens of milliseconds up to a few hundred milliseconds, they quickly spin down to periods of a few seconds unless they gain angular momentum via mass accretion from their Roche lobe filling giant companions and become fast rotators. Such spin-up is possible only for the first born neutron star in a DNS. That is why in the present article, I ignore the Lense-Thirring effect from the companions of pulsars in DNSs.
In such a case, i.e., when ω ˙ LT a + 1 is negligible, ω ˙ LT a can be extracted by subtracting ω ˙ PN a + ω ˙ Kop a from ω ˙ a , obs (Equation (14)). One can estimate the value of I a from ω ˙ LT a if all other relevant parameters like P b , e, M a , M a + 1 , sin i , χ a , λ a are known (Equations (7), (18)–(20)). Then, from the knowledge of I a , M a and P s a , a new constraint on the EoS can be placed, as it is a well-known fact that the moment of inertia of a neutron star depends on its mass, spin-period and the EoS (see [25] and the references therein). For this purpose, DNSs with large values of ω ˙ LT a (larger than the measurement uncertainty in ω ˙ a , obs ) are preferable. Large values of the orbital eccentricity e and small values of P b and P s make ω ˙ LT a larger.
The procedure is actually a bit trickier than it seems at first. Usually, timing analysis of a binary pulsar reports the most accurate measurement of ω ˙ a , obs out of all PK parameters, and the masses of the pulsar and the companion are estimated by equating ω ˙ a , obs with ω ˙ PN a and using the second most accurately-measured PK parameter. If ω ˙ LT a is large, then this procedure would become erroneous. One should instead estimate masses using two PK parameters other than ω ˙ a , obs , and these mass values should also be accurate enough to give a precise estimate of ω ˙ PN a . Accurate measurements of proper motion are also needed to evaluate ω ˙ Kop a . Only then, ω ˙ PN a + ω ˙ Kop a can be subtracted from the observed ω ˙ a , obs to get the value of ω ˙ LT a . In short, one needs very precise measurements of at least three PK parameters, ω ˙ a , obs and two more, as well as the distance and the proper motion of the pulsar (the a-th neutron star).
Because of its small values of P b and P s , PSR J0737 − 3039A is expected to have a large value of ω ˙ LT a in spite of its low eccentricity. Similarly, as PSR J1757−1854 has a slightly larger P b , slightly smaller P s and much larger e, it is also expected to have a large value of ω ˙ LT a . For this reason, in the next section, I explore ω ˙ LT a for these two DNSs in detail. Note that, although, PSR J1946 + 2052 might be a useful system for this purpose, it is not possible to extend my study for this system due to lack of knowledge of masses of the pulsar and the companion.
One should also remember the dependence of ω ˙ LT a (via g s a ) on the orientation of s a with respect to k and h a (Equation (20)). If I assume that the a-th neutron star is the pulsar, then all ‘a’s in the subscripts of the above equations will be replaced by ‘p’s, e.g., s p , h p , etc. Most of the time, the value of g s p , is estimated in the literature. However, at least for some pulsars, χ p 0 . In such cases, the values of χ p and λ p can be measured by analyzing the change of the pulse profile shape and the polarization. Conventionally, χ p is denoted by δ , and λ p is denoted by ζ . Note that, even if χ p is a constant over time, λ p would vary due to the spin precession and can be written as [26]:
cos λ p = cos χ p cos i + sin χ p sin i cos [ Ω s p ( t T p 0 ) ] .
where Ω s p is the angular spin precession frequency (amplitude of Ω s p in Equation (2)), t is the epoch of the observation and T p 0 is the epoch of the precession phase being zero. Note that i itself might change with time as given in Equation (13).
In Table 1, with other relevant parameters, I compile values of χ p wherever available. I do not compile values of λ p , although available in some cases, as it is a time-dependent parameter. From Figure 1, it is clear that the spin would precess in such a way that it would lie on the surface of a fiducial cone with the vertex at the pulsar and having the half-opening angle as χ p . The constraint on λ p is: (i) if χ p < i , then λ p , m a x = i + χ p and λ p , m i n = i χ p (Figure 1a), (ii) if χ p > i , then λ p , m a x = χ p + i and λ p , m i n = χ p i (Figure 1b).

3. Results

For the purpose of demonstration, I compute the value of the moment of inertia ( I p ) for the two most interesting pulsars, PSR J0737 − 3039A and PSR J1757 − 1854 using the Akmal-Pandharipande- Ravenhall (APR) equation of state [27] and the RNS code3 [28]. I find that, both of these pulsars have I p = 1.26 × 10 45 gm cm 2 . I use this value in my calculations, remembering the fact that the true values of I p would be different depending on the true EoS, and we are actually seeking an answer whether it would be possible to know the value of I p by singling out ω ˙ LT p from the total observed ω ˙ p , obs .
I tabulate values of ω ˙ PN p , ω ˙ LT p , and ω ˙ LT p , m a x for all pulsars in DNSs in Table 2. I exclude candidate DNSs, i.e., PSR B1820−11 and PSR J1753 − 2240 from this calculation due to poor mass constraints. For other cases, I use the limiting mass values where actual values are unavailable. Only ω ˙ LT p , can be calculated for the systems with undetermined/unpublished values of sin i . As expected, PSR J0737 − 3039A shows the largest Lense-Thirring effect, followed by PSR J1757 − 1854. In fact, | ω ˙ LT p , m a x | / | ω ˙ PN p | is slightly larger for PSR J1757 − 1854 ( 2.92 × 10 5 ) than that for PSR J0737 − 3039A ( 2.81 × 10 5 ). This fact makes it a very interesting system, and I investigate this system in detail. The discovery paper [17] also mentioned the large Lense-Thirring effect causing significant amount of d i / d t (Equation (13)).
I find that, for PSR J1757 − 1854, | ω ˙ LT p , m a x | = 3.03170 × 10 4 deg y 1 is achieved for χ p = 3 , λ p = 99 . However, as the value of χ p for this system is not yet known (although Cameron et al. [17] obtained a constraint χ p 25 based on their simulation of the kick velocity), I vary χ p in the range of 0 60 . It is obvious that ω ˙ LT p depends more on χ p than on λ p as the absolute value of the factor (see Equation (20)) with cos χ p , i.e., 3 sin 2 i 1 , is larger than that of the factor with cos λ p , i.e., cos i . Furthermore, the measured value of sin i gives two values of i (as sin θ = sin ( π θ ) ), and I use the larger one to get a bigger range of λ p . The minimum value of | ω ˙ LT p | = 1.37597 × 10 4 deg y 1 is achieved for χ p = 60 , λ p = 36 . If I fix χ p strictly as 25 , then the maximum and minimum values of | ω ˙ LT p | are 2.81108 × 10 4 deg y 1 and 2.67663 × 10 4 deg y 1 , for λ p = 121 and λ p = 71 , respectively. In Figure 2, I show the variation of ω ˙ LT p for χ p in the range of 0 60 , and for each value of χ p , λ p varies in the range of i χ p , i + χ p . On the other hand, as PSR J0737 − 3039A has χ p < 6 [20], ω ˙ LT p always remains close to ω ˙ LT p , . The maximum and minimum values of | ω ˙ LT p | are 4.74488 × 10 4 deg y 1 at χ p = 0.65 , λ p = 91.95 ; and 4.712931 × 10 4 deg y 1 at χ p = 6.0 , λ p = 85.31 . To check the validity of my assumption that the Lense-Thirring effect from the companion could be ignored, I find that ω ˙ LT c , = 3.792899 × 10 6 deg y 1 using I c = 1.16 × 10 45 gm cm 2 , where the subscript ‘c’ stands for the companion, i.e., the slow PSR J0737 − 3039B. Varying χ c in the range of 130 150 [20], I find that the minimum value of | ω ˙ LT c | is 2.60052 × 10 6 deg y 1 . These results support my argument that the Lense-Thirring effect from the companion of a pulsar in a DNS would have a negligible contribution to the rate of the periastron advance of the binary.
As Table 2 shows, | ω ˙ LT p , m a x | for PSR J0737 − 3039A is around 1.4-times smaller than the presently published accuracy, ω ˙ p , obs = 16.89947 ± 0.00068 deg y 1 where the uncertainty is twice the formal 1 σ value obtained in the timing solution [9]. The uncertainties in other PK parameters are even poorer (Table 1 of [9]), so if we exclude ω ˙ p , obs , the mass estimates would be less accurate and should not be used to calculate ω ˙ PN p . Therefore, to achieve the goal of estimating ω ˙ LT p , lowering only the uncertainty in ω ˙ p , obs would not be enough; one needs to improve the accuracy of at least two other parameters that can be used to get masses at least as precise as the ones already published. It is impossible to improve the accuracy of the Keplerian parameter R used by Kramer et al. [9] in combination with ω ˙ p , obs to report the values of the masses, because R = M p / M c = x c / x p involves x c , the projected semi-major axis of the companion (PSR B), which is not visible presently and even when it was visible, it was not a good timer like PSR A. Solving the expressions for PK parameters, I find that the combination of s and P ˙ b is the best choice left when ω ˙ p , obs and R are excluded, i.e., gives the narrowest ranges of masses M p = 1.34 ± 0.02 M , M c = 1.251 ± 0.007 M using the published values of the uncertainties. I also find the fact that, to get masses as accurate as the ones reported by Kramer et al. [9], the uncertainty in P ˙ b should be reduced at least to 8.0 × 10 16 assuming that s and relevant Keplerian parameters ( P b , e, x p ) would also be improved by at least an order of magnitude. The published uncertainty in P ˙ b is 1.7 × 10 14 . However, the question arises whether for such an improved measurement of P ˙ b , the dynamical contribution coming from the proper motion of the pulsar (Shklovskii effect) and the relative acceleration between the pulsar and the solar system barycenter could be ignored. Using the proper motion reported by Kramer et al. [9] and the distance (1.1 kpc) estimated by Deller et al. [29] using the Very-long-baseline interferometry (VLBI) parallax measurement, I find that P ˙ b , dyn = 4.91 × 10 17 , where I have used a realistic model of the galactic potential [30]. The improvement required in the timing solution of PSR J0737 − 3039A is expected by 2030 [31,32] accompanied by an improvement in ω ˙ p , obs . Fortunately, ω ˙ Kop p is smaller than | ω ˙ LT p , | ; using the proper motion reported in Kramer et al. [9], I find that it can be at most 1.18 × 10 6 deg y 1 (I could not calculate the actual value, as φ p is unknown). However, one will still need better estimates of χ p and λ p to be able to extract the value of the moment of inertia of the pulsar from the ω ˙ LT p .
On the other hand, for PSR J1757 − 1854, both | ω ˙ LT p , m a x | and | ω ˙ LT p , | are about 1.5-times larger than the presently-published accuracy, ω ˙ p , obs = 10.3651 ± 0.0002 deg y 1 . However, even this system cannot be used at present to estimate ω ˙ LT p , not only because of the unknown values of χ p and λ p , but also because of the less accurate values of other PK parameters. If one uses γ and P ˙ b , then the uncertainties in these parameters should be improved by one and two orders of magnitudes respectively to get the uncertainty in masses of about 0.0008 M . This seems plausible, as the present accuracy in P ˙ b is only 0.2 × 10 12 [17]. However, due to the lack of measurements of the proper motion and the parallax distance, I could not calculate P ˙ b , dyn for this pulsar and only can give rough estimates based on the distance guessed from the dispersion measure using two different models of the galactic electron density, e.g., the NE2001 model [33,34] and the YMW16 model [35]. The contribution due to the relative acceleration is 2.11 × 10 15 for the NE2001 distance (7.4 kpc) and 1.65 × 10 14 for the YMW16 distance (19.6 kpc). Both of the values seem to be too large to be canceled out by the Shklovskii term. Measurements of the proper motion and the distance (using parallax) are essential to estimate P ˙ b , dyn . The measurement of the proper motion will also help estimate ω ˙ Kop p .

4. Summary and Conclusions

Measurement of the moment of inertia of a pulsar will provide yet another constraint on the dense matter equation of state. This issue has been discussed in the literature in the past (see [32,51,52] and the references therein), but unfortunately, the effect of the misalignment between the spin and the orbital angular momentum vectors has not been explored in detail mainly because these two vectors are almost parallel for PSR J0737 − 3039, the only suitable system existing until recently.
However, even a small misalignment angle would have significant consequences, because for a particular neutron star; the theoretical values of the moment of inertia using different EsoS are sometimes very close (depending on the stiffness of the EoS). One can see Figure 1 of Bejger et al. [53] for predicted values of the moment of inertia for PSR J0737 − 3039A using a sample of EsoS. The present article extensively studied the effects of such misalignment.
I have demonstrated the fact that the recently-discovered DNS, PSR J1757 − 1854, is almost as good as PSR J0737 − 3039A regarding the prospect of determining the moment of inertia of the pulsar from ω ˙ LT p . For PSR J0737 − 3039A, ω ˙ LT p lies in the range of 4.757 × 10 4 4.713 × 10 4 deg y 1 , whereas for PSR J1757 − 1854 , ω ˙ LT p lies in the range of 1.376 × 10 4 3.032 × 10 4 deg y 1 . Future measurement of χ p for PSR J1757 − 1854 will narrow down the range of ω ˙ LT p .
For both of the systems, significant improvements in the timing solution are needed. Moreover, the proper motion and the parallax measurements for PSR J1757 − 1854 are crucial. Additionally, it is essential to measure χ p and λ p for PSR J1757 − 1854; hence, long-term studies of polarization properties and profile variation would be very useful.
Note that, although all of the calculations presented in this article are within the framework devised by Barker & O’Connell [19] and Damour & Schäfer [23], alternative formalisms to incorporate the Lense-Thirring effect in ω ˙ p , obs exist. One example is Iorio [54], but for systems with i 90 , χ p 0 , these two formalisms take similar forms. Precisely, Iorio [54] reports ω ˙ LT p , = 3.7 × 10 4 deg y 1 for PSR J0737 − 3039A using I p = 1.00 × 10 45 gm cm 2 , which is close to what I get, i.e., ω ˙ LT p , = 3.766 × 10 4 deg y 1 if I use the same value of I p instead of 1.26 × 10 45 gm cm 2 used in this article. A fundamentally different idea, i.e., to incorporate the Lense-Thirring effect into the delays in the pulse arrival times, has been proposed recently [55]. If this can be implemented in future algorithms for pulsar timing, the effect of the Lense-Thirring effect could be measured directly.
Finally, if the measurement accuracy of ω ˙ p , obs improves by a few orders of magnitudes, one will need to subtract the Kopeikin term (the secular variation due to the proper motion). The measurement of the proper motion of PSR J1757 − 1854 will help us calculate this term. Moreover, in order to use P ˙ b to estimate masses accurately, one needs to subtract the contributions from the proper motion and the relative acceleration between the pulsar and the solar system barycenter. A package to perform these tasks based on the realistic potential of the Milky Way has been recently developed [30], and gradual improvement in the model of the galactic potential is expected with the help of Gaia data.

Acknowledgments

The author thanks the organizers of ‘Compact stars in the QCD phase diagram VI’ (especially David Blashcke) for their excellent hospitality, all three reviewers for extensive suggestions on the earlier version of the manuscript and ANI Technologies Pvt. Ltd. and Uber Technologies Inc. for helping to survive in Chennai.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Van Straten, W.; Bailes, M.; Britton, M.; Kulkarni, S.R.; Anderson, S.B.; Manchester, R.N.; Sarkissian, J. A test of general relativity from the three-dimensional orbital geometry of a binary pulsar. Nature 2001, 412, 158–160. [Google Scholar] [CrossRef] [PubMed]
  2. Deller, A.T.; Boyles, J.; Lorimer, D.R.; Kaspi, V.M.; McLaughlin, M.A.; Ransom, S.; Stairs, I.H.; Stovall, K. VLBI Astrometry of PSR J2222 − 0137: A Pulsar Distance Measured to 0.4% Accuracy. Astrophys. J. 2013, 770, 145. [Google Scholar] [CrossRef]
  3. Deller, A.T.; Vigeland, S.J.; Kaplan, D.L.; Goss, W.M.; Brisken, W.F.; Chatterjee, S.; Cordes, J.M.; Janssen, G.H.; Lazio, T.J.W.; Petrov, L.; et al. Microarcsecond VLBI Pulsar Astrometry with PSRπ. I. Two Binary Millisecond Pulsars with White Dwarf Companions. Astrophys. J. 2017, 828, 8. [Google Scholar] [CrossRef]
  4. Lorimer, D.; Kramer, M. Handbook of Pulsar Astronomy; Cambridge Observing Handbooks for Research Astronomers; Cambridge University Press: Cambridge, UK, 2005; pp. 32–58. [Google Scholar]
  5. Freire, P.C.C.; Wex, N. The orthometric parametrization of the Shapiro delay and an improved test of general relativity with binary pulsars. Mon. Not. R. Astron. Soc. 2010, 409, 199–212. [Google Scholar] [CrossRef]
  6. Kramer, M.; Stairs, I.H.; Manchester, R.N. Strong-field tests of gravity with the double pulsar. Annalen der Physik 2006, 15, 34–42. [Google Scholar] [CrossRef]
  7. Stairs, I. Testing General Relativity with Pulsar Timing. Living Rev. Relat. 2003, 6, 5. [Google Scholar] [CrossRef] [PubMed]
  8. Kramer, M. Pulsars as probes of gravity and fundamental physics. Int. J. Mod. Phys. 2016, 25, 1630029–1630061. [Google Scholar] [CrossRef]
  9. Kramer, M.; Stairs, I.H.; Manchester, R.N. Tests of General Relativity from Timing the Double Pulsar. Science 2006, 314, 97–102. [Google Scholar] [CrossRef] [PubMed]
  10. Demorest, P.B.; Pennucci, T.; Ransom, S.M.; Roberts, M.S.E.; Hessels, J.W.T. A two-solar-mass neutron star measured using Shapiro delay. Nature 2010, 467, 1081–1083. [Google Scholar] [CrossRef] [PubMed]
  11. Antoniadis, J.; Freire, P.C.C.; Wex, N.; Tauris, T.M.; Lynch, R.S. A Massive Pulsar in a Compact Relativistic Binary. Science 2013, 340, 1233232. [Google Scholar] [CrossRef] [PubMed]
  12. Abbott, B.P. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [PubMed]
  13. Annala, E.; Gorda, T.; Kurkela, A.; Vuorinen, A. Gravitational-wave constraints on the neutron-star-matter Equation of State. arXiv 2017, arXiv:1711.02644. [Google Scholar]
  14. Ayriyan, A.; Bastian, N.U.; Blaschke, D.; Grigorian, H.; Maslov, K.; Voskresensky, D.N. How robust is a third family of compact stars against pasta phase effects? arXiv 2017, arXiv:1711.03926. [Google Scholar]
  15. Kurkela, A.; Romatschke, P.; Vuorinen, A. Cold quark matter. Phys. Rev. D 2010, 81, 105021. [Google Scholar] [CrossRef]
  16. Perera, B.B.P.; McLaughlin, M.A.; Kramer, M.; Stairs, I.H.; Ferdman, R.D.; Freire, P.C.C.; Possenti, A.; Breton, R.P.; Manchester, R.N.; Burgay, M.; et al. The Evolution of PSR J0737 − 3039B and a Model for Relativistic Spin Precession. Astrophys. J. 2010, 721, 1193–1205. [Google Scholar] [CrossRef]
  17. Cameron, A.D.; Champion, D.J.; Kramer, M.; Bailes, M.; Barr, E.D. The High Time Resolution Universe Pulsar Survey—XIII. PSR J1757 − 1854, the most accelerated binary pulsar. 2018, 475, L57–L61. [Google Scholar] [CrossRef]
  18. Stovall, K.; Freire, P.C.C.; Chatterjee, S.; Demorest, P.B.; Lorimer, D.R. PALFA Discovery of a Highly Relativistic Double Neutron Star Binary. Astrophys. J. 2018, 231, 453. [Google Scholar]
  19. Barker, B.M.; O’Connell, R.F. The gravitational interaction: Spin, rotation, and quantum effects—A review. Gen. Relat. Gravit. 1979, 11, 149–175. [Google Scholar] [CrossRef]
  20. Perera, B.B.P.; Kim, C.; McLaughlin, M.A.; Ferdman, R.D.; Kramer, M.; Stairs, I.H.; Freire, P.C.C.; Possenti, A. Realistic Modeling of the Pulse Profile of PSR J0737 − 3039A. Astrophys J. 2014, 787, 51. [Google Scholar] [CrossRef]
  21. Pol, N.; McLaughlin, M.; Kramer, M.; Stairs, I.; Perera, B.B.P.; Possenti, A. A direct measurement of sense of rotation of PSR J0737 − 3039A. Astrophys. J. 2018, 853, 73. [Google Scholar] [CrossRef]
  22. Breton, R.P.; Kaspi, V.M.; Kramer, M. Relativistic Spin Precession in the Double Pulsar. Science 2008, 321, 104–107. [Google Scholar] [CrossRef] [PubMed]
  23. Damour, T.; Schafer, G. Higher-order relativistic periastron advances and binary pulsars. Nuovo Cim. B 1988, 101, 127–176. [Google Scholar] [CrossRef]
  24. Kopeikin, S.M. Proper Motion of Binary Pulsars as a Source of Secular Variations of Orbital Parameters. Astrophys. J. 1996, 467, L93. [Google Scholar] [CrossRef]
  25. Bagchi, M. Rotational parameters of strange stars in comparison with neutron stars. New Astron. 2010, 15, 126–134. [Google Scholar] [CrossRef]
  26. Ferdman, R.D.; Stairs, I.H.; Kramer, M.; Breton, R.P.; McLaughlin, M.A.; Freire, P.C.C.; Possenti, A.; Stappers, B.W.; Kaspi, V.M.; Manchester, R.N.; et al. The Double Pulsar: Evidence for Neutron Star Formation without an Iron Core-collapse Supernova. Astrophys. J. 2013, 767, 85. [Google Scholar] [CrossRef]
  27. Akmal, A.; Pandharipande, V.R.; Ravenhall, D.G. Equation of state of nucleon matter and neutron star structure. Phys. Rev. C 1998, 58, 1804–1828. [Google Scholar] [CrossRef]
  28. Nozawa, T.; Stergioulas, N.; Gourgoulhon, E.; Eriguchi, Y. Construction of highly accurate models of rotating neutron stars—Comparison of three different numerical schemes. Astron. Astrophys. Suppl. Ser. 1998, 132, 431–454. [Google Scholar] [CrossRef]
  29. Deller, A.T.; Bailes, M.; Tingay, S.J. Implications of a VLBI Distance to the Double Pulsar J0737-3039A/B. Science 2009, 323, 1327–1329. [Google Scholar] [CrossRef] [PubMed]
  30. Pathak, D.; Bagchi, M. GalDynPsr: A package to estimate dynamical contributions in the rate of change of the period of radio pulsars. arXiv 2017, arXiv:1712.06590. [Google Scholar]
  31. Kramer, M.; Wex, N. TOPICAL REVIEW: The double pulsar system: A unique laboratory for gravity. Class. Quantum Gravity 2009, 26, 073001. [Google Scholar] [CrossRef]
  32. Kehl, M.S.; Wex, N.; Kramer, M.; Liu, K. Future measurements of the Lense-Thirring effect in the Double Pulsar. In Proceedings of the Fourteenth Marcel Grossmann Meeting, Rome, Italy, 12–18 July 2015. [Google Scholar]
  33. Cordes, J.M.; Lazio, T.J.W. NE2001.I. A New Model for the Galactic Distribution of Free Electrons and its Fluctuations. arXiv 2001, arXiv:astro-ph/0207156. [Google Scholar]
  34. Cordes, J.M.; Lazio, T.J.W. NE2001. II. Using Radio Propagation Data to Construct a Model for the Galactic Distribution of Free Electrons. arXiv 2003, arXiv:astro-ph/0301598. [Google Scholar]
  35. Yao, J.M.; Manchester, R.N.; Wang, N. A New Electron-density Model for Estimation of Pulsar and FRB Distances. Astrophys. J. 2017, 835, 29. [Google Scholar] [CrossRef]
  36. Weisberg, J.M.; Huang, Y. Relativistic measurements from timing the binary pulsar PSR B1913 + 16. Astrophys. J. 2016, 829, 55. [Google Scholar] [CrossRef]
  37. Martinez, J.G.; Stovall, K.; Freire, P.C.C.; Deneva, J.S.; Jenet, F.A.; McLaughlin, M.A.; Bagchi, M.; Bates, S.D.; Ridolfi, A. Pulsar J0453+1559: A Double Neutron Star System with a Large Mass Asymmetry. Astrophys. J. 2015, 812, 143. [Google Scholar] [CrossRef]
  38. Martinez, J.G.; Stovall, K.; Freire, P.C.C.; Deneva, J.S.; Tauris, T.M.; Ridolfi, A.; Wex, N.; Jenet, F.A.; McLaughlin, M.A.; Bagchi, M. Pulsar J1411+2551: A Low-mass Double Neutron Star System. Astrophys. J. 2017, 851, L29. [Google Scholar] [CrossRef]
  39. Janssen, G.H.; Stappers, B.W.; Kramer, M.; Nice, D.J.; Jessner, A.; Cognard, I.M.B.; Purver, M.B. Multi-telescope timing of PSR J1518 + 4904. Astron. Astrophys. 2008, 490, 753–761. [Google Scholar]
  40. Fonseca, E.; Stairs, I.H.; Thorsett, S.E. A Comprehensive Study of Relativistic Gravity Using PSR B1534 + 12. Astrophys. J. 2014, 787, 82. [Google Scholar] [CrossRef]
  41. Keith, M.J.; Kramer, M.; Lyne, A.G.; Eatough, R.P.; Stairs, I.H.; Possenti, A.; Camilo, F.; Manchester, R.N. PSR J1753 − 2240: A mildly recycled pulsar in an eccentric binary system. Mon. Not. R. Astron. Soc. 2009, 393, 623–627. [Google Scholar] [CrossRef]
  42. Ferdman, R.D.; Stairs, I.H.; Kramer, M.; Janssen, G.H.; Bassa, C.G.; Stappers, B.W.; Demorest, P.B.; Cognard, I.; Desvignes, G.; Theureau, G.; et al. PSR J1756 − 2251: A pulsar with a low-mass neutron star companion. Mon. Not. R. Astron. Soc. 2014, 443, 2183–2196. [Google Scholar] [CrossRef] [Green Version]
  43. Lynch, R.S.; Freire, P.C.C.; Ransom, S.M.; Jacoby, B.A. The Timing of Nine Globular Cluster Pulsars. Astrophys. J. 2012, 745, 109. [Google Scholar] [CrossRef]
  44. Corongiu, A.; Kramer, M.; Stappers, B.W.; Lyne, A.G.; Jessner, A.; Possenti, A.; D’Amico, N.; Löhmer, O. The binary pulsar PSR J1811 − 1736: Evidence of a low amplitude supernova kick. Astron. Astrophys. 2007, 462, 703–709. [Google Scholar] [CrossRef]
  45. Hobbs, G.; Lyne, A.G.; Kramer, M.; Martin, C.E.; Jordan, C. Long-term timing observations of 374 pulsars. Mon. Not. R. Astron. Soc. 2004, 353, 1311–1344. [Google Scholar] [CrossRef]
  46. Champion, D.J.; Lorimer, D.R.; McLaughlin, M.A.; Cordes, J.M.; Arzoumanian, Z.; Weisberg, J.M.; Taylor, J.H. PSR J1829 + 2456: A relativistic binary pulsar. Mon. Not. R. Astron. Soc. 2004, 350, L61–L65. [Google Scholar] [CrossRef]
  47. Van Leeuwen, J.; Kasian, L.; Stairs, I.H.; Lorimer, D.R.; Camilo, F. The Binary Companion of Young, Relativistic Pulsar J1906+0746. Astrophys. J. 2015, 798, 118. [Google Scholar] [CrossRef]
  48. Weisberg, J.M.; Nice, D.J.; Taylor, J.H. Timing Measurements of the Relativistic Binary Pulsar PSR B1913 + 16. Astrophys. J. 2010, 722, 1030–1034. [Google Scholar] [CrossRef]
  49. Swiggum, J.K.; Rosen, R.; McLaughlin, M.A.; Lorimer, D.R.; Heatherly, S. PSR J1930 − 1852: A Pulsar in the Widest Known Orbit around Another Neutron Star. Astrophys. J. 2015, 805, 156. [Google Scholar] [CrossRef]
  50. Jacoby, B.A.; Cameron, P.B.; Jenet, F.A.; Anderson, S.B.; Murty, R.N.; Kulkarni, S.R. Measurement of Orbital Decay in the Double Neutron Star Binary PSR B2127 + 11C. Astrophys. J. 2006, 644, L113–L116. [Google Scholar] [CrossRef]
  51. Iorio, L. Prospects for measuring the moment of inertia of pulsar J0737-3039A. New Astron. 2009, 1, 40–43. [Google Scholar] [CrossRef]
  52. O’Connell, R.F. Proposed New Test of Spin Effects in General Relativity. Phys. Rev. Lett. 2004, 93, 081103. [Google Scholar] [CrossRef] [PubMed]
  53. Bejger, M.; Bulik, T.; Haensel, P. Constraints on the dense matter equation of state from the measurements of PSR J0737 − 3039A moment of inertia and PSR J0751 + 1807 mass. Mon. Not. R. Astron. Soc. 2005, 364, 635–639. [Google Scholar] [CrossRef]
  54. Iorio, L. General relativistic spin-orbit and spin-spin effects on the motion of rotating particles in an external gravitational field. Gen. Relat. Gravit. 2012, 44, 719–736. [Google Scholar] [CrossRef]
  55. Iorio, L. Post-Keplerian perturbations of the orbital time shift in binary pulsars: an analytical formulation with applications to the galactic center. Eur. Phys. J. C 2017, 77, 439. [Google Scholar] [CrossRef]
1.
The remaining Keplerian parameter, the longitude of the ascending node φ , does not come in the standard pulsar timing algorithm. It can be measured via proper motion only in very special cases; see [1,2,3] for details.
2.
Although for the last few years, the slow pulsar of the system is not visible and is believed to be beaming away from the Earth due to its spin-precession [16].
3.
RNS stands for ‘Rapidly Rotating Neutron Star’, a package to calculate different properties of rotating neutron stars, freely available at http://www.gravity.phys.uwm.edu/rns/.
Figure 1. Orientation of different vectors relevant for the estimation of ω ˙ LT a where the subscript a refers to the object being observed, i.e., the pulsar (so the companion would be represented by the subscript a + 1 ). The vectors are as follow: k is the unit vector along the orbital angular momentum, h a is the unit vector along the line-of-sight and s a is the unit spin vector. The angles in the figures are: i is the inclination angle between the orbit of the pulsar and the sky-plane, as well as the angle between k and h a , χ a is the angle between k and s a and λ a is the angle between h a and s a . The left and right panels are for χ a < i and χ a > i , respectively. In both panels, two positions of s a are shown corresponding to maximum and minimum values of λ a (see the text). The pulsar, i.e., the object a, is located at the vertex of the fiducial cone (gray in color) made by the precessing s a around k .
Figure 1. Orientation of different vectors relevant for the estimation of ω ˙ LT a where the subscript a refers to the object being observed, i.e., the pulsar (so the companion would be represented by the subscript a + 1 ). The vectors are as follow: k is the unit vector along the orbital angular momentum, h a is the unit vector along the line-of-sight and s a is the unit spin vector. The angles in the figures are: i is the inclination angle between the orbit of the pulsar and the sky-plane, as well as the angle between k and h a , χ a is the angle between k and s a and λ a is the angle between h a and s a . The left and right panels are for χ a < i and χ a > i , respectively. In both panels, two positions of s a are shown corresponding to maximum and minimum values of λ a (see the text). The pulsar, i.e., the object a, is located at the vertex of the fiducial cone (gray in color) made by the precessing s a around k .
Universe 04 00036 g001
Figure 2. Variation of ω ˙ LT p with χ p (in degrees) along the horizontal axis and λ p (in degrees) along the vertical axis for PSR J1757 − 1854. The color code represents values of ω ˙ LT p in deg y 1 .
Figure 2. Variation of ω ˙ LT p with χ p (in degrees) along the horizontal axis and λ p (in degrees) along the vertical axis for PSR J1757 − 1854. The color code represents values of ω ˙ LT p in deg y 1 .
Universe 04 00036 g002
Table 1. Relevant observed parameters for known pulsars in DNSs. See the references for these parameters with more significant digits and uncertainties. The timing model (e.g., DDGR) is mentioned when more than one timing solution is available in the original reference.
Table 1. Relevant observed parameters for known pulsars in DNSs. See the references for these parameters with more significant digits and uncertainties. The timing model (e.g., DDGR) is mentioned when more than one timing solution is available in the original reference.
DNSMassese P b sin i P s 1 ω ˙ χ p Refs.
M Daysms deg y 1 deg
J0453+15591.5590.112518324.0724685880.967134345.7818161630930.0379412R1
(DDGR)1.174
J0737−3039A 1.3381 0.08777750.102251562480.9997422.7016.89947 < 6.1 or < 2.3 R2, R χ 1
J0737−3039B 1.2489 2773.46138
J1411+2551 1.62 0.16993082.6158567793962.4528955175900.0768R3
0.92
J1518+4904 1.17 0.249484518.6340050964 0.73 40.9349889080.0113725R4
1.55
B1534+12 1.3455 0.273677400.4207372988810.9749637.90444117831.755795 27.0 ± 3.0 R5
(DDGR) 1.3330 153.0 ± 3.0
J1753−2240 0.30358213.637566895.1378086771R6
0.4875
J1756−22511.3410.18056940.319633901430.9328.46158902599832.58240 < 5.9 R7
1.230( 1 σ value)
J1757−18541.33840.60581420.183537835870.994521.49723189002710.3651 25 R8
1.3946(theoretical)
J1807−2500B 1.36550.7470331989.95666815880.99564.186177202840890.0183389R9
( g 1 ) 1.2064
J1811−1736 1.64 0.82801118.7791691104.18195479680.0090R10
0.93
B1820−11 0.794608357.76199279.8286965653490.00007R11
J1829+2456 1.38 0.139141.17602841.009823580.28R12
1.22
J1906+07461.2910.08529960.16599304686 0.690882411144.073155387.5844 110 55 + 21 R13
(DDGR)1.322
B1913+16 1.4398 0.61713340.3229974489110.68 59.030003224.22659822R14, Rχ1
1.3886
J1930−1852 1.32 0.3988634045.0600007185.520160479260.00078R15
1.30
J1946+20520.0638480.0784880416.960175323025.6R16
B2127+11C1.3580.6813950.3352820482830.529296148644.4644R17
( g 2 ) 1.354
Candidate DNS; was never published, but here, I quote the value used byWeisberg & Huang [36]; pulsars denoted with g are in globular clusters; g1 is in NGC 6544, and g2 is in M15; References: R1: Martinez et al. [37], R2: Kramer et al. [9], R3: Martinez et al. [38], R4: Janssen et al. [39], R5: Fonseca, Stairs, Thorsett [40], R6: Keith et al. [41], R7: Ferdman et al. [42], R8: Cameron et al. [17], R9: Lynch et al. [43], R10 Corongiu et al. [44], R11: Hobbs et al. [45], R12: Champion et al. [46], R13: van Leeuwen et al. [47], R14: Weisberg, Nice, & Taylor [48], R15: Swiggum et al. [49], R16: Stovall et al. [18], R17: Jacoby et al. [50], Rχ1 (for χp): Perera et al. [20].
Table 2. Values of ω ˙ PN p , ω ˙ LT p , and ω ˙ LT p , m a x for known pulsars in DNSs using I p = 1.26 × 10 45 gm cm 2 and neglecting the Lense-Thirring effect of the companion. To show the difference between ω ˙ LT p , and ω ˙ LT p , m a x , I keep values up to five decimal places even for the cases where sin i is known with less accuracy.
Table 2. Values of ω ˙ PN p , ω ˙ LT p , and ω ˙ LT p , m a x for known pulsars in DNSs using I p = 1.26 × 10 45 gm cm 2 and neglecting the Lense-Thirring effect of the companion. To show the difference between ω ˙ LT p , and ω ˙ LT p , m a x , I keep values up to five decimal places even for the cases where sin i is known with less accuracy.
DNS ω ˙ PN p ω ˙ LT p , ω ˙ LT p , max
( deg y 1 ) ( deg y 1 ) ( deg y 1 )
J0453+15590.03794 1.30190 × 10 7 1.31310 × 10 7
J0737−3039A16.90312 4.74458 × 10 4 4.74489 × 10 4
J1411+25510.07681 2.32510 × 10 7
J1518+49040.01138 4.48005 × 10 8
B1534+121.75533 1.84880 × 10 5 1.86078 × 10 5
J1753−2240
J1756−22512.58363 4.01930 × 10 5 4.09703 × 10 5
J1757−185410.36772 3.02752 × 10 4 3.03170 × 10 4
J1807−2500B0.01834 8.98082 × 10 7 8.99076 × 10 7
J1811−17360.00895 1.45038 × 10 8
B1820−11
J1829+24560.29284 1.96704 × 10 6
J1906+07467.58528 2.91884 × 10 5 3.29423 × 10 5
B1913+164.22760 3.43984 × 10 5 3.90790 × 10 5
J1930−18520.00079 3.86975 × 10 10
J1946+2052
B2127+11C4.46458 8.11241 × 10 5

Share and Cite

MDPI and ACS Style

Bagchi, M. Prospects of Constraining the Dense Matter Equation of State from Timing Analysis of Pulsars in Double Neutron Star Binaries: The Cases of PSR J0737 ‒ 3039A and PSR J1757 ‒ 1854. Universe 2018, 4, 36. https://doi.org/10.3390/universe4020036

AMA Style

Bagchi M. Prospects of Constraining the Dense Matter Equation of State from Timing Analysis of Pulsars in Double Neutron Star Binaries: The Cases of PSR J0737 ‒ 3039A and PSR J1757 ‒ 1854. Universe. 2018; 4(2):36. https://doi.org/10.3390/universe4020036

Chicago/Turabian Style

Bagchi, Manjari. 2018. "Prospects of Constraining the Dense Matter Equation of State from Timing Analysis of Pulsars in Double Neutron Star Binaries: The Cases of PSR J0737 ‒ 3039A and PSR J1757 ‒ 1854" Universe 4, no. 2: 36. https://doi.org/10.3390/universe4020036

APA Style

Bagchi, M. (2018). Prospects of Constraining the Dense Matter Equation of State from Timing Analysis of Pulsars in Double Neutron Star Binaries: The Cases of PSR J0737 ‒ 3039A and PSR J1757 ‒ 1854. Universe, 4(2), 36. https://doi.org/10.3390/universe4020036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop