Could acceleration of a pulsar affect braking index?

In the coming era of multi-messenger astro-physics, pulsars might be one of the most possible electromagnetic counterparts of the gravitational wave. The braking indices, which are related closely to the electromagnetic radiation of pulsars, are shown to be larger for the pulsars with companion. It motivates us to set up a modiﬁed spin-down equation for accelerated pulsars. In this model, we attempt to ﬁgure out whether acceleration of a pulsar can cause a larger braking index.


Introduction
The gravitational wave, which is theoretically predicted over 100 years ago, has been finally detected recently [1][2][3]. In the coming era of multi-messenger astrophysics, it requires searches for electromagnetic counterpart of the gravitational wave than ever before. Moreover, the gravitational wave is related to electromagnetic radiation closely. For instance, gravitational wave event GW-170817 is produced along with gamma-ray burst [3] and pulsar timing array is used to detect gravitational wave.
Pulsars, especially binary pulsars, one of the most possible electromagnetic counterparts of gravitational wave, are known as highly magnetized rotational stars in the space. Pulsars have been observed and investigated extensively for several decades. Rotation velocity of a pulsar could be observed precisely. It shows that pulsars would rotate slower as time goes by. In the canonical model, pulsar can be simply regarded as a rotational magnetic dipole. Loss of rotation energy is due to the magnetic dipole radiation [4,5]. For a pulsar with moment of inertia I and magnetic dipole m, the evolution of rotation velocity Ω is described by the so-called spin-down equation, a e-mail: zhuqh@ihep.ac.cn (corresponding author)Ω where θ is magnetic inclination angle and c is speed of light. However, the spin-down equation does not describe the observed rotation velocity Ω well. The deviation of spindown equation can be indicated by a dimensionless quantity, the so-called braking index, In the canonical model, the braking index equals to 3, while almost all observed braking indices are beyond 3 [6,7]. The deviated braking indices indicate a modified spin-down equation.
In fact, various models were proposed to deal with the braking index problem. There are two major scenarios. The first scenario suggests that pulsars should have other kinds of radiation source besides the magnetic dipole radiation. And the second indicates that magnetic moment or moment of inertia of pulsars should evolve with time. In the first scenario, the expected radiation source might be gravitational quadrupole radiation [8][9][10], in which the pulsar is thought as an imperfect sphere. The energy radiation also can be caused by the outflow of relativistic particles in pulsar wind model [11][12][13], or fall-back disk around pulsars [14][15][16]. In the second scenario, the braking index can be given [6,17] by, The deviated braking indices turn to require physical origin forṀ,θ andİ . The magnetic moment, or the magnetic field, of pulsars might change with time due to Hall drift, Ohmic decay or other magnetic mechanism [18][19][20]. It was also studied statistically by providing phenomenological model for the evolution of magnetic field [21][22][23]. For the magnetic inclination angle, it might evolve with time tending to be in alignment or out of alignment because of the plasma in mag- The gray filled histogram is the braking indices of the rest 409 pulsars. The dark gray filled histogram is the cross region of these two kinds of pulsars netosphere [24] or the rearrangement of pulsar components [25], respectively. The inertia of moment related to equation of state of the pulsars was also considered [26]. From these models, the deviated braking indices around 3 seem to be well understood. However, as Refs. [22,27] pointed out, there are numerous braking indices beyond 3 over several orders of magnitude. Maybe, due to limitation of structure of pulsars or attributing to inaccuracy observation, the large braking indices are rarely explored. In this paper, we wish to be confronted with the braking index problem and focus on the large braking indices. In the ATFN pulsar catalogue (http://www.atnf.csiro.au/research/pulsar/ psrcat/), there are 441 pulsars whose braking indices can be obtained. 32 of them have companions. We find that the 32 pulsars tend to have larger braking indices than that of the rest. The statistical result is shown in Fig. 1. It motivates us to consider that the gravitational field around pulsars might cause the large braking indices.
In this paper, we deal with braking index problem by proposing a modified spin-down equation based on a moving dipole potential [28]. In the modified spin-down equation, the braking indices increase with the accelerations of pulsars. We expect that the gravitational field around a pulsar might be the possible origin of the acceleration. The paper is organized as follows. In Sect. 2, we show the magnetic dipole radiation of a moving pulsar. In Sect. 3, we derive the braking indices via constructing the modified spin-down equation and provide estimation of the accelerations of pulsars. The main conclusions and discussions are summarized in final Sect. 4.

Magnetic dipole radiation of a moving pulsar
In canonical model, the pulsars spin down due to the magnetic dipole radiation. Likewise, we only refer to the magnetic dipole radiation of pulsars. The difference is that we consider an accelerated pulsar based on a moving dipole potential [28]. The potential is of the form, where The p μ and m σ are electric and magnetic dipole, respectively. Source terms are functions of τ which is defined by the light cone condition, Thus, the point dipole in the retard potential also can be represented as The τ is function of x μ . One can calculate partial derivative of Eq. (5), the partial derivative of τ referred to x μ is of the form, Using Eqs. (4) and (6), we can obtain the electromagnetic field F μν of the dipole, where a μ ≡ du μ dτ ,ȧ μ ≡ da μ dτ andQ ≡ d dτ Q. The construction of spin-down equation requires the form of electromagnetic radiation. For simplicity, we calculate the radiation in the moving frame adapted to u μ . Namely, we set u = (c, 0) in the electromagnetic field F μν . And changes of space-time metric are not considered. We would discuss it in the Sect. 4. In the reference frame of pulsars, the Eq. (7) can be simplified, such as r μ u μ = r 0 u 0 = −cr and the vanished electric dipole. By making using of normalization of 4-velocity u μ u μ = −c 2 , we know that the time component of acceleration is zero.
In classical electrodynamics, radiation of electromagnetic field is derived from the electromagnetic energy-momentum tensor, And the radiation angular momentum, which is also the radiated electromagnetic torque, is calculated via the surface integral related to the spatial part of energy-momentum tensor at spatial infinity, The formula can be obtained without any ambiguities, although the calculation might be a bit cumbersome. Checking the case that a =ȧ = 0, one can obtain the canonical spin-down equation using the radiation torque. And the braking index of the canonical spin-down is always 3. In the case that the acceleration a cΩ 1 andȧ cΩ 2 1, the radiation energy of pulsars would not be affected too much. It would be shown in the next section that the accelerations of the case could still affect braking indices of pulsars in our model.

Braking index of accelerated pulsars
The spin-down equation describes time evolution of the rotation of a pulsar. In this section, we would construct a modified spin-down equation for an accelerated pulsar and study the corresponding braking indices.
Using conservation law of the energy momentum tensors, we know that loss of rotational angular momentum equals to the radiated electromagnetic torque. The modified spin-down equation is established as, For simplicity, we assume that the rotational magnetic dipole of a pulsar satisfies the equation, which means a pure rotation. For pulsars with companions, a ∼ |P b ×a| aΩ, where P b is orbit period. In the torque Eq. (9), theȧ term is shown to be unimportant. Thus, we assume the pulsars undergoing nearly uniform accelerations. The electromagnetic torque is of the form, (ṁ·(a × m)) + 4 15 a · m(a ×ṁ) Using the electromagnetic torque, we can write the modified spin-down equation as, With the modified spin-down equation, we can calculate the braking index, where |Ω| ≡ Ω. In principle, the braking index can be obtained by using Eqs. (11), (13) and (14). For the sake of intuition, we here would choose representable cases to tell the story. In these cases, the accelerations of pulsars are chosen to be parallel or vertical to the direction of rotation axes of pulsars.

Accelerations of pulsars parallel to direction of rotation axes
For the pulsars whose accelerations are parallel to the direction of Ω, the modified spin-down equations can be simplified as, a(Ω × m)·(a × m) · m)(a×(Ω × m)) . The modified spin-down equation can be rewritten as, In the case of a = 0, the Eq. (16) reduces to the canonical spin-down equation (1). The magnetic inclination angle evolves with time satisfying the equation, The differential equation can be solved analytically. Solution is of the form, where k ≡ m 2 I c 3 and b ≡ a c . They are both constant parameters in our models. The Eq. (18) shows that the magnetic inclination angle evolves towards θ = π 2 . The accelerations of this case could be the causes of oblique rotation of pulsars. Using Eqs. (14), (16) and (17), we obtain the braking indices, In this case, the braking indices can't be larger than 3. In the Fig. 2, we plot the braking index n(b, θ) as function of acceleration b for fixed magnetic inclination angle θ . The braking indices decrease with the accelerations of pulsars. And the braking indices drop to a minimum value as the acceleration tends to infinity.

Accelerations of pulsars vertical to direction of rotation axes
The other case is that the accelerations of pulsars are vertical to the direction of Ω. In the case, the modified spin-down equation can be simplified as, where α ≡ a · m, β ≡ a·(Ω × m). The modified spin-down equation can be rewritten as, whereα ≡ a·m am andβ ≡ a·(Ω×m) aΩm . For pulsars with Ω ∼ 1s −1 andΩ 10 −10 s −2 , it's reasonable to use an approximation that kΩ ≡ m 2 Ω I c 3 Ω Ω 2 1. In the approximation, the modified spin-down equation can be simplified further as, The magnetic inclination angle evolves with time satisfying the differential equation as follows, In the approximation,α andβ seem like trigonometric functions, since the differential equations referred toα andβ are given bẏ The rotation velocity Ω of pulsars can be solved numerically with Eqs. (23)- (26) in principle. Fortunately, what we focus is the braking index, which can be obtained without solving these differential equations, where f b Ω , θ ∼ 1. The function f turns to be unimportant for the most situations, since (kΩ) −1 1. The most interesting is that the positive or negative of n depends onα and β, which are related to position of pulsars and direction of the acceleration. Using Eqs. (25) and (26), we can obtain that α = sin φ andβ = cos φ approximately, where φ ≈ Ωτ . It leads to value of braking indices around the order of magnitude of b 2 kΩ 3 . We rewrite the braking index in the term of b, θ, k and φ without the unimportant terms, In Fig. 4, we plot the n(b, θ, k, φ) as function of the acceleration b for fixed θ, k and φ. The braking indices can be different from 3 when the acceleration a ≈ 1m · s −1 . It indicates that the braking indices are more sensitive to the acceleration whose directions are vertical to the rotation velocity of pulsars. And in this case, the absolute value of braking indices increase with the acceleration. In Fig. 5, we plot the n(b, θ, k, φ) as function of the magnetic inclination angle θ for fixed b, k and φ. It shows that the acceleration affect braking indices more apparently for alignment pulsars. The braking indices are also related to the k, namely, moment of inertia and magnetic dipole of pulsars. Pulsars with larger moment of inertia or smaller magnetic dipole would have larger braking indices. It has been shown that the large braking index can be caused by acceleration of a pulsar in our model. In this section, we would provide the order of magnitude estimation of the acceleration. In the estimation, we assume that deviation of braking indices is completely caused by acceleration of pulsars. The pulsar timing data comes from ATNF pulsar catalogue. When the direction of acceleration is vertical to the rotation velocity of pulsars, the braking indices can be changed more apparently than the other case. Therefore, we use Eq. (26) to estimate the accelerations of pulsars. As |α| 1, |β| 1 and | sin θ | 1, they turn to be unimportant for the braking index in the order of magnitude. The approximated Eq. (26) is reduced to the form, We consider the case that b Ω 1, which leads to kΩ 3 ≈Ω. In the approximation, we could estimate acceleration of pulsars completely with the pulsar timing data. The accelerations can be expressed as function of n and Ω, In Fig. 6, the accelerations of 32 pulsars with companions are shown. The estimated acceleration a ranges from 10 2 to 10 6 m · s −2 . For most of binary systems, pulsars have accelerations around 10 4 m · s −2 . It's the order of magnitude far away from the event horizon of black hole. In our model, the acceleration is covariant and involves the friction force.
In the frame of a moving pulsar, one can show that the friction force equal to external force always. Thus, our results on acceleration could indicate the external force, maybe gravity we expected, experienced by pulsars. The estimated accelerations for the 441 pulsars are shown in Fig. 7. For the pulsars with braking indices less than 100, the braking index seems not correlation with acceleration. It indicates that, for the small braking indices, there might be other effect for the deviated braking indices as reviewed in the introduction. For larger braking indices, the accelerations of pulsars have positive correlation with the braking indices overall. And the pulsars with companions tend to have larger estimated acceleration than that of the rest as we expected. Our estimation is meaningful for the case that a cΩ ∼ 10 8 m · s −2 , the pulsar (in the left-top of Fig. 7) with acceleration over 10 8 m · s −2 should be excluded.

Discussions and conclusions
In this paper, we discussed the possibility that acceleration motion of a pulsar can affect its magnetic dipole radiation. In detail, we derived the modified spin-down equation and showed that the accelerations of pulsars are closely related to Here, we only considered the most simple cases that the accelerations are along specific directions. As we known, it's not realistic. In this sense, our results are preliminary and just give an order of magnitude estimation. The braking indices could be affected by acceleration motions of pulsars, only if the accelerations are larger than 1 m s −2 .
We used the acceleration to indicate the external force experienced by a pulsar. However, the origin of the acceleration in our model still needs to be explored further, especially, for those pulsars without being marked as binaries. For the binaries, Keplerian parameters can be used to estimate the accelerations of pulsars via orbit period and mass of companions. The data of the accelerations can be obtained from ATNF pulsar catalogue. We found that the accelerations estimated from braking indices in our results is larger than that estimated by Keplerian parameters over several orders of magnitude. Of course, the acceleration experienced by a pulsar may be a combination of gravity from normal matter, dark matter and dark energy. In the future, we can test every component of the acceleration by more accurate observations.
As suggested by Lyne et al. [6] from observations, the braking indices that are beyond 3 over several orders of magnitude seem not very reliable. There might be problems from accuracy ofΩ or dependence of fitting models. And in the ATNF pulsar catalogue, there aren't data ofΩ for most pulsars. This limits us to the 441 pulsars. With more accuracy observations in the future, the studies about the lager braking indices would be more meaningful.
We calculated the braking indices in local inertial frames adapted to the 4-velocity for a pulsar. Firstly, This frame is comoving with a pulsar. Namely, there are not translation motions of pulsars in the calculation. In this situation, the radiating angular momentums are clear from physical point of view. Secondly, it neglects the changes of metric for the reference frames. In the uniformly accelerating frames, the only difference of the metric could be g 00 compared with the metric of Minkowski space-time, such as Rindler metric. The spatial part of energy-momentum tensors are nearly the same form as that in flat space-time. It indicates that our calculation might be valid at least in the case of uniform acceleration motion for a pulsar.