Perturbative correction terms to electromagnetic self-force due to metric perturbation : astrophysical and cosmological implications

We consider the equation of motion of a charged particle or a charged compact object in curved space-time, under the reaction of electromagnetic radiation and also consider a physical situation such that the charged particle or compact object emits gravitational radiation, thereby gravitational radiation reaction also acts on it. We investigate the effect of this metric perturbation i.e. the gravitational radiation on the electromagnetic self-force. We show that, besides the interaction terms derived by P. Zimmerman and E. Poisson (Phys. Rev. D 90, 084030, 2014), additional perturbative terms are generated, which are linear in metric perturbation and are generated due to perturbation of the electromagnetic self-force by the metric perturbation. We discuss the conditions of significance of these perturbative terms and also the interaction terms with respect to the gravitational self-force in various astrophysical and cosmological cases ; such as the motion of charged particles around black holes, some extreme mass-ratio inspirals (EMRIs) involving sufficiently accelerated motion of charged stars (specially neutron stars) or charged stellar mass black holes around supermassive black holes, and motion of charged particles around charged primordial black holes formed in the early Universe etc.. We find that in some astrophysical and cosmological cases these perturbative terms can have significant effect in comparison with the gravitational radiation-reaction term.


I. INTRODUCTION
The motion of a point charge in flat space-time was one of the main topics of research in physics from as early as 1930s. Many pioneering physicists like Lorentz, Abrahams, Poincare and Dirac contributed to the development of the subject from its early onset [1]. DeWitt and Brehme generalized Dirac's result to curved spacetimes, and, for the first time, gave a precise derivation of Electromagnetic self-force [2]. Later, Hobbs applied vierbein treatment to derive their equations and found that their results must be corrected by a term involving the Ricci tensor [3]. The rigorous derivation of the electromagnetic self-force was given by Samuel E. Gralla et al, in their work in 2009 [4]. Similar counterpart of electromagnetic self-force in gravity viz. the 'Gravitational self-force' also was derived, first by Mino, Sasaki, and Tanaka [5], and then by Quinn and Wald using a different method [6]. In this work, we consider an interesting case where both the electromagnetic and gravitational self-forces are present ; and both of them produce corresponding radiation reactions in the motion. We consider the equation of motion of a charged particle in curved space-time under electromagnetic radiaton reaction, with external Lorentz force too and consider a physical situation such that the particle emits gravitational radiation, which perturbs the surrounding space-time. We then follow the procedure of deriving the MiSaTaQuWa equation for this equation of motion. Here the term 'particle' does not strictly mean that it has very tiny size like elementary particles ; the mass should be centralized enough so that the equations of motion of pointparticles can be applied. In this sense, a compact object like a neutron star or a stellar mass black hole orbiting a supermassive black hole can also be treated by a point-particle equation of motion. Although there are problems with point-particle notion when terms with second order metric perturbations are considered in calculation, we will not be facing it as we are considering linear metric perturbations only. In this context the work by Peter Zimmerman and Eric Poisson [7] is very important, as it is probably the first work treating the case where electromagntic and gravitational self-forces act together, while deriving their interaction terms. The authors in this work have considered not only the electrovac space-time, but also a more general case of scalarvac spacetime, where the metric of the background space-time is a solution of the Einstein's equation in the presence a scalar field. We show that in comparison to the case where only gravitational self-force is present, the presence of electromagnetic self-force with it, adds not only the interaction terms but also several extra perturbative terms in the equation of motion.
If we see the order only in terms of the charge and mass of the charged particle, then the order of these terms are proportional to q 2 m for force and to q 2 for acceleration ; q being the electric charge of the particle and m being its mass. We shall discuss the issue of orders of terms in detail, in section III in this work. We interpret that these additional terms are just the perturbations to the electromagnetic self-force, due to the gravitational radiation or metric perturbations emitted by the system. We have considered these perturbative terms to the first order of the metric perturbations. We investigate some astrophysical systems and cosmological cases, where these additional perturbative terms produced from the electromagnetic self-force, are significant in comparison with the gravitational self-force. In this case it is to be noted that although the gravitational self-force and these perturbative terms have a distinct difference to the sense that the first one is a purely gravitational aspect, while the origin of the perturbative terms are electromagnetic ; yet, they have the similarity that they contain the metric perturbations (to the linear order in this case) as they are generated due to metric perturbations. The motive of this comparison is that this will help to identify the cases where the perturbative terms generated from electromagnetic self-force would have significance or would dominate over the gravitational self-force and viceversa. The equations of motion in the corresponding cases can be simplified or approximated accordingly. We find that although in a few cases it is possible that these perturbative terms can be significant in comparison with the gravitational self-force, the special interest comes out to be in the cases of charged particle's motions around primordial black holes within certain mass range, which were produced in early Universe by direct gravitational collapse of sufficiently deep density perturbations. We have discussed the orthogonality, with the four-velocity of the particle, of different terms present in the radiation reaction, in Appendix-1. There are other two Appendices. In Appendix-2, we have explained why we have neglected the perturbations originated from the term containing Riccitensors and in Appendix-3 we have discussed certain issues related to the comparison of different parts within the perturbative correction terms, with the gravitational self-force term. We have not discussed the additional perturbative terms generated from the electromagnetic 'Tail term' due to the metric perturbations, in this work, which will be left for a separate future work. In this work, we have preliminarily written the electromagnetic self-forces in Gaussian units with c (speed of light in vacuum ) = 1 ; but while estimating some numerical quantities related to it, we have converted this to S.I. units system.

II. OVERCOMING SINGULARITY OF THE RETARDED METRIC PERTURBATION ON THE WORLD LINE OF THE PARTICLE AND THE GAUGE FIXING :
One main issue of our work is that the physical or retarded metric perturbation h ret µν emitted from the charged particle is singular at the particle, or in other words, it is singular on the world-line of the particle. As our work deals with the world-line of the particle, we must address this singularity of the metric perturbation emitted form the particle. For overcoming the problem of singularity of the retarded metric perturbation on the world-line of the particle, we follow the Detweiler-Whiting formalism [8], in which Detweiler and Whiting proposed a reformulation of the perturbed motion, where instead of breaking the overall retarded perturbation into 'Direct' and 'Tail' parts, they decomposed it into 'Singular' and 'Regular' parts : where the singular part h S µν is responsible for the singular behaviour of the retarded perturbation on the world-line, while it does not generate any self-force or does not affect the motion of the particle. On the otherhand, the regular part h R µν is a smooth solution of the perturbation equations and this exerts the identical self-force, as generated by the overall retarded perturbation. The Detweiler-Whiting formalism indicates the fact that the particle effectively moves along a geodesic 1 of a smooth perturbed space-time with the metric g ′ µν = gµν + h R µν [8][9][10], where gµν is the unperturbed metric. On the world-line of the particle, the regular part of the perturbation satisfies [10] where theh T ail µν;λ is given by (in trace-reversed form) : (3) Therefore, we shall only work with the regular metric perturbation h R µν , thereby eliminating the singular behaviour of the retarded metric perturbation, emitted from the particle, on its world-line and all the consequent general relativistic perturbation quantities will be in terms of h R µν . In this work, it is pertinent to be mentioned that, although in case of electromagnetic self-force, the point particle approximation is valid, but in case of general relativity, the point particle concept fails at non-linear orders. But, as the point particle concept has no problem with the linear orders of metric perturbation, which we are considering here, we can stick to this concept. However, it is worth noting that if any calculation is necessary in the non-linear orders of metric perturbation, where instead of point-particle notion the object of interest is a compact one, one has to employ a different method, known as 'Puncture method' [11][12][13]. In this method the retarded metric perturbation is divided into two parts known as 'Puncture part' h P µν and 'Residual part' h R µν . For detailed discussion on this method, previous references [11][12][13] and the review article by L. Barack and A. Pound [14] can be consulted. We donot go into detail about this 'Puncture method' here, as it is not required for our analysis done with first-order metric perturbations. Another issue in our work is that as the gravitational selfforce and the metric perturbation are both gauge-dependent quantities, we have to fix the gauge to describe this physical effect meaningfully. We here choose it to be the Lorenzgauge, in which many preliminary and foundational results in gravitational self-force had been obtained. In terms of the trace-reversed form of the metric perturbation : (h ret = g αβ h ret αβ ), the Lorenz-gauge condition is given by : Here, a confusion may arise that we will be working with the regular part of the perturbation h R µν , while the gaugecondition is in terms of the trace-reversed version of the overall perturbation h ret µν . But, there is no problem with this issue, as what happens is actually that the singular part h S µν remains invariant under a smooth gauge-transformation, while only the regular part h R µν changes. Furthermore, in our case, as there is electromagnetic radiation emitted from the charged particle in the space-time surrounding it and in the most general case there is also the external electromagnetic field, hence the metric actually satisfies the Einstein-Maxwell equation, not the vacuum Einstein equation. If we consider an EMRI, where the smaller component is charged, then the metric outside the larger massive body in the EMRI satisfies the Einstein's equation : where the energy-momentum tensor Tµν of the source outside the larger massive body contains two components : the energy-momentum tensor of the smaller charged massive body itself and the energy momentum tensor of the electromagnetic radiation emitted from the smaller charged massive body (if there be any external electromagnetic field then its energymomentum tensor is also to be added). With the pointparticle approximation, the energy momentum tensor of the smaller massive body can be represented by a Dirac-Delta function. 2 Hence, outside both the larger and smaller massive bodies, the only energy-momentum is of the electromagnetic radiation, or in other words, the metric there satisfies the Einstein-Maxwell equation : where µ0 is the permeability in vacuum and Fµν is the electromagnetic field-strength tensor. So, the metric perturbations outside the larger massive body would satisfy the first order perturbation-equation of the background equation 7 .

III. ORDER OF THE PERTURBATIVE TERMS:
Samuel E. Gralla et al gave a rigorous derivation of the electromagnetic self-force in their work in 2009 [4]. In their work, they followed an approach called "asymptotic selfsimilar manner". In this approach, about the worldline of the charged body as λ → 0 (λ is a small parameter measuring the size of the charge and mass, not a parameter along the worldline [4]), the charge q and mass m of the charged body (of whose equation of motion is to be studied) tends to zero, but the charge-to-mass ratio q/m tends to a well-defined limit. This approach is considered mainly to tackle the difficulty with the point particle description of the charged body, when the limit is taken to zero-size in a straightforward way and also to avoid the problems associated with the body's finite-size consideration. As a result of following this asymptotic self-similar approach, in most of the cases the scaling of q and m are treated on equal footing i.e. it is presumed that q ∼ m, which may not be the case in some astrophysical and cosmological scenarios. Furthermore, when the order of a certain term in the equation of motion of the charged object is spoken about, it is often estimated and compared with the other terms only in terms of q and m. But in various cases, any other physical quantity present in that certain term, e.g. the four velocity, four acceleration and rate of change of four acceleration etc., may have such a huge order that they must have to be considered. Otherwise the estimation or comparison of the orders of these terms would turn out to be just incorrect. In our work, we discuss some astrophysical and cosmological cases where the orders of these quantities like four velocity, four acceleration or rate of change of four acceleration are so huge that the overall orders of the concerned terms can not be judged only in terms of q and m. In fact, in this work we study all the perturbative terms which are of linear order in metric perturbation h R αβ , instead of designating the order of perturbation both in terms of the electromagnetic perturbation and gravitational or metric perturbation. Generally the electromagntic self-force is seen as a perturbation over the external lorentz force and hence it may seem that the electromagntic self-force can not be compared with the gravitational force viz. the main newtonian-part of the gravitational force. But, quite recently A. Tursunov et al has shown [15] that for a charged particle with charge q and with relativistic speed, the electromagnetic self-force, which is of the order of ∼ q 4 B 2 /m 2 , can have same order of magnitude as that of Newtonian gravitational force (of the order ∼ GM m/r 2 ), when the charged particle is moving around a supermassive black hole of mass ∼ 10 9 M⊙ ( M⊙ is the usual symbol of Solar-mass), in presence of a magnetic field of B ∼ 10 4 G. 3 So, as it is possible in a practical case that the electromagnetic self-force can be of same order of magnitude with the newtonian gravitational force, it may also be possible for the perturbations in the eletromagntic self-force caused by metric perturbations to have similar orders of magnitude as that of the gravitational self-force. At last it is of utter importance to remember the fact that we have considered here metric perturbations originated due to the motion of the charged body itself. If there is external metric perturbation, when we speak about the orders in terms of charge and mass only, then the perturbative terms of the electromagnetic self-force originated from that external metric fluctuations would not have the order q 2 m, instead they would have the order q 2 , same as that of the electromagnetic self-force (As the h αβ then would not be of the order of m). So, in that case these terms would have similar order with the gravitational self-force, which is of the order of m 2 . Although in that case, for determination of the correct order would require the knowledge of the mass of that external souce of the gravitational radiation. In the case of external gravitational wave, we shall not need the Detweiler-Whiting reformulation to break the metric perturbation into singular and regular parts, as the external metric perturbation will not be singu-lar on the world-line of the particle. However, in this work we shall stick to the case where the gravitational radiation is generated from the charged particle itself.

IV. THE EQUATION OF MOTION OF A CHARGED PARTICLE IN CURVED SPACE-TIME UNDER ELECTROMAGNETIC RADIATION REACTION AND APPLYING METRIC PERTURBATION TO IT :
Here we consider the explicit form of the equation of motion of a charged particle or compact object in curved spacetime under the reaction of the electromagnetic radiation emitted by itself [10,15] and also consider the reaction of gravitational radiation, generated due to the motion of the particle around a comparatively very bigger massive compact object, preferably a black hole. The gravitational radiation emitted from the system creates a perturbation of space-time and that would have an effect on the motion of the particle, which is the reason behind the gravitational radiation reaction. Let τ ′ be the proper time associated with this perturbed metric of the particle, and τ the proper time for the unperturbed metric, without the reaction of gravitational radiation. We denote the unperturbed and perturbed metric as gµν and g ′ µν respectively, and physical or 'retarded' part of the metricperturbation (without the 'advanced' part of it), i.e. here the gravitational radiation emitted from the system, as h ret µν . To tackle the singularity of the metric perturbation on the world-line of the particle, we take the perturbed metric on the world-line as the effective metric : g ′ µν = gµν + h R µν , as we have already explained in the previous section 2. The equation of motion of the charged particle in the perturbed metric is given by : Here, Du ′µ dτ ′ is the covariant derivative of the particle's 4velocity with respect to τ ′ , given by On the right-hand side of Eqn. (8), the first term is the Lorentz force acting on the particle, the second term is the electromag-netic radiation reaction in curved space time, the third term is due to the interaction of the particle with the surrounding matter (if there be any) and the fourth term is the 'tail term' of the electromagnetic radiation in curved space-time. The quantity D 2 u µ dτ 2 can be expanded as [15]: Now, we designate the contribution of the metric perturbation by an additive vector a µ 4 in the equation of motion of the particle in unperturbed metric, in the following way (for the method used here, Ref. [9] may be consulted): Further substituting the expression of the first and second order covariant derivatives of the four-velocity of the particle, into the above Eqn. (11), we obtain Next, we substitute the operators d dτ , d 2 dτ 2 and d 3 dτ 3 in the above Eqn. (12) with the similar ones with respect to τ ′ ; and obtain : 4 It is to be noted very carefully that it can not be called solely a gravitational radiation reaction term. The reason for this will be clear at last, when we shall get its expression.
We substitute the expression of the d 2 x µ dτ ′2 in the LHS of the above Eqn.(13) from the Eqn. (8) and hence obtain : Now we arrange the equation in such a way that the similar terms in the perturbed and unperturbed metric come together so that it can be identified. Doing so Eqn. (14) can be written in the following form : From now, we shall neglect the terms containing Ricci tensors, due to interaction with surrounding matter. The reason, for which we neglect these terms, is explained in detail in Appendix-2.
We simplify different terms as differences between quantities in the unperturbed metric i.e. with respect to proper time τ and in the perturbed metric i.e. with respect to proper time τ ′ as follows : where the terms − 1 2 qF µ ν u ν u α u β h R αβ and −q(g µν + u µ u ν )h R να F α β u β have already been derived by P. Zimmerman and E. Poisson in an earlier work [7].
For brevity, we introduced the quantity In the above equation the quantity ∆Γ µ νρ is the part of the perturbation in the Christoffel symbol tensor, caused by the regular part h R µν of the metric perturbation and is given by : The above expression of ∆Γ µ νρ can be further simplified as : V. EXTRA TERMS GENERATED DUE TO PERTURBATION OF THE ELECTROMAGNETIC SELF-FORCE AND THEIR SIGNIFICANCE : In this section, we investigate the significance of the additional perturbative terms generated due to metric perturbations from the electromagnetic self-force in comparison with the gravitational self-force. The gravitational radiation reaction on the motion of the particle in curved space time in the absence of electromagnetic radiation reaction is given by : It is to be noted that the above simplified expression of the gravitational radiation reaction term can be obtained by applying the orthogonality property of the reaction, in the case where there is no electromagnetic radiation reaction [9]. Now, the additional perturbative terms, in linear order of the metric perturbation h R αβ , generated from the electromagnetic radiation reaction due to the metric perturbations emitted from the particle, which are absent when there is only one among electromagnetic self-force and metric-perturbations, are : and In this case, there may be a confusion that why we have written these correction terms in equations 23 to 26 separately, although their basic-source is same : the second term on the R.H.S. of the equation 8 i.e. the Abraham-Lorentz-Dirac term. The simple reason behind this is that despite having identical source, these terms originate from four different kind of terms within the Abraham-Lorentz-Dirac term. This can be clearly checked from the equation 10, where we have written the detailed expression of D 2 u µ dτ 2 , expanding the covarinat-derivatives within it.
In the next subsections we analyze the significance of the terms a µ 3 , a µ 4 and a µ 5 with respect to the gravitational radiation reaction term a µ 1 . We shall avoid discussing the significance of the term a µ 2 , as it contains time-rate of change of acceleration and hence this is very complicated to compare for practical astrophysical and cosmological phenomena.
A. The significance of the term a µ 3 : Let us now analyze the ratio : We see that for a µ 3 ∼ a µ 1 , one of the requirement is: After substituting u µ uη ≡ dτ ′ dτ 2 u ′µ u ′ η and cancelling out the u ′µ from both sides of the above condition, we obtain : It is to be noticed that the above relation is a tensorial one where on both sides three indices α, β and η are repeated indices and they are contracted among the tensors in such a way that we can not cancel out the term ∆Γ η αβ u ′α u ′ η from both sides of Eqn. (29), although that is common. For that, we write the expanded expression of the quantity ∆Γ η αβ du ′β dτ ′ u ′α u ′ η with respect to the repeated index β : Where the indices r, θ, φ and t denote the four coordinates of the coordinate system. In a similar way, we expand the quantity ∆Γ η αβ u ′β u ′ η u ′α with respect to β, and then substituting the expanded forms of these two quantities in both sides of the relation (29) to obtain the condition to be satisfied in coordinate-wise manner, we get 5 : Let us consider any one of the above set of relations, say the first one based on radial coordinate r. It gives : provided that the quantity ∆Γ η αr u ′ η u ′α is non-zero, which should be obvious for any curved space time as any of the spatial components of the four-velocity must be non-zero for the motion of the particle or compact object, as well as any of the perturbed component of the Christoffel. Thus, the condition (34) in S.I. units reads: 6 We describe the fulfillment of the condition (35) in two different classes of charged objects : (i) charged sub-atomic particles (e.g. we estimate it numerically for a proton) and (ii) charged neutron stars or stellar mass black holes. If we consider the orbiting particle to be a proton, the above condition yields du ′r dτ ′ u ′r ∼ 10 26 s −1 .
Therefore, for a proton orbiting around a black hole the condition for the significance of the perturbative term a µ 3 is that the acceleration of the proton has to be 10 26 order larger than the speed 7 . It would be interesting to test the significance of the term a µ 3 in the cases of charged stars, specially charged neutron stars or white dwarfs, or charged stellar mass black holes revolving around a supermassive black hole. Although, still there is no distinct astronomical evidence for any compact object containing a significant amount of net electric charge, specially for stars and black holes, many researchers have been working on theoretical models of stars containing a significant amount of net electric charge. For instance, in the reference [18] the authors have discussed a class of static stellar equilibrium configurations of relativistic spheres made of charged perfect fluids, where they have analyzed the physical acceptability of their theoretical model for some compact star candidates like SAX J1808.4-3658, 4U 1538-52, PSR J1903+327, Vela X-1 and 4U1608-52. They have concluded that their results strongly suggest that a class of compact stellar models with charged perfect fluid matter distribution is permitted with the new solution discussed in their work.
According to the references [19][20][21], and [22], the global balance of forces allows a net charge as large as 10 20 C in neutron stars, producing a very high electric field of order ∼ 10 21 V /m. Then the condition (35) for that would be : 6 As, we have considered the metric perturbation to be sufficiently small, so that the perturbation terms are retained up to linear order only ; hence the difference between the corresponding proper times of perturbed and unperturbed metrics viz. τ ′ and τ should also be small enough such that dτ ′ dτ ∼ 1 ). 7 Here we use proton just to get an idea on the order of acceleration required for such sub-atomic particles to satisfy the condition 35.
where MNS and M⊙ denote the mass of the neutron star and the solar mass, respectively. We know that usually neutron stars and white dwarfs have mass of the order of solar mass M⊙ i.e. MNS/M⊙ ∼ 1. Here, we consider the observational work in the reference [23], which reports that the speed of a star around the supermassive black hole at the center of Milky-way reaches approximately 8 × 10 6 ms −1 , and hence, in such cases, for satisfying the condition derived above, the acceleration of the star around the supermassive black hole has to be ∼ 10 12 ms −2 . So, achieving acceleration of this order would be quite difficult if a single neutron star or stellar mass black hole revolves around a supermassive black hole.
To satisfy the condition (37), a different type of astrophysical configuration is required. When a stellar mass black hole binary or a neutron star binary or a neutron star-black hole binary would be revolving around a supermassive black hole, then the smaller components in binary formation within the three-body system should achieve the required acceleration to satisfy the condition (37). It has been shown in the work by Xian Chen et al [24], that such EMRIs are expected to be produced by tidal capture of smaller binaries by a supermassive black hole. Again, at the late inspiral stage or at the merging stage, this binary becomes sufficiently compact with respect to the supermassive black hole such that it can be treated with the point-particle equations in our work.
In support of the fact that the smaller components of these EMRIs can achieve such order of acceleration required for satisfaction of the condition (37), we give an example of the acceleration in the binary black hole candidate GW150914, from which first direct detection of gravitational waves by aLIGO had been done [25]. For this candidate, during 0.2 second time interval of the detectable gravitational wave signal, the relative orbiting velocity of the black holes increased from 30% to 60% of the speed of light, and hence in this case the order of acceleration was approximately 10 8 m.s −2 . Although the merger-stage dynamics of the black holes can be accurately determined by numerical general relativistic techniques only, yet from this estimation of acceleration, we get an intuitive idea that in case of typical binaries of stellar mass black holes or neutron stars or binaries of black hole-neutron star, the acceleration achieved in the late inspiral stage and merger stage would be of similar order or even more. So, if the smaller component of an EMRI be such a binary, where the stars or stellar mass black holes contain sufficient net electric charge, then satisfying the condition (37) is clearly possible.
Therefore, binary formations of neutron stars or stellar mass black holes containing a net electric charge of order 10 20 C and inspiralling around supermassive black holes, are expected to satisfy the condition of significance (37). This type of extreme mass-ratio inspirals (EMRIs) are expected to be detected by the upcoming space-based gravitational wave detector LISA. Hence, if the smaller mass components of such EMRIs contain a significant amount of net charge, then neglecting the term a µ 3 , generated due to perturbation of electromagnetic radiation reaction by metric fluctuations or gravitational radiation, may lead to theoretically wrong estimation of the parameters related to these sources of gravitational waves.
Next, we consider a period of early universe where primordial black holes (PBHs) are expected to be produced by direct gravitational collapse of sufficiently deep density perturbations.
Furthermore, during this epoch, as the universe was full of charged particles (i.e. atom formation did not start yet), it is unlikely that the PBHs would be neutral. Therefore, in this case charged particles would inspiral around charged PBHs and ultimately fall into the PBHs, emitting gravitational and electromagnetic radiation. The interesting fact here is that many of these systems could have the size of atoms. For instance, the PBHs which had mass smaller than ∼ 10 20 kg, their Schwarzschild radius would be less than 10 −7 m, and thus, we expect that they can be treated as quantum particles. Therefore, in such systems we expect that the charged particles orbiting around charged PBHs would have huge acceleration, as required by the condition (36), and consequently, the perturbative term a µ 3 is expected to be significant 8 .
Another astrophysical phenomenon where the condition (36) may be satisfied is in 'Relativistic Astrophysical Jets'. In these case, accelerated ionized matter are emitted in the form of a beam from some high-energy astrophysical sources and usually the magnitude of their acceleration is huge. If in any of such astrophysical jet, the ions are accelerated during sufficiently small time to relativistic speeds, then the condition (36) is expected to be satisfied in the part of the jet closest to the source, specially a supermassive black hole at the center of an active galaxy. Even if the accelerated ions of such a relativistic astrophysical jet passes through the vicinity of another black hole or compact object, we may also expect that the condition (36) to be satisfied. 9 For a µ 3 ∼ a µ 1 , another requirement is: Analyzing the condition in coordinate-wise manner as we did before , we obtain for the radial coordinate: At this point, we make use of Eqn. (18), and substitute just the first term in the expression of ∆Γ µ rβ into Eqn.(39), and we get 2q 2 m du ′β dτ ′ ∼ u ′β .
(40) 8 An example of a system where we can have a huge acceleration is the revolving of an electron around a nucleus in the Bohrmodel of atom. The order of magnitude of such acceleration is ∼ 10 22 ms −2 . Hence, as in the early Universe the charged particles revolving around charged PBHs constituted systems of atomicsize, emitting both gravitational and electromagnetic radiation, there also the acceleration of the revolving particle would be quite similar, even expected to be larger due to the curvature of the PBH in comparison with nucleus of a typical atom. 9 In any astrophysical scenario, if the the Plasma acceleration of ions can be achieved, then that would be an ideal case for satisfying the condition (36). Indeed, in plasma acceleration of ions, the magnitude of acceleration as high as 10 22 − 10 23 ms −2 can be reached [26,27].
It is interesting to note that finally we obtain the identical condition given in Eqn.(34) (or, equivalently, in Eqn.(35) ), and hence the practical cases where this would be satisfied are also same.
Next, we compare the left hand side of Eqn.(39) with the part of ∆Γ µ rβ involving derivatives of h R µν , and we find It is to be noted that whether the above condition would be satisfied in any case, would depend on the associated components of the regular part of gravitational radiation h R µκ and metric g µκ ( more specifically saying it would depend on the index µ, as the rest of the indices are repeated indices). Here, we are giving a qualitative discussion on satisfying the condition, rather than a quantitative analysis. Actually it depends on the fact that how does the regular part of the gravitational radiation vary with radial distance from the source. On the other hand, if we consider spherically symmetric metrics, then the partial derivative with respect to r of g θθ and g φφ would give a factor of 2r, while those of gtt and grr would depend on the particular type of that metric. However, it is well known that if the radial distance r is not too small, then the amplitude of gravitational radiation (and the regular part of this gravitational radiation too) acting here as the metric perturbation, is very much less than that of the metric components. Hence, for typical astrophysical systems, the values of quantities h R µκ and (∂rh R κβ + ∂ β h R κr − ∂κh R rβ ) should be very small. While, the values of quantities g µκ and (∂rg κβ +∂ β gκr−∂κg rβ ) are relatively very large than them. Hence, depending on the situation it is possible to get some cases where the condition (41) is satisfied.
B. The significance of the term a µ 4 : Next, we take the ratio of the perturbative term a µ 4 with the gravitational radiation reaction term a µ 1 : and one of the conditions for a µ 4 ∼ a µ 1 reads Following the same steps as in the previous subsection, we find 2q 2 3m u ′γ ∂γg µr ∼ g µr .
If we consider the metric to be Reissner-Nordstrom metric, then as the metric is spherically symmetric and timeindependent, the above condition (44) yields : It is very interesting to find that not only for g rr , but for all non-zero components of the Reissner-Nordstrom metric, the above condition gives in S.I. units: where r is the radial distance from the black hole. For protons, the above condition yields |u ′r | ∼ 10 26 s −1 r(in m).
Therefore, we expect the term a µ 4 to be significant for ultrarelativistic motion of protons around PBHs of mass less than 10 9 kg, which is equivalent to Schwarzschild length scale of 10 −18 m.
C. The significance of the term a µ 5 : Now, we check the significance of the term a µ 5 with respect to the gravitational radiation reaction term a µ 1 . Like the previous two terms, we first observe the ratio : So, one of the conditions for a µ 5 ∼ a µ 1 is given by: Following the similar steps we carried out previously, we explore this condition in coordinate-wise manner, for the Reissner-Nordstrom black hole carrying constant charge 10 , and after some algebra we obtain Taking ξ ≈ 1, and using grr = (1 − rs r + r 2 Q r 2 ) −1 , the above condition becomes 11 Where, rs is the Schwarzschild radius of the black hole given by rs = 2GM/c 2 and rQ is the length scale associated with the electrical charge Q of the black hole, given by rQ = Q 2 G/4πǫ0c 4 . For practical cases of charged particles or compact objects moving around charged black holes, 12 r is larger than both these length scales. However, the situation can be classified into two different cases : (i) r > rs, rQ, but yet r ∼ rs, rQ, and (ii) r >> rs, rQ. In the former case, the quantity rs Q r 2 is of order unity, whereas in the later case it would be very large. We here focus on case (i), because if the particle moves very far away from the black hole then the effect of gravitational wave on its motion would be very little. Thus, condition (51) reduces (in S.I. units) to This condition is similar to the one in Eqn.(46), and, hence, the cases of satisfaction of this condition would also be identical as discussed earlier.
Another condition of significance of a µ 5 with respect to a µ 1 can be given by : We expand the repeated index β on both sides of the above condition, and follow similar steps using the condition coordinate-wise, as described in previous cases. Then we take the condition for the radial coordinate for the repeated-index β and from that, cancelling the quantity u ′α ∆Γ µ αr from both sides, we get 2q 2 3m u ′ρ u ′σ Γ r ρσ ∼ u ′r .
If we expand the sum denoted by contracted repeated index ρ, in the LHS of above condition (54), then we obtain : If we consider the Reissner-Nordstrom metric, then on the LHS of condition (55), the only non-zero Christoffel symbol tensor component would be Γ r rr and hence, For Reissner-Nordstrom metric, the Γ r rr is given by : So, to find out the practical cases where the condition (56) holds, we need to first analyze the value of the Christoffel symbol component Γ r rr around typical black holes. The value of the quantity 3m 2q 2 , which is actually 6πǫ 0 c 3 m q 2 in S.I. units, on the RHS of condition (56), is of the order of 10 26 for proton. We have already calculated and used the similar quantity in condition (35). There would be a limitation on the radial component of the four velocity |u ′r |, as it can not exceed the speed of light in vacuum i.e. c. Due to this limitation on |u ′r |, the order of the value of Γ r rr must be higher than 10 18 for satisfying the condition for the case of proton.
Note that, in the expression of Γ r rr given in the Eqn.(57), for any charged particle or charged compact object orbiting around a typical Reissner-Norstdorm black hole, the radial distance r must be greater than the outer horizon, which is given by r+ = 1 2 (rs + r 2 s − 4r 2 Q ). For stellar mass black holes, the value of Γ r rr would decrease with increasing r and the same happens for supermassive black holes. But, for primordial black holes (PBHs), specially the ones having mass less than 10 20 kg, exactly the opposite happens. For instance, a PBH with mass ∼ 10 20 kg would have Schwarzschild length scale ∼ 10 −7 m, and thus for any charged particle revolving around the PBH at a radial distance not exceeding 10 −3.5 m, the value of Γ r rr would be less than unity and this would increase with the decrease of the radial distance. Therefore, for systems where relativistic charged particles revolve around charged PBHs of sufficiently smaller masses, created in early Universe, the condition (56) could be easily satisfied resulting in significance of the term a µ 5 . In those cases even when the speeds of the charged particles is about 10 2 ms −1 (while it is expected to be near the speed of light for smaller particles like proton), it would also not be difficult for the system to satisfy the condition (56).
D. The significance of the term a µ int 1 : The ratio of the interaction term a µ int 1 to the gravitational radiation reaction term a µ 1 is : So, for a µ int 1 ∼ a µ 1 , the requirement is : Using the first part of the expression of ∆Γ η αβ given in the equation 18, in the above condition 59, we have : 60) As we did in previous cases, writing the above condition 60 in coordinate-wise manner for α, β = r, we have : For Reissner-Nordstrom metric, the above condition 61 can be further simplified to : If we expand the repeated index η on both sides of the condition 62 ; then again following the coordinate-wise manner and chosing the case of η = r only, it reduces to : For the Reissner-Nordstrom metric, we have already evaluated the quantity g rr ∂rgrr in the previous sub-section with grr = (1 − rs r + r 2 Q r 2 ) −1 . Inserting the expression of g rr ∂rgrr in the above condition 63, we obtain : We have already discussed about the quantity in the previous subsection, while discussing the condition of significance of the term a µ 5 . For the motion of a charged particle or a compact object near a Reissner-Nordstrom black hole, if we consider r ∼ rs, rQ, the the above quantity is ∼ 1. In that case, the above condition 64 simplifies to : In S.I. system of units we write this condition 65 as : For protons, the value of the quantity q 4πǫ 0 c 3 m is ≈ 3.19 × 10 −8 C −1 s and hence the condition 66 gives : So, for protons or charged particles of similar category, it is expected that the above condition 67 is satisfied, if the charged particles move in the vicinity of a stellar-mass black hole with sufficient speed. Although, till now there is not probably any distinct observational evidence of electrically charged black holes, yet there can be external electromagnetic field, generated from any other source, near some astrophysical black holes and if this external electromagnetic field is sufficiently strong to satify the condition 67, then the interaction term a µ int 1 will be significant.
E. The significance of the term a µ int 2 : From the ratio : we can say that one of the conditions for a µ int 2 ∼ a µ 1 is : Using the first part of the expression of ∆Γ µ αβ given in the equation 18, in the above condition 69, we have : Now similar to the previous cases, considering the above condition 70 in coordinate-wise manner for α, β = r, we get : For diagonal metrics like Reissner-Nordstrom metric, the above condition 71 will be further simplified to : or, q m Frr ∼ u r ∂rgrr .
So, for Reissner-Nordstrom metric, the above condition 73 becomes : For charged particles or charged compact objects moving near the black hole in such a way that r ∼ rs, rQ, the quantity rs r − in the R.H.S. of the above condition 74 is ∼ 1. Then in that case, the above condition simplifies to : or in S.I. system of units, it can be written as : Similarly as before, substituting the value of the quantity q 4πǫ 0 c 3 m for proton, we write the condition 76 as : which is quite similar to the previous case of the condition 67. So, for the same cases, as discussed in the prvious subsection, this condition too would be satisfied.

VI. CONCLUSION AND DISCUSSION :
In this work we have shown that coexistence of metric perturbations and electromagnetic self-force can lead to an effect, in the motion of charged particles in curved space-time, which does not exist when any one among these two is absent. In most of the physical situations, the metric perturbation is the gravitational radiation emitted from the charged particle itself, which causes the gravitational self-force. The perturbative terms of the electromagnetic self-force, which we have derived, are different from the interaction terms of electromagnetic and gravitational self-forces, given in the work of P. Zimmerman and E. Poisson [7]. We have analyzed different conditions for which these perturbative terms generated from the electromagnetic self-force due to its perturbation by the gravitational radiation would be significant in comparison with the gravitational self-force. We have also analyzed the conditions of significance of the interaction terms of electromagnetic and gravitational selfforces in comparison with the gravitational self-force. It is interesting to find that there are astrophysical phenomena and cosmological cases where these perturbative terms can play a significant role. The physical interpretation of these perturbations to the electromagnetic self-force by the gravitational radiation can be understood as the fact that, when electromagnetic self-force is acting in curved space-time, where there is gravitational wave emission from the system, then the electromagnetic wave produced due to the motion of the charged particle has to traverse through the ripples in the curved space-time due to the gravitational radiation. However, in the absence of the gravitational radiation, the electromagnetic wave propagates through the curved space-time but it does not face ripples in space-time. It is this difference with the case where gravitational radiation is present, that manifests in the form of these perturbative terms generated in the equation of motion of the charged particle. It is important to note that not taking into account these perturbative terms in the specified astrophysical or cosmological phenomena involving relativistic charged particles or compact objects, and considering only the gravitational radiation reaction can lead to incorrect estimation of their motions. In this way, we have also demonstrated that it can be misleading if we estimate or compare magnitude of terms in the equation-of-motion of the charged particle, only in terms of q and m. The other physical quantities present in the terms also matter. They can play a significant role in determining the overall order of certain term. We have not dicussed the case where the source of gravitational radiation perturbing the system is external. But, from the study in this work we can say, that for a suitable external source of gravitational radiation too, the perturbative terms of the electromagnetic self-force can be significant in comparison with the gravitational self-force.

VII. ACKNOWLEDGEMENT
Arnab Sarkar thanks S. N. Bose National Centre for Basic Sciences, Kolkata 700106, under Department of Science and Technology, Govt. of India, for funding through institutefellowship. Amna Ali is thankful to UGC for providing financial support under the scheme Dr. D.S. Kothari postdoctoral fellowship. Arnab Sarkar is grateful for the valuable advice and comments of Dr. Adam Pound, Royal Society University Research Fellow within Mathematical Sciences at the University of Southampton, the UK. Dr. Adam Pound's suggestions have helped greatly to improve this work and it is a privilege to get his opinions. Also, Arnab Sarkar wants to thank International Centre for Theoretical Sciences (ICTS) TIFR, Bengaluru 560089, India, as the motivation for this work came after attending a lecture series on 'Self-force and radiation reaction in general relativity' by Dr. Adam Pound, as a part of the programme -" Summer School on Gravitational Wave Astronomy 2019 " (Code: ICTS/gws2019/07).

VIII. APPENDIX-1 : ORTHOGONALITY OF THE RADIATION REACTION TERMS WITH FOUR-VELOCITY :
A. The orthogonality properties of different terms in a µ : Now, we test the orthogonality of the radiation reaction term a µ . We just write the expression of a µ from Eqn. (11) in the following way : It is already known that the part Du µ dτ is perpendicular to the four-velocity u µ [9] i.e.
Next, we test the orthogonality with the electromagnetic radiation reaction term 2q 2 3m (g µ ν + u µ uν) D 2 u ν dτ 2 , using the relation u α uα = −1 : The orthogonality of the Lorentz force term q m F µ ν u ν results from the antisymmetry of the field strength tensor under the exchange of the space-time indices. Then, we test the orthogonality of the tail term 2q 2 m f µν T ail uν . The f µν T ail in the tail term is the 'tail integral' given by [10,15]: +λ ′′ (z(τ ), z(τ ′′ ))u λ ′′ dτ ′′ .
Where G µ +λ is the retarded Green's function associated with the vector potential of the electromagnetic field. Hence, con-tracting the term f µν T ail uν with uµ we have: f µν T ail uν uµ = uν uµ +λ ′′ (z(τ ), z(τ ′′ ))u λ ′′ dτ ′′ = uν uµ (82) In this case also, in the two parts on the RHS of the above Eqn.(82), the indices µ and ν are repeated and in the similar way, as can be done for the previous case of Lorentz force, in this case also we interchange the indices (µ ↔ ν) for the first term on the RHS of Eqn.(82) and obtain : Hence, the overall radiation reaction term a µ is orthogonal to the four-velocity uµ : a µ uµ = 0. In the next sub-section we shall show the utility of this orthogonality property.
B. Utility of the Orthogonality Property of the reaction a µ with the four-velocity : As we have the orthogonality property of the overall radiation reaction, we can use it to have constraints or relations between different coefficients and terms present in the reaction. For doing this, we contract the radiation reaction a µ with unperturbed four-velocity uµ in the Eqn. (15). Thus we get : We have already shown that F µν uν uµ = 0 and f µν T ail uµuν = 0 in the previous sub-section. The similar results are also valid for the perturbed external lorentz force and the perturbed tail term : F ′µν uν uµ = 0 and f ′µν T ail uµuν = 0. Using those results and 1 4 (4)F 2 ) So, we have shown that for electromagnetic field the Ricciscalar R is zero (0). Now, we substitute R = 0 in the equation 88 and get : Or, in mixed-tensorial form, The term containing Ricci-tensors, appearing in the equation of motion of a charged particle in curved space-time, is q 2 3m (R µ λ u λ + R ν λ u λ u µ uν ). We first evaluate the first part R µ λ u λ :