# Innermost stable circular orbits and epicyclic frequencies around a magnetized neutron star/Orbitas circulares marginalmente estables y frecuencias epiciclicas alrededor de una estrella de neutrones magnetizada/Orbitas circulares marginalmente estave is e frequencias epicicloidais nas proximidades de uma estrela de neutrons magnetizada.

IntroductionStellar models are generally based on the Newtonian universal law of gravitation. However, considering the size, mass or density, there are five classes of stellar configurations where one can recognize significant deviations from the Newtonian theory, namely, white dwarfs, neutron stars, black holes, supermassive stars and relativistic star clusters (Misner et al. 1973). In the case of magnetized objects, such as white dwarfs (~[10.sup.5] T) or neutron stars (~ [10.sup.10] T), the Newtonian theory not only fails in describing the gravitational field generated by the matter distribution, but also in accounting for the corrections from the energy stored in the electromagnetic fields. Despite this fact and due to the sheer complexity of an in-all-detail calculation of the gravitational and electromagnetic fields induced by these astrophysical objects, one usually appeals to approaches based on post-Newtonian corrections (Aliev & Ozdemir 2002, Preti 2004), which may or may not be enough in order to provide a complete and accurate description of the space-time around magnetized astrophysical objects. In particular, these approaches consider, e.g., that the electromagnetic field is weak compared to the gravitational one and therefore, the former does not affect the space-time geometry (Aliev & Ozdemir 2002, Preti 2004, Mirza 2005, Bakala et al. 2010, 2012). That is, it is assumed that the electromagnetic field does not contribute to the space-time curvature, but that the curvature itself may affect the electromagnetic field. Based on this approximation, the space-time around a stellar source is obtained from a simple solution of the Einstein's field equations, such as Schwarzschild's or Kerr's solution, superimposed with a dipolar magnetic field (Aliev & Ozdemir 2002, Preti 2004, Mirza 2005, Bakala et al. 2010, 2012). Although, to some extend this model may be reliable for weakly magnetized astrophysical sources, it is well known that in presence of strong magnetic fields, non-negligible contributions to the space-time curvature are expected, and consequently on the physical parameters that describe the physics in the neighbourhood of these objects.

In particular, one expects contributions to the radius of the innermost stable circular orbit (ISCO) (Sanabria-Gomez et al. 2010), the Keplerian frequency, frame-dragging frequency, precession and oscillation frequencies of the radial and vertical motions of test particles (Stella & Vietri 1999, Bakala et al. 2010, 2012), and perhaps to other physical properties such as the angular momentum of the emitted radiation (Tamburini et al. 2011), which could reveal some properties of accretion disks and therefore of the compact object (Bocquet et al. 1995, Konno et al. 1999, Broderick et al. 2000, Cardall et al. 2001). In other words, to construct a more realistic theoretical description that includes purely relativistic effects such as the modification of the gravitational interaction by electromagnetic fields, it is necessary to use a complete solution of the full Einstein-Maxwell field equations that takes into account all the possible characteristics of the compact object. In this paper, we use the six parametric solution derived by Pachon et al. (2006) (hereafter PRS solution), which provides an adequate and accurate description of the exterior field of a rotating magnetized neutron star (Pachon et al. 2006, 2012), to calculate the radius of the ISCO, the Keplerian, frame-dragging, precession and oscillation frequencies of neutral test-particles orbiting the equatorial plane of the star. The main purpose of this paper is to show the influence of the magnetic field on these particular quantities.

The paper is organized as follows. In section Description of the space-time around the source, we briefly describe the physical properties of the PRS solution, the general formulae to calculate the parameters that characterize the dynamics around the star are presented in section Characterization of the dynamics around the source. Sections Influence of the dipolar magnetic field in the ISCO radius and Keplerian and epicyclic frequencies are devoted to the study of the influence of the magnetic field on the ISCO radius and on the Keplerian and epicyclic frequencies, respectively. The effect of the magnetic field on the energy E and the angular momentum L are outlined in section Energy and angular momentum. Finally, the conclusions of this paper are given in the Concluding remarks.

Description of the space-time around the source

According to Papapetrou (1953), the metric element d [s.sup.2] around a rotating object with stationary and axially symmetric fields can be cast as

d[s.sup.2] = -f [(dt - [omega]d[phi]).sup.2] + [f.sup.-1] [[e.sup.2[gamma]] (d [[rho].sup.2] + d[z.sup.2]) + [[rho].sup.2]d[[phi].sup.2]],

where f, [gamma] and [omega] are functions of the quasiclyndrical Weyl-Papapetrou coordinates ([rho], z). The non-zero components of metric tensor, which are related to the metric functions f, [omega] and [gamma], are

[g.sub.[phi][phi]] [[rho].sup.2]/f([rho],z) - f([rho],z)[omega][([rho],z).sup.2], [g.sub.tt] = -f ([rho],z), (1)

and

[g.sub.t[phi]] = f([rho],z)[omega]([rho],z), [g.sub.zz] = [g.sub.[rho][rho]] [e.sup.2[gamma]([rho],z)]/f([rho],z) = 1/[g.sup.zz] = 1/[g.sup.[rho][rho]]. (2)

By using the Ernst procedure and the line element in equation (1), it is possible to rewrite the Einstein-Maxwell equations in terms of two complex potentials E([rho], z) and [PHI]([rho], z) [see Ernst (1968) for details]

(ReE + [[absolute value of ([PHI]].sup.2])) [[nabla].sup.2]E = ([nabla]E + 2[PHI]*[nabla][PHI]) x [nabla]E, (ReE + [[absolute value of ([PHI]].sup.2])) [[nabla].sup.2][PHI] = ([nabla]E + 2[PHI]*[nabla][PHI]) x [nabla][PHI]), (3)

where * stands for complex conjugation. The above system of equations can be solved by means of the Sibgatullin integral method (Sibgatullin 1991, Manko & Sibgatullin 1993), according to which the Ernst potentials can be expressed as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],

where [xi] = z + i[rho][sigma] and e(z) = E(z, [rho] = 0) and f (z) = [PHI](z, [rho] = 0) are the Ernst potentials on the symmetry axis. As shown below, these potentials contain all information about the multipolar structure of the astrophysical source [see also Pachon & Sanabria-Gomez (2006) for a discussion on the symmetries of these potentials]. The auxiliary unknown function [mu]([sigma]) must satisfy the integral and normalization conditions

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

Here [eta] = z + i[rho][tau] and [sigma],[tau] [member of] [-1,1], [??]([eta]) = [e.sup.*]([[eta].sup.*]), [??]([eta]) = [f.sup.*]([[eta].sup.*])

For the PRS solution (Pachon et al. 2006), the Ernst potentials were chosen as:

e(z)=[z.sup.3] - [z.sup.2] (m + ia) - kz + is/[z.sup.3] + [z.sup.2] (m - ia) - kz + is, f(z) = q[z.sup.2] + i[mu]z/[z.sup.3] + [z.sup.2] (m - ia) - kz + is. (5)

The electromagnetic and gravitational multipole moments of the source were calculated by using the Hoenselaers & Perjes method (Hoenselaers & Perjes 1990) and are given by (Pachon et al. 2006):

[M.sub.0] = m, [M.sub.2] = (k-[a.sup.2])m, ... [S.sub.1] = a, [S.sub.3] = - m([a.sup.3] - 2ak + s), ... (6)

[Q.sub.0] = q, [Q.sub.2] = -[a.sup.2]q - a[mu] + kq, ... [B.sub.1] = [mu] + aq, ... [B.sub.3] = -a[mu] + [mu]k - [a.sup.3]q + 2akq - qs, ... (7)

where the [M.sub.i]s denote the moments related to the mass distribution and [S.sub.i]s to the current induced by the rotation. Besides, the [Q.sub.i]s are the multipoles related to the electric charge distribution and the [B.sub.i]s to the magnetic properties. In the previous expressions, m corresponds to the total mass, a to the total angular moment per unit mass (a = J/[M.sub.0], being J the total angular moment), q to the total electric charge. In our analysis, we set the electric charge parameter q to zero because, as it is case of neutron stars, most of the astrophysical objects are electrically neutral. The parameters k, s and [mu] are related to the mass quadrupole moment, the current octupole, and the magnetic dipole, respectively.

Using Eqs. (4) and (5), the Ernst potentials obtained by Pachon et al. (2006) are

E = [A+B/A-B], [PHI] = C/A-B, (8)

which leads to the following metric functions:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

[omega] = [Im[(A+B)[bar.H] - ([bar.A] + [bar.B]) G-C[bar.I]]/A[bar.A] - B[bar.B] + C[bar.C]]. (10)

The analytic expressions for the functions A, B, C, H, G, K, and I can be found in the original reference (Pachon et al. 2006) or in Appendix of Pachon et al. (2012). A Mathematica 8.0 script with the numerical implementation of the solution can be found at http://gfam.udea.edu.co/~lpachon/scripts/nstars.

Characterization of the dynamics around the source

In the framework of general relativity, the dynam-ics of a particle may be analyzed via the Lagrangian formalism as follows. Let us consider a particle of rest mass [m.sub.0] = 1 moving in a space-time characterized by the metric tensor [g.sub.[mu]v], thus the Lagrangian of the particle is given by

L = 1/2 [g.sub.[mu]v][[??].sup.[mu]][[??].sup.v], (11)

where the dot denotes differentiation with respect to the proper time [tau], [x.sup.[mu]]([tau]) are the coordinates. Since the fields are stationary and axisymmetric, there are two constants of motion related to the time coordinate t and azimuthal coordinate [phi] [see Ryan (1995) for details], these are given by:

E = -[partial derivative]L/[partial derivative][??] = [g.sub.tt] (dt/d[TAU]) - [g.sub.t[phi]] (d[phi]/d[tau]), (12)

L = [partial derivative]L/[partial derivative][??] = [g.sub.t[phi]] (dt/d[TAU]) - [g.sub.[phi][phi]] (d[phi]/d[tau]), (13)

where E and L are the energy and the canonical angular momentum per unit mass, respectively. For real massive particle, the four velocity is a time-like vector with normalization [g.sub.[mu]v] [u.sup.[mu]] [u.sup.v] = -1. If the motion takes place in the equatorial plane of the source z = 0, this normalization condition leads to

[g.sub.[rho][rho]] [(d[rho]/d[tau]).sup.2] = -1 + [g.sub.tt] [(dt/d[tau]).sup.2] + [g.sub.[phi][phi]] [(d[phi]/d[tau]).sup.2] + [2g.sup.t[phi]] (dt/d[tau]) (d[rho]/d[tau]). (14)

From equation (14), one can identify an effective potential that governs the geodesic motion in the equatorial plane [see e.g. Bardeen et al. (1972)]

[v.sub.eff]([rho]) = 1 - [[E.sup.2][g.sub.[phi][phi]] + 2EL[g.sub.t[phi]] + [L.sup.2][g.sub.tt]/[g.sup.2.sub.t[phi]] - [g.sub.[phi][phi]][g.sub.tt], (15)

For circular orbits, the energy and the angular momentum per unit mass, are determined by the conditions [V.sub.eff]([rho]) = 0 and d[V.sub.eff]/d[rho] = 0. Based on these conditions, one can obtain expressions for the energy and the angular momentum of a particle that moves in a circular orbit around the star [see e.g. Stute & Camenzind (2002)], namely,

E = [square root of f]/[square root of 1 - [f.sup.2] [chi square]/[[rho].sup.2]], L = E([omega] + [chi]), (16)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (17)

where the colon stands for a partial derivative respect the lower index. The radius of innermost stable circular orbit (ISCO's radius) is determined by solving numerically for [rho] the equation

[d.sup.2][V.sub.eff]/d[[rho].sup.2] = 0,

which arising from the marginal stability condition. This condition together with the equations (16) and (17) can be write in terms of the metric functions as [see Stute & Camenzind (2002)]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (18)

As is usual in the literature the physical ISCO radius reported here corresponds to evaluation of [square root of [g.sub.[phi][phi]]] at the root of equation (18). This equation is solved for fixed total mass of the star M, the dimensionless spin parameter j = J/[M.sup.2] (being J the angular momentum), the quadrupole moment [M.sub.2] and the magnetic dipolar moment [mu] [see Table VI of Pappas & Apostolatos (2012)].

Keplerian, oscillation and precession frequencies

The Keplerian frequency [[OMEGA].sub.K] at the ISCO can be obtained from the equation of motion of the radial coordinate [rho]. This equation is easily obtained by using the Lagrangian (11),

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (19)

By imposing the conditions of circular orbit or constant orbital radius, d[rho]/d[tau] = 0 and [d.sup.2][rho]/d[[tau].sup.2] = 0 and taking into account that d[phi]/d[tau] = [[OMEGA].sub.K]dt/d[tau], one gets [see e.g. Ryan (1995)]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (20)

where "+" and "-" denotes the Keplerian frequency for corotating and counter-rotating orbits, respectively.

The epicyclic frequencies are related to the oscillations frequencies of the periastron and orbital plane of a circular orbit when we apply a slightly radial and vertical perturbations to it. According to Ryan (1995), the radial and vertical epicyclic frequencies are given by the expression

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII],

where [alpha] = {[rho],z}. The periastron [v.sup.p.sub.[rho]] and the nodal [v.sup.p.sub.z] frequencies are defined by:

[v.sup.p.sub.[rho]] = [[OMEGA].sub.k]/2[pi] - [v.sub.[alpha]], (21)

which are the ones with an observational interest (Stella & Vietri 1999). The radial oscillation frequency vanishes at the ISCO radius and therefore, the radial precession frequency equals to the Keplerian frequency.

Finally, the frame dragging precession frequency or Lense-Thirring frequency [V.sub.LT] is given by [see e.g. Ryan (1995)]

[v.sub.LT] = - 1/2[pi] [g.sub.t[phi]]/[g.sub.[phi][phi]]. (22)

[v.sub.LT] is related to purely relativistic effects only (Misner et al. 1973). In this phenomenon, the astrophysical source drags the test particle into the direction of its rotation angular velocity. In absence of electromagnetic contributions, as in the case of the Kerr solution, the frame dragging comes from the non-vanishing angular momentum of the source. In the presence of electromagnetic contributions in non-rotating sources, it was shown by Herrera et al. (2006) that the non-zero circulation of the Poynting vector is able to induced frame dragging, in which case it is electromagnetically induced. In the cases discussed below, due to the presence of fast rotations and magnetic fields, the frame dragging will be induced by a combination of these two processes.

Influence of the dipolar magnetic field in the ISCO radius

As discussed in the Introduction, contributions from the energy stored in the electromagnetic fields come via equivalence between matter and energy, E = m[c.sup.2]. For the particular case of a magnetar (~[10.sup.10]T), the electromagnetic energy density is around 4 x [10.sup.25] J/[m.sup.3], with an E/[c.sup.2] mass density [10.sup.4] times larger than that of lead. Hence, the relevant physical quantities should depend upon the magnetic field, and in particular on the magnetic dipole moment [mu].

Figure 1 shows the ISCO radius as a function of the parameter [mu] for four particular realistic numerical solutions for rotating neutron stars models derived by Pappas & Apostolatos (2012). The models used coincide with the models 17,18,19 and 20 of Table VI in that reference and correspond to results for the Equation of State L (see Table 1). The lowest multipole moments of the PRS solution, namely, mass, angular moment and mass quadrupole were fixed to t he n umerical o nes o btained by P appas & Apostolatos (2012) (see Table 1). Since the main objective here is to analyze the influence of the magnetic dipole, the current octupole parameter s was set to zero. Note that this does not mean that the current octupole moment vanishes (see in Eq. (6). The parameter [mu] was varied between 0 and 6.4 [km.sup.2], which corresponds to magnetic fields dipoles around 0 and 6.3 x [10.sup.31] [Am.sup.2], respectively. In all these cases, the ISCO radius decreases for increasing [mu], this can be understood as a result of the dragging of inertial frames induced by the presence of the magnetic dipole, we elaborate more on this below.

Keplerian and epicyclic frequencies

Some of the predictions from General Relativ-ity (RG) such as the dragging of inertial systems (frame-dragging or Lense-Thirring effect) (Everitt et al. 2011), the geodesic precession (geodesic effect or de Sitter precession) (Everitt et al. 2011) or the analysis of the periastron precession of the orbits (Lucchesi & Peron 2010) have been experimental verified. However, since they have observed in the vicinity of the Earth, these observations represent experimental support to RG only in the weak field limit. Thus, it is fair to say that the RG has not been checked in strong field limit [see e.g. Psaltis (2008)]. In this sense, the study of compact objects such as black holes, neutron stars, magnetars, etc., is certainly a topic of great interest, mainly, because these objects could be used as remote laboratories to test fundamental physics, in particular, the validity of general relativity in the strong field limit yet (Stella & Vietri 1999, Pachon et al. 2012).

The relevance of studying the Keplerian, epicyclic and Lense-Thirring frequencies relies on the fact that they are usually used to explain the quasiperiodic oscilation phenomena present in some Low Mass X-ray Binaries (LMXRBs). In this kind of systems a compact object accretes from another one (which is usually a normal star) and the X-ray emissions could be related to the relativistic motion of the accreted matter e.g. rotation, oscillation and precession. Stella & Vietri (1999) showed that in the slow rotation regime, the periastron precession and the Keplerian frequencies could be related to the phenomena of the kHz quasiperiodic oscillations (QPOs) observed in many accreting neutron stars in LMXRBs (Stella & Vietri 1999). This model is known as the Relativistic Precession Model (RPM) (Stella & Vietri 1999) and identifies the lower and higher QPOs frequencies with the periastron precession and the Keplerian frequencies respectively. The RPM model has been used to predict the values of the mass and angular momentum of the neutron stars in this kind of systems [see e.g. Stella & Vietri (1999) and Stella et al. (1999), for details].

In the framework of the RPM, it is assumed that the motion of the accretion disk is determined by the gravitational field alone and thereby, it is normal to assume that the exterior gravitational field of the neutron star is well described by the Kerr metric. However, most of the neutron stars have (i) quadrupole deformations that significantly differ from Kerr's quadrupole deformation (Laarakkers & Poisson 1999) and (ii) a strong magnetic field. The influence of the non-Kerr deformation was discussed, e.g., by Johannsen & Psaltis (2010) and Pachon et al. (2012). Based on observational data, it was found that non-Kerr deformations dramatically affect, e.g., predictions on the mass of the observed sources. Below, we consider the influence of the magnetic field and find that it introduces non-negligible corrections in the precession frequencies and therefore, on the predictions made by the Kerr-based RPM model.

In the previous section, it was shown that the magnetic field affects the value of the ISCO radius of a neutral test particle that moves around the exterior of a magnetized neutron star. Below, it is shown that other physicals quantities such as the Keplerian and epicyclic frequencies, which are used to describe the physics of accretion disk, are also affected. The influence on the Lense-Thirring frequency is also discussed.

Figure 2 shows the influence of the magnetic field in the Keplerian [Fig. 2(a)], nodal precession [Fig. 2(b)] and Lense-Thirring [Fig. 2(c)] frequencies. In all cases, the frequencies are plotted as a function of the parameter [mu] while fixing the others free param-eters (mass, angular moment and mass quadrupole) according to the realistic numerical solutions in Table 1. As it is expected for shorter ISCO radius (see previous section), all frequencies increase with increasing [mu]. The changes in the Lense-Thirring fre-quency come from the electromagnetic contribution discussed above.

Energy and angular momentum

In order to understand the results presented above, we consider that it is illustrative to consider first the effect that the magnetic field has on the energy E and the angular momentum L needed to described marginally stable circular orbits [see equations (16) and (17)]. In doing so, we fix the total mass [M.sub.0], the spin parameter j and the quadrupole moment [M.sub.2], according to the values in Table 1. Figures 3(a) and 3(b) show E and L at the ISCO as a function of the dipolar moment [mu]. By contrast to the case of the Keplerian frequency, an increase of [mu] decreases the value of the energy and the angular momentum needed to find an ISCO. Complementarily, in figures 4(a) and 4(b), E and L are depicted as a function of [rho] for various values of the dipolar moment. Since the energy and the angular momentum associated to the ISCO correspond to the minima of the curves E([rho]) and L([rho]) (indicated by triangles), figure shows that an increase of the magnetic dipole moment induces a decrease of the ISCO radius (see Fig. 1 above) and simultaneously a decrease of the energy and the angular momentum.

At a first sight, for an increasing magnetic field, it may seem conspicuous that the angular momentum, of co-rotating test particles, decreases [see Fig. 4(b)] whereas the Keplerian frequency increases [see Fig. 2(a)]. However, this same opposite trend is already present in the dynamics around a Kerr source when the angular momentum of the source is increased [cf. equations (2.13) for the angular momentum and (2.16) for the Keplerian frequency in Bardeen et al. (1972)]. In Kerr's case, the co-rotating test particles are dragged toward the source thus inducing a shorter ISCO radius [cf. equation (2.21) in Bardeen et al. (1972)], and since the leading order in the Keplerian frequency goes as ~ 1/[[rho].sup.2], a larger frequency is expected. By contrast, the contra-rotating test particles are "repelled" by the same effect thus resulting in an increase of the ISCO radius and in a decrease of the Keplerian frequency [see also Fig. 2 in Pachon et al. (2012)].

Having in mind the situation in Kerr's case and by noting that the frame dragging can be induced by current multipoles of any order (Herrera et al. 2006), the goal now is to compare the characteristics of the contributions from the angular momentum and the dipole moment to the Keplerian frequency based on the approximate expansions derived by Sanabria-Gomez et al. (2010) and appeal then the general theory of multipole moments to track the contribution of dipole moment to the current multipole moments. If the expression for the Keplerian frequency [Eq. (21)] and for the angular momentum [Eq. (23)] in Sanabria-Gomez et al. (2010) are analyzed in detail, one finds that the signs of the contributions of the angular momentum and the magnetic dipole of the source coincide, this being said, one could argue that the contribution of the dipole moment is related to an enhancement of the current multipole moment of the source. This remark is confirmed by the general multipole expansion discussed by Sotiriou & Apostolatos (2004). In particular, the influence of the dipole moment in the higher current-multipole moments is clear from Eq. (23)-(25) of this reference.

Hence, the role of the dipole moment of the source in the dynamics of neutral test particles is qualitatively analogous to the role of the angular momentum, albeit it induces corrections of higher orders in the multipole expansion of the effective potential. This is a subject that deserves a detailed discussion and will be explained somewhere else.

Conclusion

The influence of weak electromagnetic fields on the dynamics of test particles around astrophysical objects has mainly been studied for the case of orbiting charged particles. Based on a analytic solution of the Einstein-Maxwell field equations, here it is shown that a strong magnetic field, via the energy mass relation, modifies the dynamics of neutral test particles. In particular, it is shown that an intense magnetic field induces corrections in the Keplerian, precession and oscillation frequencies of the radial and vertical motions of the test particles as well as in the dragging of inertial frames.

In particular, it was shown that if the angular momentum and the dipole moment are parallel (note that j and [mu] have the same sign), then the ISCO radius of co-rotating orbits decreases for increasing dipole moment, this leads to an increase of the Keplerian, and precession and oscillation frequencies of the radial and vertical motions frequency. The angular momentum and the energy of the ISCO decrease for increasing dipole moment as a consequence of the dragging of inertial frames.

The kind of geodetical analysis performed here is widely used for instance, in the original RPM model (Stella & Vietri 1999, Stella et al. 1999), and int its subsequent reformulations, of the HF QPOs observed in LMXBs. However, Lin et al. (2011) and Torok et al. (2012) have concluded that although these models of HF QPOs, which neglect in the influence of strong magnetic fields, are qualitatively satisfactory, they do not provide satisfactory fits to the observational data. Hence, in order to improve (i) the level of physical description of these models and (ii) the fit to observational data, a more detailed analysis on the role of the extremely strong magnetic field in the structure of the spacetime is necessary and will be performed in the future.

Acknowledgements

This work was partially supported by Fundacion para la Promocion de la Investigacion y la Tecnologia del Banco la Republica grant number 2879. C.A.V.-T. is supported by Vicerrectoria de Investigaciones (UniValle) grant number 7859. L.A.P. ac-knowledges the financial support by the Comite para el Desarrollo de la Investigacion -CODI- of Universi-dad de Antioquia.

Conflict of interest

The authors declare no conflict of interest.

References

Aliev AN, Ozdemir N (2002) Motion of charged particles around a rotating black hole in a magnetic field. Mon. Not. Roy. Astr. Soc. 336(1):241-248 doi:10.1046/j.1365-8711.2002.05727.x

Bakala P, Urbanec M, Sramkova E, Stuchlik Z, Torok G (2012) On magnetic-field-induced corrections to the orbital and epicyclic frequencies: paper II. Slowly rotating magnetized neutron stars. Class. Quantum Grav. 29(6):065012 doi:10.1088/0264-9381 /29/6/065012

Bakala P, Sramkova E, Stuchlik Z, Torok G (2010) On magnetic-field-induced non-geodesic corrections to relativistic orbital and epicyclic frequencies. Class. Quantum Grav. 27(4):045001 doi:10.1088/0264-9381 /27/4/045001

Bardeen JM, Press WH, Teukolsky SA (1972) Rotating Black Holes: Locally Nonrotating Frames, Energy Extraction, and Scalar Synchrotron Radiation. Astrophys. J. 178:347-370 doi:10.1086/151796

Bocquet M, Bonazzola S, Gourgoulhon E, Novak J (1995) Rotating neutron star models with a magnetic field. Astron Astrophys 301:757

Broderick A, Prakash M, Lattimer JM (2000) The Equation of State of Neutron Star Matter in Strong Magnetic Fields. Astrophys. J. 537:351-367 doi:10.1086/309010

Cardall CY, Prakash M, Lattimer JM (2001) Effects of strong magnetic fields on neutron star structure. Astrophys. J. 554:322-339 doi:10.1086/321370

Ernst FJ (1968) New Formulation of the Axially Symmetric Gravitational Field Problem. II. Phys Rev 168:1415-1417 doi:10.1103/PhysRev.168.1415

Everitt CWF, Debra DB, Parkinson BW, Turneaure JP, Conklin JW, et al. (2011) Gravity Probe B: Final Results of a Space Experiment to Test General Relativity. Phys. Rev. Lett. 106(22):221101 doi:10.1103/PhysRevLett.106.221101

Herrera L, Gonzalez GA, Pachon LA, Rueda JA (2006) Frame dragging, vorticity and electromagnetic fields in axially symmetric stationary spacetimes. Class Quant Grav 23(7):2395-2408 doi:10.1088/0264-9381/23/7/011

Hoenselaers C, Perjes Z (1990) Multipole moments of axisymmetric electrovacuum spacetimes. Class. Quantum Grav. 7(10):1819-1825 doi:10.1088/0264-9381/7/10/012

Johannsen T, Psaltis D (2010) Testing the NoHair Theorem with Observations in the Electromagnetic Spectrum: I. Properties of a Quasi-Kerr Spacetime. Astrophys. J. 716:187-197 doi:doi:10.1088/0004-637X/716/1/187

Konno K, Obata T, Kojima Y (1999) Deformation of relativistic magnetized stars. Astron. Astrophys. 352:211-216

Laarakkers WG, Poisson E (1999) Quadrupole Moments of Rotating Neutron Stars. Astrophys. J. 512(1):282-287 doi:10.1086/306732

Lin YF, Boutelier M, Barret D, Zhang SN (2011) Studying Frequency Relationships of Kilohertz Quasi-periodic Oscillations for 4U 163653 and Sco X-1: Observations Confront Theories. Astrophys. J. 726(74):1-12 doi:10.1088/0004637X/726/2/74

Lucchesi DM, Peron R (2010) Accurate Measurement in the Field of the Earth of the General-Relativistic Precession of the LAGEOS II Pericenter and New Constraints on Non-Newtonian Gravity. Phys. Rev. Lett. 105(23):231103 doi:10.1103/PhysRevLett.105.231103

Manko VS, Sibgatullin NR (1993) Construction of exact solutions of the Einstein-Maxwell equations corresponding to a given behaviour of the Ernst potentials on the symmetry axis. Class. Quantum Grav. 10:1383-1404 doi:10.1088/02649381/10/7/014

Mirza BM (2005) Charged particle dynamics in the field of a slowly rotating compact star. Int. J. Mod. Phys D14:609-620 doi:10.1142/S0218271805006250

Misner C, Thorne K, Wheeler J (1973) Gravitation. W.H. Freeman, San Francisco

Pachon LA, Rueda JA, Sanabria-Gomez JD (2006) Realistic exact solution for the exterior field of a rotating neutron star. Phys. Rev. D 73(10):104038 doi:10.1103/PhysRevD.73.104038

Pachon LA, Rueda JA, Valenzuela-Toledo CA (2012) On the Relativistic Precession and Oscillation Frequencies of Test Particles around Rapidly Rotating Compact Stars. Astrophys. J. 756:82 doi:10.1088/0004-637X/756/1/82

Pachon LA, Sanabria-Gomez JD (2006) Note on reflection symmetry in stationary axisymmetric electrovacuum spacetimes. Class. Quantum Grav. 23:3251-3254 doi:10.1088/0264-9381/23/9/N01

Papapetrou A (1953) Eine rotationssymmetrische Losung in der allgemeinen Relativitatstheorie. Annalen der Physik 447:309-315 doi:10.1002/andp.19534470412

Pappas G, Apostolatos TA (2012) Revising the multipole moments of numerical spacetimes, and its consequences. Phys. Rev. Lett. 108:231104 doi:10.1103/PhysRevLett.108.231104

Preti G (2004) On charged particle orbits in dipole magnetic fields around Schwarzschild black holes. Class. Quantum Grav. 21:3433

Psaltis D (2008) Probes and Tests of Strong-Field Gravity with Observations in the Electromagnetic Spectrum. Living Rev. Relativity 11(9) doi:10.12942/lrr-2008-9

Ryan FD (1995) Gravitational waves from the inspiral of a compact object into a massive, axisymmetric body with arbitrary multipole moments. Phys. Rev. D 52:5707-5718 doi:10.1103/PhysRevD.52.5707

Sanabria-Gomez JD, Hernandez-Pastora JL, Dubeibe FL (2010) Innermost stable circular orbits around magnetized rotating massive stars. Phys. Rev. D 82(12):124014 doi:10.1103/PhysRevD.82.124014

Sibgatullin NR (1991) Oscillations and Waves in Strong Gravitational and Electromagnetic Fields. Springer-Verlag

Sotiriou TP, Apostolatos TA (2004) Corrections and comments on the multipole moments of axisymmetric electrovacuum spacetimes. Class. Quantum Grav. 21:5727-5733 doi:10.1088/02649381/21/24/003

Stella L, Vietri M (1999) kHz Quasiperiodic Oscillations in Low-Mass X-Ray Binaries as Probes of General Relativity in the Strong Field Regime. Phys. Rev. Lett. 82:17-20 doi:10.1103/PhysRevLett.82.17

Stella L, Vietri M, Morsink SM (1999) Correlations in the Quasi-periodic Oscillation Frequencies of Low-Mass X-Ray Binaries and the Relativistic Precession Model. Astrophys. J. Lett. 524:L63-L66 doi:10.1086/312291

Stute M, Camenzind M (2002) Towards a self consistent relativistic model of the exterior gravitational field of rapidly rotating neutron stars. Mon. Not. Roy. Astr. Soc. 336:831-840 doi:10.1046/j.1365-8711.2002.05820.x

Tamburini F, Thide B, Molina-Terriza G, Anzolin G (2011) Twisting of light around rotating black holes. Nature Phys. 7:195-197 doi:10.1038/nphys1907

Torok G, Bakala P, Sramkova E, Stuchlik Z, Urbanec M, et al. (2012) Mass-Angular-momentum Relations Implied by Models of Twin Peak Quasiperiodic Oscillations. Astrophys. J. 760(2):138

Edited by Humberto Rafeiro & Alberto Acosta

Andres F. Gutierrez-Ruiz (1), Leonardo A. Pachon (1), Cesar A. Valenzuela-Toledo (2) ([mail])

(1) Grupo de Fisica Atomica y Molecular, Instituto de Fisica, Facultad de Ciencias Exactas y Naturales, Universidad de Antioquia UdeA; Calle 70 No. 52-21, Medellin, Colombia.

(2) Departamento de Fisica, Universidad del Valle, A.A. 25360, Santiago de Cali, Colombia.

Received: 25-06-2013 Accepted: 13-09-2013 Published online: 23-12-2013

doi: 10.11144/Javeriana.SC19-1.isco

Funding: Fundacion para la Promocion de la Investigacion y la Tecnologia del Banco de la Republica; Vicerrectoria de Investigaciones, Universidad del Valle; Comite para el Desarrollo de la Investigacion--CODI, Universidad de Antioquia.

Electronic supplementary material: N/A

Table 1. Realistic numerical solutions for ratating neutron stars derived by Pappas & Apostolatos (2012). Here, [M.sub.0] is the total mass of the star, j is the dimensionless spin parameter: j= J/[M.sup.2.sub.0] (being J the angular momentum), [M.sub.2] is the quadupole moment and [S.sub.3] is the current octupole moment [see Talbe VI of Pappas & Apostolatos (2012)]. [M.sub.0] [M.sub.2] [S.sup.3] Model [km] j [[km.sup.3]] [[km.sup.4]] M17 4.120 0.588 -51.8 -210.0 M18 4.139 0.635 -62.6 -279.0 M19 4.160 0.682 -74.9 -365.0 M20 4.167 0.700 -79.8 -401.0

Printer friendly Cite/link Email Feedback | |

Title Annotation: | ORIGINAL ARTICLE |
---|---|

Author: | Gutierrez-Ruiz, Andres F.; Pachon, Leonardo A.; Valenzuela-Toledo, Cesar A. |

Publication: | Revista Universitas Scientarum |

Date: | Jan 1, 2014 |

Words: | 5566 |

Previous Article: | Chemical and microbiological comparison of biodeterioration in Colombian heritage constructions/Comparacion quimica y microbiologica del biodeterioro... |

Next Article: | Aves de Quindicocha--Valle de Sibundoy, Putumayo--Colombia: potencial area de conservacion. |