Pulse Profiles with NICER

# Measuring Neutron Star Radii via Pulse Profile Modeling with NICER

Feryal Özel11affiliation: Astronomy Department, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721, USA; fozel@email.arizona.edu , Dimitrios Psaltis11affiliation: Astronomy Department, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721, USA; fozel@email.arizona.edu , Zaven Arzoumanian22affiliation: Center for Research and Exploration in Space Science and Technology & X-Ray Astrophysics Laboratory, NASA Goddard Space Flight Center, Code 662, Greenbelt, MD 20771, USA , Sharon Morsink33affiliation: Department of Physics, University of Alberta, 11455 Saskatchewan Drive, Edmonton, AB T6G 2E1, Canada , and Michi Bauböck11affiliation: Astronomy Department, University of Arizona, 933 N. Cherry Ave., Tucson, AZ 85721, USA; fozel@email.arizona.edu
###### Abstract

The Neutron-star Interior Composition Explorer (NICER) is an X-ray astrophysics payload that will be placed on the International Space Station. Its primary science goal is to measure with high accuracy the pulse profiles that arise from the non-uniform thermal surface emission of rotation-powered pulsars. Modeling general relativistic effects on the profiles will lead to measuring the radii of these neutron stars and to constraining their equation of state. Achieving this goal will depend, among other things, on accurate knowledge of the source, sky, and instrument backgrounds. We use here simple analytic estimates to quantify the level at which these backgrounds need to be known in order for the upcoming measurements to provide significant constraints on the properties of neutron stars. We show that, even in the minimal-information scenario, knowledge of the background at a few percent level for a background-to-source countrate ratio of 0.2 allows for a measurement of the neutron star compactness to better than 10% uncertainty for most of the parameter space. These constraints improve further when more realistic assumptions are made about the neutron star emission and spin, and when additional information about the source itself, such as its mass or distance, are incorporated.

stars: neutron — dense matter – equation of state – pulsars: general — gravitation

## 1. Introduction

Periodic brightness oscillations originating from the surface of a neutron star provide a powerful means of measuring its mass and radius (see Strohmayer 2004 for an early review). In the case of rotation-powered pulsars, the temperature non-uniformity that causes the brightness oscillations is due to the magnetic field topology of the stellar surface and the associated energetic particle flux from the magnetosphere, which heats the polar caps relative to the rest of the star such that they emit in the soft X-ray band (e.g., Ruderman & Sutherland 1975; Arons 1981; Harding & Muslimov 2001, 2002; Bai & Spitkovsky 2010a, b). The trajectories of the photons emitted from the surface are bent by the strong gravitational field of the neutron star, imparting on the pulse shapes and amplitudes information about the stellar compactness (e.g., Pechenick et al. 1983; Weinberg et al. 2001; Poutanen & Beloborodov 2006; Cadeau et al. 2007; Psaltis & Özel 2014).

The degree to which this information can be decoded to infer the mass and radius, either jointly or separately, depends on a number of factors such as the X-ray brightness of the neutron star, its spin period, , the beaming pattern of the radiation emitted from the surface, as well as the spot colatitude and the inclination of the observer with respect to the stellar spin axis. This is because all of these variables, along with the stellar compactness, determine the pulse amplitude and its shape (i.e., the amplitudes of the harmonics). In addition, there are several sources of background that affect the measurement of the properties of the pulse profiles, such as X-ray emission from the magnetosphere or the rest of the neutron star surface, X-ray emission from other sources in the field of view, and instrumental background. The resulting signal-to-background ratio ultimately determines the accuracy with which the stellar compactness, , can be measured for a star of mass and radius .

NICER is a NASA mission, currently in development, designed primarily to observe thermal X-rays originating from the polar caps of a number of rotation-powered millisecond pulsars. Its goal is to analyze pulse profiles to measure the stellar compactness and radius (Gendreau et al. 2012; Arzoumanian et al. 2014). Being a dedicated mission, it will have the ability to observe selected sources for long periods to accumulate the necessary number of counts to measure the properties of the pulse profile. Nevertheless, the accuracy it will achieve in the measurement of the various harmonic amplitudes, and, thus, in the compactness, will depend on the level and our knowledge of the several sources of background. Because it is not an imaging instrument and because of the International Space Station environment, it will register both time-varying and constant sources of background, arising from the various origins discussed above. In this paper, we formulate the effect of the background on the pulse profiles and investigate the background tolerances that will still enable 5-10% accuracy in the measurements of the neutron star properties, which is necessary to pin down the dense matter equation of state. We provide approximate formulae for the uncertainties in the inferred stellar compactness and radius as a function of the signal-to-noise ratio, the uncertainty in the measurement of the background, and the uncertainty in the system geometry, described by the spot colatitude and the observer’s inclination.

Thermal emission from the polar caps of rotation powered pulsars has been observed with previous generations of X-ray satellites (Ögelman 1991; Finley et al. 1992; Becker & Trümper 1993; Pavlov & Zavlin 1997). Pulsed X-ray emission has been detected and the thermal X-ray brightness has been determined for a number of sources, such as PSR J04374715, PSR J00300451, PSR J21243358, and PSR J10240719, most recently by targeted Chandra and XMM-Newton observations (e.g., Zavlin 2006; Bogdanov et al. 2007; Bogdanov & Grindlay 2009; Bogdanov 2013). In addition, the number of rotation-powered millisecond pulsars has grown exponentially in recent years, thanks to discoveries by the Fermi satellite (e.g., Abdo et al. 2013) and radio surveys1, some of which are likely to exhibit pulsed soft X-ray emission. The spectra are usually described by a thermal component originating from the surface, typically modeled by unmagnetized hydrogen atmospheres or blackbodies, as well as a non-thermal component, likely of magnetospheric origin, which manifests itself as a power-law tail at higher energies. Using this pulsed surface emission to constrain the neutron star properties has so far led to broad constraints on masses and radii (Zavlin 2006; Bogdanov et al. 2007; Bogdanov & Grindlay 2009; Bogdanov 2013).

The rotation-powered pulsars that are the primary targets for NICER have spin frequencies of  Hz. This range of relatively low frequencies (compared to neutron stars in some accreting systems) minimizes higher order spin effects such as the stellar oblateness and the spacetime quadrupole (Morsink et al. 2007; Cadeau et al. 2007; Psaltis & Özel 2014; Psaltis et al. 2014). The first order Doppler effects can still be appreciable in this range of spin frequencies. This simplifies the modeling by allowing us to treat the stellar spacetime in the Schwarzschild+Doppler approximation. A large amount of work has been carried out in this regime, starting with the numerical calculations of self-lensing in the Schwarzschild metric by Pechenick et al. (1983) and its application to the millisecond brightness oscillations observed by the Rossi X-ray Timing Explorer (Strohmayer et al. 1997; Nath et al. 2002). Furthermore, Poutanen & Beloborodov (2006) devised an analytic expression that we use in the present study when calculating the amplitudes of the harmonics.

In this paper, we investigate the effect of background on the determination of the neutron star compactness given measurements of the shapes of pulse profiles from rotation-powered pulsars. Specifically, we make use of the analytic formulae for the pulse shapes of moderately spinning neutron stars in the Schwarzschild metric to calculate the number of observables one can determine from the profiles in different situations. We then use these to quantify the level at which the various sources of background must be independently known in order for the upcoming NICER measurements to provide significant constraints on the properties of neutron stars.

## 2. Measuring the Properties of Pulse Profiles

In the following, we will assume that the surface emission of a slowly spinning, rotation powered pulsar originates in one small circular polar cap on its surface. Bauböck et al. (2015) estimated the typical sizes of polar caps, defined as the radius within which the open magnetic field lines originate from the neutron star surface, for the NICER targets. The analytic estimates in that study of the size of the polar cap (Bauböck et al. 2015), as well as numerical simulations of surface heating caused by return magnetospheric currents (e.g., Bai & Spitkovsky 2010), support the assumption that the polar cap size is too small for its detailed emission structure to affect the pulse profile observed at infinity. Because of this, we make the small spot approximation when calculating the harmonic amplitudes of the pulse profiles. However, observations of the prototypical NICER target, PSR J04374715, suggest that its surface emission may arise from two small but non-antipodal polar caps (Bogdanov 2014). This latter case is more complicated than what we explore here and will change the actual shape of the pulse profile, but its complexity does not severely affect our conclusions in the uncertainties achieved in the various quantities as a function of count rate and backgrounds.

For a slowly spinning neutron star, the pulse profile observed by NICER will be characterized by the following quantities (see Figure 1): the background countrate , the constant component of the pulse profile , the amplitudes and the phases of the various harmonics, and, if present, the phase beyond which the flux from the polar cap is blocked by the star. The general form of the profile will be

 F(t)=A+{Q+∑n=1Vncos(2πntP+ϕn),if |2πt/P|<ϕp0,otherwise (1)

If the spin of the neutron star is significant, then additional harmonics appear in the profile with orthogonal phases.

The total number of counts collected when periods are accumulated is equal to

 C = N∫P0F(t)dt (2) = (A+Qϕpπ)NP+2N∑n=1Vn∫Pϕp/2π0cos(2πntP+ϕn)dt = (A+Qϕpπ)NP+∑n=1VnNPπ[sin(nϕp+ϕn)−sinϕn].

The uncertainty in the measurement of the total number of counts will simply be

 σC=√C≃√(A+Q)NP, (3)

where, in the last equality, we assumed that the polar cap is visible throughout the pulse period (i.e., ). From this last equation, we also get the uncertainty in measuring the sum as

 σQ+A=σcNP=√A+QNP. (4)

In the following calculation, we will need primarily the constant component of the source countrate and the associated uncertainty . In order to do this, we will need a measurement of the background countrate, , and its associated uncertainty, . The latter will have a Poisson component that will depend on the total time used to measure the background, as well as a possible systematic component, such as

 σA=√AT+λA (5)

with an appropriate constant. With these definitions, we obtain

 σQQ = (σ2Q+A+σ2A)1/2Q (6) = [(AQ+1)1QNP+(σAA)2(AQ)2]1/2.

We can also write this equation in terms of the total source counts and the total background counts as

 σQQ=[(BS+1)1S+(σBB)2(BS)2]1/2. (7)

Figure 2 shows the predicted fractional uncertainties in the measurement of the quantity as a function of the total number of source counts collected. In making this figure, we have assumed a constant fractional uncertainty in the a priori knowledge of the background, i.e., , and a constant ratio of background to source counts, i.e., . Under these assumptions, the asymptotic value of the uncertainty (i.e., for a large number of source counts) is dominated by the uncertainty in the knowledge of the background and scales as

 limS→∞σQQ=(σBB)(BS). (8)

In other words, the fractional uncertainty in the measurement of is limited by the product of the fractional uncertainty in the knowledge of the background times the ratio of the background to source counts.

The additional quantities that we can measure from the pulse profile are the fractional r.m.s. amplitudes of the various harmonics. Considering, for example, the th harmonic, the measured fractional r.m.s. amplitude (including the background) is

 robs=Vn√2(Q+A). (9)

Using the result in Psaltis et al. (2014), the uncertainty in the measurement of scales as

 σr,obs=√S+BS=√(A+Q)NPQNP, (10)

from which we can write

 σ2V = (∂Vn∂robs)2σ2r,obs+[∂Vn∂(Q+A)]2σ2Q+A (11) = 2(Q+A)3Q2NP+V2n(Q+A)NP.

The fractional uncertainty in the measurement of the amplitude is

 σVVn=[2(1+AQ)3(QVn)21QNP+(1+AQ)−11QNP]1/2. (12)

Writing this equation in terms of the total number of source and background counts, we obtain

 σVVn=[2(1+BS)3(QVn)2+(1+BS)−1]1/21S1/2. (13)

Figure 3 shows the predicted fractional uncertainty in the measurement of the quantity , as a function of the total number of source counts. As the last equation shows, the uncertainty in the measured amplitude of pulsations simply scales as the inverse square root of the source counts.

Finally, the last quantity to be measured, when present in the pulse profile, is the phase , beyond which the polar cap is blocked by the star. The main uncertainty in this measurement will arise from the fact that the evolution of the pulse profile with phase is masked by the uncertainty in the background, i.e.,

 V1cos(ϕP+δϕP)≃σA. (14)

Setting , we can use this equation to find

 σϕ=σAV1sinϕp. (15)

and write this last equation in the more useful form

 (16)

For the remainder of our discussion, we will only consider pulse shapes in which the polar cap is visible throughout the pulse period.

## 3. From pulse profiles to stellar compactness

Having measured the characteristic properties of a pulse profile, i.e., its constant component and the amplitudes and phases of its various harmonics, we can use theoretical models of the gravitationally lensed emission from a neutron-star hotspot to obtain the compactness and the radius of the neutron star.

In the following, we will make use of the approximate analytic expressions and terminology of Poutanen & Beloborodov (2006). These expressions take into account light bending around a slowly spinning star and yield the bolometric flux observed at infinity under the assumption that the polar cap is infinitesimally small and that no emission emerges from the rest of the star. We will denote by and the mass and radius of the neutron star, by its spin period, by the inclination of the observer with respect to the spin axis, and by the colatitude of the hot spot.

Model atmospheres of neutron stars, in general, and of magnetic polar caps, in particular, show that the emission from the neutron-star surface is not isotropic but is beamed (see, e.g., the discussion in Özel 2013). We will represent this beaming with a functional form

 I(α)=I0(1+hcosα), (17)

where is the emission angle with respect to the normal to the surface and is a quantity that is calculated from model atmospheres and measures the degree of beaming ( being the isotropic limit).

We will also define the compactness of the neutron star as

 u≡2GMRc2 (18)

as well as the following two quantities

 q ≡ u+(1−u)cosicosθ v ≡ (1−u)sinisinθ. (19)

The relations between the characteristics of the pulse profile (see eq. [1]) and the properties of a slowly spinning neutron star are (Poutanen & Beloborodov 2006, eq. [28]-[38])

 Q = I0GdSD2[q+h(q2+v22)] V1 = I0GdSD2(1+2hq)vcosϕ1 V2 = I0GdSD2hv22cosϕ2 tanϕ1 = (3+Γ)q+(4+Γ)h(q2+v2/4)u(1+2hq)(2πRcP)sinisinθ tanϕ2 = (4+Γ)(1+2hq)−1hu(2πRcP)sinisinθ. (20)

In these equations, is the surface area of the polar cap, is the distance to the source, is the photon index of the spectrum, and is a (calculable) factor that depends on the spectrum of the source and the photon-energy range of the detector.

### 3.1. The Limiting Case of Minimal Information

To calculate the background tolerances in the measurements of the neutron star properties from the characteristics of its pulse profile, we will first consider the worst possible case, i.e., one in which the emission from the neutron star surface is isotropic and the Doppler effects on the pulse profile are not measurable. Although this is not expected in reality, this situation both reduces the amplitude of the fundamental (exacerbating background effects) and the amplitudes of the harmonics, thus minimizing the number of measurable quantities in the profile. This can be seen by taking the limit and in the above equations, which simplifies them to

 Q = I0GdSD2q V1 = I0GdSD2v V2 = 0. (21)

One of the most significant benefits of neutron-star radius measurements with pulse profile modeling is that, in principle, they can provide results that are independent of any prior knowledge of the distance to each source (although, as we will discuss below, when precise distance measurements are available, they can be folded in to increase the number of independently measured quantities from the pulse profile). In order not to rely on the distance information, we cannot separately use the absolute constant component and fundamental amplitude of the pulse profile (both of which depend on the size of the polar cap and on the distance), but we can only consider the fractional amplitude of the fundamental, i.e.,

 r1≡V1Q=vq=(1−u)sinisinθu+(1−u)cosicosθ. (22)

In this last step, we used the definitions (19) to connect the pulse profile properties to the neutron-star compactness. Solving this last equation for the compactness , we obtain

 u=r1cosicosθ−sinisinθ−r1+r1cosicosθ−sinisinθ. (23)

Clearly, in order to measure the neutron-star compactness in this case, we need independent knowledge of the inclination of the observer, , and of the colatitude of the polar cap, . Such prior knowledge is not unreasonable for many rotation-powered pulsars, the geometries of which can be inferred by fitting model -ray pulse profiles to Fermi observations of individual sources (Johnson et al. 2014) and, to a certain extent, by modeling polarization angle swings in the radio (Yan et al. 2011; Dai et al. 2015).

The uncertainty within which the fractional amplitude of the fundamental can be measured is

 σr1r1=(σ2V1V21+σ2QQ2)1/2, (24)

where and are given by equations (7) and (13), respectively. At the limit of a very large number of counts,

 limS→∞σr1r1=(σBB)(BS). (25)

Using these expressions, we can finally write the uncertainty in the measurement of the neutron-star compactness as

 σu = [(∂u∂r1)2σ2r1+(∂u∂i)2σ2i+(∂u∂θ)2σ2θ]1/2 (26) = [σ2r1sin2isin2θ+12r31sin(2i)sin(2θ)(σ2i+σ2θ) +r21cos2isin2θ(σ2i+r21σ2θ) +r21sin2icos2θ(r21σ2i+σ2θ)]1/2 (r1−r1cosicosθ+sinisinθ)−2.

Here, and are the uncertainties (in radians) in the prior knowledge of the observer inclination and of the colatitude of the polar cap.

Figure 4 shows the predicted fractional uncertainties in the measurement of the compactness of a neutron star, for different parameters of the system as well as for different degrees of knowledge of the background to the source. When the geometric angles and are reasonably well known via other measurements, these figures show that knowledge of the background at a few percent level allows a measurement of the neutron star compactness to better than uncertainty for most of the parameter space. This result is obtained for a background-to-source ratio of 0.2, which is achievable for the brightest NICER targets, and a total accumulation of counts. It is only when the pulse amplitudes plummet for the cases where either the observer or the polar cap are nearly aligned with the spin axis that the uncertainties become larger than (see the left panel of Fig. 4). When the system geometry is such that the pulse amplitudes are reasonably large (as is the case with prototypical NICER targets), the right panel of Figure 4 shows the dependence of the fractional error in the neutron star compactness for different background-to-signal ratios and different levels of uncertainty in the determination of the background. Even when the background is comparable to the signal and is measured only to within , the fractional error in the compactness remains at or below . Note that this is for the minimal-information case of isotropic beaming, where only one fractional amplitude can be obtained from the pulse profiles.

We also calculate the limiting case of negligible uncertainties in the knowledge of these angles and of an infinite number of source counts. The uncertainty in the measurement of the neutron-star compactness then becomes

 limS→∞σu=(σBB)(BS)|sinisinθ|r1(−r1cosicosθ+sinisinθ+r1)2. (27)

Using equation (22) and after some rearranging of terms, we reach the final result for the limiting case

 (28)

This expression demonstrates that, even in the case of perfect prior knowledge of the geometry of the system and of an infinite number of accumulated photons, the accuracy in the fractional measurement of the neutron star compactness is proportional to the uncertainty in our measurement of the background () and to the ratio of the background-to-source counts ().

In the previous section, we considered the case of minimal information that can be obtained from the pulse profiles, assuming isotropic emission, slow rotation, and no distance or mass information for the pulsar. We now discuss the additional constraints that can be achieved on the neutron star radius and compactness, if any of these assumptions is relaxed.

We first explore the more favorable case when the emission is anisotropic, i.e., . This is also a more realistic consideration because of the strong temperature gradients in the bombarded atmospheres of the polar cap regions. The uncertainties discussed in the previous section all improve because of the presence of one additional constraint related to the measurement of the second harmonic of the profile. The two independently measurable quantities then become

 V1Q = (1+2hq)vq+h(q2+v2/2) V2Q = hv22[q+h(q2+v2/2)]. (29)

These two independent pieces of information allow for a unique determination of the auxiliary quantities and , which can be used to determine the compactness and the combination of the angles and (since these angles enter and only as a product of their sines and cosines, see eq.[19]). The fractional uncertainties and , then, will give rise to similar uncertainties in the two ratios in equation (29). Individually, these translate into uncertainties in the compactness in much the same way as in the previous section, following Equations (22)-(26). However, combining them allows for smaller uncertainties on the compactness and the combination of the angles and . More harmonics can be measured also in the case of a moderately spinning star where Doppler effects are detectable (e.g., Psaltis & Özel 2014), leading to similar tightening of the uncertainties.

A similar situation arises when the polar-cap size and the measured intensity of radiation are measured spectroscopically and the distance to the source is known a priori. This is achievable for some of the representative NICER targets where a precise measurement of the distance is possible via pulsar timing parallax or very long baseline interferometry (e.g., PSR J04374715, see Deller et al. 2008; Verbiest et al. 2008). In this case, instead of taking the ratios of the various amplitudes, it becomes possible to directly invert equations (20) to obtain the auxiliary quantities and . As before, these additional constraints lead to much tighter constraints on the compactness and the combination of the angles and . The independent knowledge of one of these angles then completely solves the system of equations.

Finally, we consider in Figure 5 the case where the mass of the pulsar is known through other means, such as a dynamical measurement of the pulsar mass in a binary system. This is indeed a likely situation, as a large fraction of millisecond pulsars are in binary systems with white dwarf companions. Currently, there are 17 millisecond pulsars with precisely measured masses, including some of the NICER targets, and this number is growing rapidly (see a recent compilation in Özel & Freire 2016). This figure shows that combining the measurement of the compactness with a prior knowledge of the mass with an uncertainty that is comparable to that of PSR J04374715 (; Reardon et al. 2016) is sufficient to distinguish between different equations of state.

## 4. Conclusions

In this paper, we investigated the effect of instrument or sky background on the determination of neutron star compactness and radius from a measurement of the shapes of pulse profiles from rotation-powered pulsars. Specifically, we made use of analytic formulae for the pulse shapes of moderately spinning neutron stars in the Schwarzschild metric to connect the observables from the profile to the properties of the neutron star. We showed that, in the case of isotropic emission and not utilizing any distance or mass information for the pulsar, the ratio of harmonics can be used to constrain the neutron star compactness. Furthermore, we showed that knowledge of the background at a few percent level for a background-to-source countrate ratio of 0.2 and a total accumulation of counts allow for a measurement of the neutron star compactness to better than uncertainty for most of the parameter space. These uncertainties become even smaller for more favorable geometries where the harmonic amplitudes are higher, as is the case for prototypical NICER targets such as PSR J04374715. For these geometries, a uncertainty in the compactness is achievable even when the background-to-signal ratio becomes comparable to (this is, e.g., a factor of two worse than what is expected for PSR J00300451).

Polar caps of rotation powered pulsars emit thermally in the soft X-rays because they are heated from above by particle bombardment from the magnetosphere. In this case, it is expected from model atmosphere calculations that the emission from the surface will be beamed. This anisotropy allows for more fractional harmonic amplitudes to be measured compared to the isotropic case and can reduce the uncertainty in the compactness even further for the same requirements on the background. Similarly, utilizing distance information to the pulsar, when it is known precisely, and when the instrument’s effective area is sufficiently well calibrated, allows for a measurement of the absolute, rather than fractional, amplitudes, increasing the number of directly and independently determined observables to two or three, depending on the beaming of emission. Finally, a mass measurement for the pulsar obtained through pulsar timing enables a precise measurement of the neutron star radius when combined with the compactness measured from the pulse profile. All of these will lead to a more accurate determination of the dense matter equation of state through means that are different from the earlier spectroscopic measurements (see Özel et al. 2015).

## References

• Abdo et al. (2013) Abdo, A. A., Ajello, M., Allafort, A., et al. 2013, ApJS, 208, 17
• Arons (1981) Arons, J. 1981, ApJ, 248, 1099
• Arzoumanian et al. (2014) Arzoumanian, Z., Gendreau, K. C., Baker, C. L., et al. 2014, Proc. SPIE, 9144, 914420
• Bai & Spitkovsky (2010a) Bai, X., & Spitkovsky, A. 2010a, ApJ, 715, 1270
• Bai & Spitkovsky (2010b) Bai, X., & Spitkovsky, A. 2010b, ApJ, 715, 1282
• Bauböck et al. (2015) Bauböck, M., Psaltis, D., Özel, F. 2015, ApJ, 811, 144
• Becker & Trümper (1993) Becker, W., & Trümper, J. 1993, Nature, 365, 528
• Bogdanov & Grindlay (2009) Bogdanov, S., & Grindlay, J. E. 2009, ApJ, 703, 1557
• Bogdanov et al. (2007) Bogdanov, S., Rybicki, G. B., & Grindlay, J. E. 2007, ApJ, 670, 668
• Braje et al. (2000) Braje, T. M., Romani, R. W., & Rauch, K. P. 2000, ApJ, 531, 447
• Cadeau et al. (2007) Cadeau, C., Morsink, S. M., Leahy, D., & Campbell, S. S. 2007, ApJ, 654, 458
• Dai et al. (2015) Dai, S., Hobbs, G., Manchester, R. N., et al. 2015, MNRAS, 449, 3223
• Deller et al. (2008) Deller, A. T., Verbiest, J. P. W., Tingay, S. J., & Bailes, M. 2008, ApJ, 685, L67
• Finley et al. (1992) Finley, J. P., Ogelman, H., & Kiziloglu, U. 1992, ApJ, 394, L21
• Gendreau et al. (2012) Gendreau, K. C., Arzoumanian, Z., & Okajima, T. 2012, Proc. SPIE, 8443, 844313
• Harding & Muslimov (2001) Harding, A., & Muslimov, A. 2001, ApJ, 556, 987
• Harding & Muslimov (2002) Harding, A., & Muslimov, A. 2002, ApJ, 568, 862
• Johnson et al. (2014) Johnson, T. J., Venter, C., Harding, A. K., et al. 2014, ApJS, 213, 6
• Morsink et al. (2007) Morsink, S. M., Leahy, D. A., Cadeau, C., & Braga, J. 2007, ApJ, 663, 1244
• Nath et al. (2002) Nath, N. R., Strohmayer, T. E., & Swank, J. H. 2002, ApJ, 564, 353
• Ögelman (1991) Ögelman, H. 1991, NATO Advanced Science Institutes (ASI) Series C, 344, 87
• Özel (2013) Özel, F. 2013, Reports on Progress in Physics, 76, 016901
• Ozel (2016) Özel, F., & Freire, P. 2016, ARA&A, in press
• Ozel et al. (2015) Ozel, F., Psaltis, D., Guver, T., et al. 2015, ApJ, in press (arXiv:1505.05155)
• Pavlov & Zavlin (1997) Pavlov, G. G., & Zavlin, V. E. 1997, ApJ, 490, L91
• Pechenick et al. (1983) Pechenick, K. R., Ftaclas, C., & Cohen, J. M. 1983, ApJ, 274, 846
• Poutanen & Beloborodov (2006) Poutanen, J., & Beloborodov, A. M. 2006, MNRAS, 373, 836
• Psaltis & Ozel (2014) Psaltis, D., Özel, F. 2014, ApJ, 792, 87
• Psaltis et al. (2014) Psaltis, D., Özel, F., & Chakrabarty, D. 2014, ApJ, 787, 136
• Reardon et al. (2016) Reardon, D. J., Hobbs, G., Coles, W., et al. 2016, MNRAS, 455, 1751
• Ruderman & Sutherland (1975) Ruderman, M., & Sutherland, P. 1975, ApJ, 196, 51
• Strohmayer (2004) Strohmayer, T. E. 2004, X-ray Timing 2003: Rossi and Beyond, 714, 245
• Strohmayer et al. (1997) Strohmayer, T. E., Zhang, W., & Swank, J. H. 1997, ApJ, 487, L77
• Verbiest et al. (2008) Verbiest, J. P. W., Bailes, M., van Straten, W., et al. 2008, ApJ, 679, 675
• Weinberg et al. (2001) Weinberg, N., Miller, M. C., & Lamb, D. Q. 2001, ApJ, 546, 1098
• Yan et al. (2011) Yan, W. M., Manchester, R. N., van Straten, W., et al. 2011, MNRAS, 414, 2087
• Zavlin (2006) Zavlin, V. E. 2006, ApJ, 638, 951
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters