Orbital and epicyclic frequencies around neutron and strange stars in $R^2$ gravity

According to various models, the orbital and the epicyclic frequencies of particles moving on a circular orbit around compact objects are related to the quasi-periodic oscillations observed in the X-ray flux of some pulsars or black hole candidates. It is expected that they originate from the inner edge of the accretion discs, deep into the gravitational field of the compact objects. Considering the planned new generation X-ray timing observatories with large collective areas, the quasi-periodic oscillations might be an excellent tool for testing gravity in strong field regime and respectively alternative gravitational theories. We examine the orbital and the epicyclic frequencies of a particle moving on a circular orbit around neutron or strange stars in $R^2$ gravity. The case of slow rotation is considered too. The $R^2$ gravity results are compared to the General Relativistic case. We comment the deviations from General Relativity, as well as the deviations due to rotation in both theories.


Introduction
In the last few years the interest in alternative theories of gravity was significantly increased. Major role for this has the experimental confirmation of the accelerated expansion of the universe and the fact that it does not fit in the predictions of General Relativity (GR) without the introduction of the so-called dark energy, a matter with exotic properties which interacts with the a kstaykov@phys.uni-sofia.bg b daniela.doneva@uni-tuebingen.de c yazad@phys.uni-sofia.bg visible matter only gravitationally. A class of viable alternative theories of gravity are the so-called f (R) theories. In this case we exchange the Lagrangian of the General Relativistic Einstein-Hilbet action, namely the Ricci scalar, with a more general one -a function of the Ricci scalar, hence f (R) theories [1][2][3][4]. These theories allow us to resolve the problem with the accelerating universe without the necessity of introducing such an exotic matter. However, the theory should be tested not only on cosmological but on astrophysical scales too. The predictions, concerning astrophysical effects, should be very close to the GR ones in the weak field regime, but it is expected to deviate for strong fields. As a source of strong gravitational fields neutron stars (NS) and black holes (BH) are natural laboratories for testing alternative theories of gravity. In this work we are concentrating our efforts on the former. Future X-ray timing observatories with large collecting area like SKA [5], NICER [6], LOFT [7], and AXTAR [8] may give us a chance to test gravity in strong field regime. A promising way to do that around compact objects are the so-called quasi-periodic oscillations (QPOs). QPOs are Hz to kHz oscillations in the X-ray flux of pulsars or black hole candidates. The latter QPOs are supposed to occur only in the presence of strong gravitational fields, with their origin in the inner edge of the accretion disc, deep into the gravitational field of the star. However, the exact source of these oscillations is unknown so there are different models for their explanation. In the beat frequency models, some connection between the orbital frequency and the spin frequency of the central compact object is suggested. These models require some azimuthally non-uniform structure co-rotating with the central object. In the relativistic resonance models a resonance may occur at a particular radii of the disc at which the orbital and the 2 epicyclic frequencies have integer ratios. It is suggested that some form of resonance may occur due to the interaction between the accretion disc and the central body too. The preferred radii models suggest that mechanism, that chooses some radii, exists. The relativistic precession models, investigated in some resent papers [9][10][11][12][13], are based on the assumption that the QPOs are directly related to the orbital (Kepler) Ω p and to the epicyclic frequencies. In the above-mentioned models all these frequencies or some of them are involved. The radial ν r and the vertical ν θ epicyclic frequencies occur if a particle on a stable circular orbit is perturbed. It will start to oscillate in a radial and in a vertical direction with some stable frequencies. This frequencies are the radial and the vertical epicyclic ones. A review of these models can be found in [14]. In some papers QPOs models concerning oscillations of the accretion disk itself are examined too [15][16][17].
The epicyclic frequencies of rapidly rotating strange stars are examined in the resent paper [18]. Some of the above mentioned models were examined in alternative theories of gravity through the years [12,[19][20][21]. The observed deviations from GR in these papers were quite low, except the case of rapidly rotating neutron stars in scalar-tensor theories [21].
The structure of this paper is as follow: In Section II we present the basic steps for deriving the radius of the innermost stable circular orbit (ISCO), the orbital frequency and the radial and the vertical epicyclic frequencies. In Section III we present and discus the results for neutron and strange stars. The paper ends with conclusions.

The epicyclic frequencies of a rotating neutron star
In this section we briefly present the basic steps in the derivation of the equations for the radius of the innermost stable circular orbit (ISCO), the equations for the radial and for the vertical epicyclic frequencies and for the orbital frequency [12,[22][23][24]. The equations describing stable stationary neutron star models and additional mathematical details concerning R 2 gravity can be found in [25][26][27].
We are considering a stationary and axisymmetric spacetime with a metric ds 2 = g tt dt 2 + g rr dr 2 + g θθ dθ 2 + 2g tϕ dtdϕ + g ϕϕ dϕ 2 , where all the metric functions depend only on the coordinates r and θ. The massive particles subject to the gravitational force only move on timelike geodesics of the metric (1). The stationary and axial Killing symmetries of metric, generated by the Killing vectors ∂ ∂t and ∂ ∂ϕ , give rise to two constants of motion, namely E = −u t and L = u ϕ . The first one corresponds to the energy per unit mass and the second one to the angular momentum along the axis of symmetry, and u µ =ẋ µ = dx µ /dτ is the four-velocity of the particle. It is not difficult one to show that the two conservation laws can be casted in the form where we defined g2 = g 2 tϕ − g tt g ϕϕ for simplicity. Here t denotes the coordinate time, and τ the proper time. From the normalization condition g µν u µ u ν = −1, we also have with U (r, θ) = g ϕϕ + 2lg tϕ + l 2 g tt g2 (5) and l = L/E being the proper angular momentum.The derivatives in the above equations are with respect to the proper time τ . For θ = π 2 the problem reduces to an effective one dimensional probleṁ with an effective potential The stable circular orbit with a radiusr is determined by the conditions V (r c ) = 0 = V ′ (r c ) and V ′′ (r c ) > 0, where with prime we denote the derivative with respect to r . The condition V ′′ (r c ) = 0 gives the ISCO radius. The angular velocity Ω p of a particle moving on a circular equatorial orbit can be found from the geodesic equations in the following way. We write down the geodesic equations in the form which for the radial coordinate gives Taking into account that the angular velocity is defined by Ω p = u ϕ u t = dϕ dt we obtain from the above equation To derive the epicyclic frequencies we should investigate small perturbations of a stable circular orbit. The perturbations are written in the form where δr(t) and δθ(t) are perturbations to the stable circular orbit with coordinate radiusr in the equatorial plane. The perturbations could be written explicitly in the form δr(t) ∼ e 2πiνr t and δθ(t) ∼ e 2πiν θ t . Substituting (11) in eq. (4) and after some calculations, and a change of the proper time τ with the coordinate one t, we obtain the expressions for the radial and the vertical epicyclic frequencies: For the case of static neutron stars the orbital frequency and the vertical epicyclic frequency coincide, i.e. ν θ = ν p , for f = 0, where f is the rotational frequency of the star (f = Ω 2π ). At the ISCO the square of the radial epicyclic frequency is equal to zero, and for smaller radius it is negative, which shows a radial instabilities for orbits with radius smaller than the ISCO.
For simple accretion disc models the inner edge of the disc is defined by the ISCO. For models with smaller masses the ISCO is in the interior of the star and for massive ones the radius of the ISCO is bigger than the radius of the star itself. The inner edge of the accretion disc, therefore, can reach down to the surface of the star in the former case and to the ISCO in the latter one.

Numerical results
We investigate what are the changes in the radius of the ISCO and in the orbital and epicyclic frequencies in f (R) gravity with Lagrangian f (R) = R + aR 2 , the socalled R 2 gravity. The results are compared to pure GR. Deviations due to rotation, in slow rotation approximation, are examined too. We consider two hadronic equations of state (EOS) and a quark one. For the hadronic EOS we are using peacewise polytropic approximation [28]. The ones we choose are APR4 and MS1. APR4 has maximal mass not much bigger than the observational limit of two solar masses (medium stiffness), and MS1 is a stiffer one with higher masses and radii. We examined a soft EOS too but the results are qualitatively the same, and because of that, graphs are not presented. The quark EOS has the analytical form where the constants b and ρ 0 are taken from [29] for EOS SQS B60. This EOS leads to maximal masses slightly bellow two solar masses, but we find it to be a good representative.
In this section we examine static, f = 0 Hz, and slowly rotating, f = 80 Hz and f = 160 Hz, models of neutron and strange stars. First we will discus the deviations from pure GR due to R 2 gravity and as a next step we will examine the differences of the results for different rotational frequencies. The presented results in all figures are up to the maximal mass. For the presented results in R 2 gravity, the maximal deviation from GR is for the maximal adopted value of the parameter a = 10 4 . This number is close to the maximal values allowed by the observations a ∼ 10 5 or in dimensional units -a 5 × 10 11 m 2 [30].
In Fig. 1 sequences of models, representing the radius of the ISCO as a function of the stellar mass are plotted. A wide range of values for the parameter a is examined. In the left panel the results for EOS APR4 are plotted, in the middle one for EOS MS1, and in the right one for EOS SQS B60. If the ISCO is in the interior of the star, that happens for models with low masses, we are plotting the radius of the star instead. The points in the different sequences of models, where the radius of the ISCO gets equal to the radius of the star are marked with an asterisk. For all examined EOS the results are qualitatively the same. For models with maximal masses the largest deviation is for the maximal values of the parameter a. It is around 10 % and decreases with the decrease of a. In the limiting case of a → 0 the solutions converge to the GR ones.
At this point a comment concerning the presented results should be done. In Fig. 1, in the present paper, the following behaviour for the models for which the radius of the ISCO gets bigger than the radius of the star can be observed. Let us start the discussion with the models marked with an asterisk. For small values of the parameter a they have smaller radius and lower mass than the GR case. The a = 10 models have parameters which are close to the GR ones and for high values of a the radii and the masses are higher compared to the GR case. When the ISCO is outside the star its radius is bigger compared to the GR case for models with equal masses and small values of a. For a = 10 the re-4 sults are comparable to the GR ones. For bigger values of a the radius of the ISCO is smaller than the GR one. The explanation of this behaviour we find in the nonmonotonous behaviour of some stellar parameters as a function of a. This is demonstrated, for example, in Fig.  3 in [25] for the maximal mass. It was shown that for small values of a the mass decreases compared the GR one. At some value of the parameter a there is a turning point and the mass starts to increase, reaching maximal masses higher than the GR one. Our investigations show that this is the case for neutron stars as well as for strange stars. Such behavior can be observed also for other stellar parameters, although in the different cases the minimum will be shifted to some other values of the parameter a. If one, for example, calculates the mass and the radius for models with fixed central energy density in GR and in R 2 gravity, the same pattern can be found. It is not hard to see that there is a similar dependence for the models marked with an asterisk as well as for the case when the ISCO is outside the star. This leads us to the conclusion that this behaviour is characteristic one for the theory. Naturally, similar patterns can be found in the other graphs presented below.
In Fig. 2 we plot the orbital frequency ν p = Ω p /2π, in kHz, at the ISCO as a function of the stellar mass. If the ISCO is in the interior of the star, we calculated ν p at the surface of the star. For all EOS we see qualitatively the same results. The deviations from GR for models with a = 10 4 is in average around 15 % and it does not change significantly with the mass. For models with small value of a the deviation from GR, when the radius of the ISCO get bigger than the radius of the star, is close to the maximal one, but it rapidly decreases with the increase of the mass, and converge to GR for maximal masses.
In Fig. 3 we plot the maximal value of the radial epicyclic frequency, in kHz, as a function of the stellar mass. If the ISCO is in the interior of the star, we take the maximal value on the surface or outside the star. In the examined interval of masses the frequency monotonically decreases. For higher values of the parameter a the maximal radial epicyclic frequency is higher, compared to the GR one. For small values of a the frequencies are lower than the GR ones. For smaller masses there are higher deviations, but for maximal masses the plots converge to GR.
The deviations due to rotation are qualitatively and quantitatively the same for the hadronic and for the quark EOS. This, combined with the small magnitude of the deviations is the reason why we chose to include only the favoured by the observations hadronc EOS in the following graphs.
In Fig. 4 we plot the difference between the orbital frequency ν p and the vertical epicyclic frequency ν θ , i.e. the nodal precession frequency ν n , as a function of the mass. In the left panel we plot the results for EOS APR4 and in the right panel for EOS MS1. The plots are for three different rotational rates: f = 0 Hz (the static case), f = 80 Hz, and f = 160 Hz. In black continuous lines are the result for GR and in colour dashed line -the results for a = 0.3 and a = 10 4 . As we said before, ν p and ν θ coincide in the static case. With the increase of the rotational frequency of the star the differences between the two frequencies increase.The difference is calculated on the surface of the star if the ISCO is inside the star, and at the ISCO in the opposite case. The transition models are marked with asterisks, in consistence with the previous graphs. The nodal precession frequency tends to decrease with the increase of the mass for models having ISCO in the interior of the star. It increases with the increase of the mass for models having ISCO outside the star. In the case of R 2 gravity we plot the results for a = 10 4 and a = 0.3. The former value of the parameter gives the maximal deviation from GR, and the latter we choose because of the interesting behaviour demonstrated in the graphs above. If the ISCO is in the interior of the star, the results are the same for all rotational rates. In the opposite case, when the ISCO is outside the star, the behaviour is the following. For a = 0.3, ν n is slightly higher than the GR case, and rapidly converge to GR with the increase of the mass. For a = 10 4 , ν n has lower values than the GR one. The deviation from GR is more or less constant with the increase of the mass, and it is about 4 %.
In Figs. 5, 6, 7 we plot the radius of the ISCO, the orbital frequency and the maximal radial epicyclic frequency as a function of the mass in the case of GR and the maximal deviation in R 2 gravity (a = 10 4 ) for the aforementioned values of f . The rotation have the same effect for all these quantities in GR and in R 2 gravity. The reason why we are skipping the interestingly behaving case of a = 0.3 is the following. The changes of the ISCO and the examined frequencies for different rotational rates causes the graphs for a = 0.3 and a = 1 to overlap with the GR ones. We find that case to be quite uninformative so only the maximal deviation is presented, namely a = 10 4 . Due to rotation the radius of the ISCO decreases with about 4 %, and the frequencies increase with the same magnitude. The percentage deviations are the same for GR and for R 2 gravity.

Conclusions
In this paper we examined the orbital and the epicyclic frequencies of particles moving in circular orbits around neutron and strange stars. Various models relate all these frequencies or some of them to QPOs. Taking into consideration the new generation observatories for timing of X-ray pulsars with large collective area which are under construction now, QPOs are expected to have major role in testing strong field regime of gravity.
R 2 gravity is a viable alternative to GR. We examined how changing the parameter of the theory effects the radius of the ISCO, the orbital and the epicyclic frequencies. Most of the observed neutron stars are slowly rotating, so we concentrated our efforts on the static and slowly rotating solutions. We investigated the de-viations of the slow rotating solutions from the static case too.
We concerned two hadronic and one quark equations of state and for both cases we observed qualitatively the same behaviour with the change of the parameter of the theory. The largest deviations from GR occurs for the maximal adopted value of the parameter a. The radius of the ISCO decreases with about 10% and the orbital and the maximal radial epicyclic frequencies increases with 15 -20 %. However, the effect due to slow rotation is expectedly small. For the highest examined rotational frequency (f = 160 Hz) the deviation from the static case is about 4% for all examined quantities.
The deviations for slowly rotating stars, up to 160 Hz, are much smaller than the deviations due to the R 2 gravity. Because of that static theoretical and numerical       models can be used for comparison with observational data even for slowly rotating stars.
The frequencies examined in this paper provide us with an opportunity to test GR in strong field regime, as well as alternative theories. Results from the numerical models could be compared to the observational data expected from the new generation X-Ray timing observatories. However, even if the observational data do not provide a concrete information about the correct gravitational theory, it will provide us with an unique opportunity to restrict the wide range of possible values for the free parameter of the theory.