Orbital motion and quasi-periodic oscillations with periastron and Lense–Thirring precession of slowly rotating Einstein–Æther black hole

We study the orbital and oscillatory motion of test particles moving around slowly rotating first and second kinds of Einstein–Æther black holes. In relation to the black hole parameters, we find analytical solutions for the radial profiles of specific energy and specific angular momentum of the equatorial stable circular orbits. The properties of the co-rotating as well as contra-rotating innermost stable circular orbits are analyzed. We examine the radial profiles of the frequencies of latitudinal and radial harmonic oscillations as a function of the black hole mass and dimensionless coupling constants of the theory. The key features of quasi-periodic oscillations of test particles near the stable circular orbits in an equatorial plane of the black hole are discussed. We investigate the positions of resonant radii for high-frequency quasi-periodic oscillations models, namely epicyclic resonance and its variants, relativistic precession and its variants, tidal disruption, as well as warped disc models, considered in the background of slowly rotating first and second kinds of slowly rotating Einstein–Æther black holes. Furthermore, Periastron and Lense–Thirring precessions have been discussed. We demonstrate that the dimensionless coupling parameters of the theory have a strong influence on particle motion around Einstein–Æther black holes.


Introduction
General relativity (GR) is a theory of gravitation that was proposed by Albert Einstein in 1915 [1]. GR has served as the most successful theory of gravitation to explain and understand various mysteries of the astrophysical as well as the cosmological realm. One needs to modify GR with the aim to avoid the fundamental issues of the theory. These issues are related to the existence of a singularity at the origin of most vacuum solutions of Einstein's equation, the inconsistency of GR with the quantum field, etc. On the other hand, the modifications to GR and the alternative theories of gravity may be considered a step forward in developing a unified theory of the interactions. One of the key principles of both modern physics and Einstein's GR is Lorentz's invariance. Another piece of evidence for the significance of the Lorentz invariance is considered to be the capability of GR to describe all observed gravitational events and its natural mathematical elegance [2]. Lorentz invariance, on the other hand, may not be a precise symmetry at all energies [3]. Any successful description must fail at some point, indicating the appearance of new physical degrees of freedom beyond that point. Lorentz invariance also causes divergences in quantum field theory, which can be rectified by a short cutoff distance that breaks it [4].
Attempts to solve problems in modern physics that are outside the scope of local Lorentz symmetry are fascinating [5]. Einstein AEther (EAE) theory is considered a covariant-based modified theory of gravitation, which violates the Lorentz symmetry locally. Its action consists of the Einstein-Hilbert term coupled with a dynamical, unit timelike vector field, v α , named as AEther. The rotational symmetry in a preferred frame is retained, while local boost invariance v α is broken [6]. Thus AEther is a type of locally favored state of rest at each location in spacetime as a result of unexplained physics. The Lorentz symmetry breaking in the gravity sector can be effectively described by the EAE theory, which has been extensively used to figure out the quantitative restrictions on Lorentz-violating gravity. In contrast to general relativity, the presence of AEther field in AEther theory establishes a preferred timelike direction that violates Lorentzian symmetry [7].
The motion of charged or neutral particles around black holes (BHs) is one of the most fascinating issues in BH astrophysics. It is essential for determining spacetime's geometric structure. The study of the general relativistic motion of particles and the electromagnetic fields near to BHs are now motivated by new observational and theoretical evidence for BHs. Astronomical investigations over the past ten years have shown the existence of supermassive BHs as well as stellar mass in galactic centers and X-ray binary systems. It has long been known that stellar-mass BH binaries exhibit quasi-periodic oscillations in the X-ray flux light curves, and this phenomenon is regarded as one of the most effective tests of strong gravity models. These fluctuations are closely proximity to the BH, and expressed frequencies scale inversely with the BH mass. According to the current developments, we can accurately measure the frequencies of QPOs, at center and its surroundings. It is possible to classify quasi-periodic oscillations (QPOs) into subclasses based on their measured frequencies less than 0.5 kHz. Most of these are high-frequencies (HF) and low-frequencies (LF) QPOs with frequencies up to 500 Hz and up to 30 Hz, respectively. The twin peaks with a frequency ratio 3:2 close to are often steady and detectable indicators of HF QPO oscillations in BH microquasars [8]. The QPOs around BHs (nonrotating:rotating) and wormholes have been discussed in several papers [9][10][11][12][13].
It is quite intriguing to investigate the properties of BHs in the presence of an AEther field. In the EAE theory, two spherically symmetric and static BH solutions have recently been presented, with the two combinations in coupling constants [14,15]. Additional spherically symmetric BH solutions in the background of a class of coupling constants have been examined via numerical computing [7], and their analytical description in polynomial form has been used in the exploration of quasi-normal modes in EAE theory [16,17]. The neutral particle dynamics surrounding non-spinning EAE BHs have been investigated in [18], while the shadow and deflection angles for slowly rotating EAE BHs have been probed in [19].
A productive method for exploring the phenomena surrounding BHs involves studying QPOs observed in microquasars. Precise frequency measurements of QPOs can provide valuable insight into the central object. Various types of QPOs have been categorized based on their frequency, which can range from just a few millihertz to 0.50 kHz. The concept of QPOs are characterized by low and high-frequency ranges. BH microquasars typically exhibit high-frequency QPOs with dual peaks, whose frequency ratio closely corresponds to the ratio [20]. Multiple studies have investigated high-frequency QPOs of neutral [21], spinning [22] and charged test particles [23,24] in the vicinity of rotating and non-rotating BHs. Recently, alternative rotating spacetimes have been explored using hot-spot data surrounding Sgr A* [25].
In the present paper, we study the orbital and epicyclic motion of neutral test particles in the background of the first and second kinds of slowly rotating EAE BHs. We obtain the analytical expressions for specific energy and specific angular momentum of equatorial circular orbits and investigate the properties of both co-rotating as well as contrarotating ISCOs. We explore the perturbed motion of stable circular orbits located in an equatorial plane and examine the radial profiles of the radial, vertical, and axial frequencies in dependence on the BH mass and parameters of the BHs. Furthermore, we examine the position of resonant radii of HF QPOs models, i.e., epicyclic resonance (ER) and its variants, relativistic precession (RP) and its variants, tidal disruption (TD), as well as warped disc (WD) models. The Periastron and Lense-Thirring Precession have also been discussed.
Greek indices are taken to run from 0 to 3. However, for expressions having astrophysical relevance we use the physical constants explicitly.

Slowly rotating Einstein-AEther black hole
The slowly rotating EAE BH solution is an asymptotically flat solution of the field equations in EAE theory. The line element describing the geometry of this BH with mass M, and electric charge Q can be given as [26,27] Two types of exact solutions exist for slowly rotating EAE BHs. The first solution (termed as the first kind of slowly rotating EAE BH) corresponds to the special choice of coupling constant, i.e., c 14 = 0, c 123 = 0, where c 123 = c 1 + c 2 + c 3 , c 14 = c 1 + c 4 , and the metric function f (r ) for this case can be written as where c 13 = c 1 +c 3 . For c 13 = Q = 0, and a = c 13 = 0, the first kind of slowly rotating EAE BH reduces to the classical slowly rotating Kerr BH, and RN BH, respectively, while a = c 13 = Q = 0, leads to the Schwarzschild BH. The outer horizon for the first kind of slowly rotating EAE BH can be found by solving f (r ) = 0. The second solution (termed as the second kind of slowly rotating EAE BH) corresponds to c 123 = 0, and the metric function f (r ) for this case can be written as For c 13 = c 14 = Q = 0, the second kind of slowly rotating EAE BH reduces to the classical slowly rotating Kerr BH. Moreover, c 13 = c 14 = a = 0, leads to the RN BH, while c 13 = c 14 = a = Q = 0, corresponds to the Schwarzschild BH. The outer horizon for the second kind of slowly rotating EAE BH is situated at There are the number of observational and theoretical bounds on coupling constants c i . In the present work, we impose the following constraints [28] 0 ≤ c 14 < 2, 0 ≤ c 13 < 1, 3c 2 + c 13 + 2 > 0.
The positions of horizons for the first and second kinds of slowly rotating EAE BHs have shown in Figs. 1 and 2, respectively. Note that the horizons of slowly rotating EAE BHs do not depend on the rotation parameter a of BH. The horizon of the first kind of BH decreases with the increase of charge Q, while it increases as the coupling parameter c 13 increases. However, for the second kind of BH, the radii of horizons decrease as the coupling parameter c 14 increases. The EAE BHs have a greater radius of the horizon as compared to RN BHs.

Circular orbits around EAE BH
The motion of a neutral particle can be described by the Hamiltonian given by where m is the mass of the particle, p α = mu α represents the four-momentum, u α = dx α /dτ denotes the four-velocity, and τ is the proper time of the test particle. The Hamilton equations of motion can be written as where ζ = τ/m is the affine parameter. Due to the symmetries of the BH geometry, there exist two constants of motion, namely specific energy E and specific angular momentum L, given by where E = E/m, L = L/m and the equations of motion in an equatorial plane can be written in the form [29,30] where R(r ) takes the form Using normalization condition g νσ u ν u σ = −1, one can write whereṙ = dr/dτ ,θ = dθ/dτ , and V e f f denotes the effective potential given by a relation The effective potential V e f f (r, θ) plays an important role to illustrate the motion of test particles. One can describe the motion of a particle with the help of V e f f (r, θ) without using the equations of motion. The circular orbits for equatorial plane θ = π/2 are given by simultaneous conditions In order to find the specific energy E and specific angular momentum L of circular orbits, we follow the formalism of forces presented in [31,32], and the corresponding expressions for the case of slowly rotating EAE BHs are given by the relations where the subscripts (I ) and (I I ) stand for the first kind and second kind of slowly rotating EAE BHs, respectively. The expressions for P (I ) (r ) and P (I I ) (r ) are given by ). (23) The locality of stable or unstable circular orbits is consistent with the minimum or maximum of the effective potential correspondingly. In Newtonian theory, the effective potential has a minimum for any value of the angular momentum, and then it has no minimum radius of a stable circular orbit (I SC O) [33]. But this position is altered when the effective potential has a difficult form liable to the particle angular momentum and other parameters. Therefore, in GR and for the particles moving near the Schwarzschild BH, the effective potential has two extrema for any value of angular momentum. But, only for a particular value of angular momentum do the two points happen together. This point presents I SC O where is placed at r = 3r g [33][34][35][36][37][38] where r g is the Schwarzschild radius.
The graphical behaviour of energy E and angular momentum L of equatorial circular orbits around the first and second kind of slowly rotating EAE BHs is depicted in Figs. 3 and 4 respectively. The energy and angular momentum of circular orbits around the first and second kinds of slowly rotating EAE BHs increases as the coupling parameter c 13 increases, while they decrease with the increase of coupling parameter c 14 , charge Q or rotation of BH. We compare the energy as well as the angular momentum of circular orbits around Schwarzschild, RN, Solwly rotating Kerr and slowly rotating EAE BHs, and observe that the circular orbits around slowly rotating EAE BHs have greater energy and angular momentum among all the cases.
The ISCO is located at d 2 V e f f /dr 2 = 0 and the position of co-rotating and counter-rotating ISCOs for first and second kinds of slowly rotating EAE BHs is illustrated in Figs. 5 and 6,    respectively. The co-rotating ISCOs shift towards the BHs, while the contra-rotating ISCOs move away from the BH with the increase of the coupling parameter c 14 or electric charge Q. In the case of co-rotating particles, the ISCOs shift away from the BH as the coupling parameter c 13 is increased but move towards the BH when the rotation of BH is increased. Smaller radii of ISCOs can be observed when the BHs rotate rapidly. We see that the ISCOs around the first and second kinds of slowly rotating EAE BHs are smaller as compared to slowly rotating Kerr BH.

Harmonic oscillations as perturbation of circular orbits
In order to investigate the oscillatory motion of neutral particles, we perturb the equations of motion in the vicinity of stable circular orbits. If a test particle is slightly shifted from the equilibrium position related to a stable circular orbit situated in an equatorial plane, it will undergo epicyclic motion, characterized by linear harmonic oscillations. The frequencies of harmonic oscillatory motion measured by the local observer are given by [40,41] The radial (ω r ), latitudinal (ω θ ), and orbital/axial (ω φ ) frequencies of the neutral test particle for the first and second kinds of slowly rotating EAE BH takes the form where the introduced coefficients read α = (σ Y +r 2 (Q 2 +3Q(r −2)r +2r 2 (r (3r −8)+6))), , respectively. The identification of different shapes of charged particle epicyclic orbits near a stable circular orbit can be facilitated by observing the behavior of the fundamental frequencies, namely ω r , ω θ , and ω φ , as well as their respective ratios. The Newtonian theory of gravitation predicts that all frequencies are identical, resulting in Fig. 7 Radial profiles of frequencies of small harmonic oscillations radial ν r , vertical ν θ and axial ν φ of a neutral particle around first kind of slowly rotating EAE BH with BH mass M = 10M , measured by a static distant observer elliptical trajectories for particles orbiting spherically symmetric bodies. However, for Schwarzschild BH, the frequencies follow the relation ω r < ω θ = ω φ , causing a periapsis shift and inducing relativistic precession as the radius of the orbit decreases and the strong gravity region is approached.

Frequencies measured by distant observer
The locally measured angular frequencies are related to the angular frequencies measured by the static distant observer ( ), by the gravitational redshift transformation where dτ/dt is the redshift coefficient. If the frequencies of small harmonic oscillations measured by the distant observer are expressed in physical units, one needs to extend the corresponding dimensionless form by the factor c 3 /G M. Thus the frequencies of the neutral particles measured by distant observers are given by [42] ν j = 1 2π where j ∈ {r, θ, φ}; r , θ , and φ denotes the dimensionless radial, latitudinal, and axial angular frequencies measured by a distant observer. The Figs. 7 and 8 display the radial profiles of frequencies ν j for small harmonic oscillations of neutral particles, as measured by an observer at a distance, around slowly-rotating EAE BHs of both the first and second kinds, with varying spin parameter a, charge Q, and dimensionless coupling parameters c 13 and c 14 . As any of these parameters increase, the radial profiles shift towards the black hole. However, coupling parameter c 13 contributes to moving the radial profiles away from the BH. The particles moving around slowly rotating Kerr BH have high frequencies as compared to moving around slowly rotating EAE BHs.

Quasi-periodic oscillations models
The previous section's exploration of neutral particle oscillations on circular orbits has potential ramifications in astrophysics, particularly about HF QPOs detected in microquasars. Hot spot models hypothesize that radiating spots exist in thin accretion discs and follow nearly circular geodesic trajectories. The ER model, on the other hand, posits resonance of axisymmetric oscillation modes of accretion discs. The frequencies of these oscillations are linked to the orbital and epicyclic frequencies of circular geodesic motion. Within the ER model, which includes axisymmetric oscillatory modes at frequencies ν θ and ν r , uniform radiation is assumed to emanate from the oscillating torus (or circle). The presence of a sufficiently large irregularity on the orbiting torus, which has a frequency of ν φ , permits the establishment of the nodal frequency associated with this irregularity. Table 1 contains the variations of the ER model, for details, see [39]. In the standard RP model [43], the higher of the paired frequencies corresponds to the orbital (azimuthal) frequency, ν U = ν φ , while the lower frequency represents the periastron precession rate, ν L = ν φ − ν r . From the vari- Table 1 The upper and lower frequencies of HF QPOs models [39] Model ants of the RP model, we select the RP1 model introduced in [44], where ν U = ν θ and ν L = ν φ − ν r , and the "total precession model" RP2 introduced in [9], where ν U = ν φ and ν L = ν θ − ν r (see Table 1). Both the RP1 and RP2 models predict frequencies ν U and ν L close to those of the RP model. The TD model, which has ν U = ν φ +ν r and ν L = ν φ , may bear some resemblance to hot spot models. This is because numerical simulations have shown that the BH tidal forces that disrupt inhomogeneities (such as asteroids) can create a ring-like structure that includes an orbiting radiating core. As for the WD oscillation model of twin HF QPOs, it assumes that there are non-axisymmetric oscillatory modes of a thin disc. In order to account for the frequencies in Table 1 of the WD model, we must make the assumption that there are also vertical axisymmetric oscillations of the thin disc, which we will call ν θ . The radial profile of the frequencies ν U and ν L of different HF QPOs models, i.e., ER (ER0, ER1, ER2, ER3, ER4, ER5), RP (RP1, RP2, RP3), TD and WD models for both first and second kind of slowly rotating EAE BHs are illustrated in Figs. 9 and 10, respectively. We compare the frequencies in the background of HF QPOs models modified for the second kind of slowly rotating EAE BH, while gray lines are for Kerr BH. BH spin parameter a = 0.3 has been used in all cases by the particle motion around a slowly rotating EAE BHs with the motion around a slowly rotating Kerr BH and see that the particle motion deviates from the Kerr limit.

Resonant radii
The quasi-harmonic character of the motion of test particles trapped in a toroidal space around the equatorial plane of BHs gives an interesting astrophysical application concerning the HF QPOs observed in the LMXB systems containing a neutron star or BH. The HF QPOs usually come in pairs of the upper and lower frequencies of twin peaks in the Fourier power spectra. The peaks of HFs are close to the axial frequency of the stable circular orbit denoting the inner edge of Keplerian discs orbiting BHs, thus the strong gravity effects must be relevant in explaining HF QPOs [45].
The HF QPOs usually appear in a pair of two peaks with upper ν U and lower ν L frequencies in the timing spectra, and the ratio of frequencies ν U : ν L is close to the ratio 3:2. Observation of this effect in different non-linear systems shows the existence of the resonances between two modes of oscillations. In the case of slowly rotating EAE BHs, the upper and lower frequencies of a neutral test particle around the first and second kinds of slowly rotating EAE BHs are functions of parameters c 13 , c 14 resonance position r , and a BH spin a, given by M, a, c 13 , c 14 ), ν L = ν L (r, M, a, c 13 , c 14 ). (36) These frequencies ν L and ν U are inversely proportional to the mass M of the BH and given by Eq. (35). In order to fit the frequencies observed in HF QPOs with the BH parameters, one needs to find the resonant radii, given by a relation In general, the resonant radii r 3:2 are given by a numerical solution of higher order polynomial in r , for given values of the BH parameters. The resonant radius solution has no dependence on the mass of the BH as Eq. (37) is independent of the mass of the BH explicitly. The position of resonant radii r 3:2 for HF QPO models in the background of slowly rotating first and second kinds of EAE BHs have been shown in Figs. 9, and 10, respectively. We compare the case of slowly rotating Kerr BH with the slowly rotating EAE BHs and see that for all considered BH QPOs models, the position of resonant radii r 3:2 for slowly rotating EAE BHs is at smaller radii as compared to slowly rotating Kerr BH. The resonant radius r 3:2 for the ER4 HF QPOs model is closer to the BH among all considered models.

Periastron and Lense-Thirring precession
In this section, we discuss the periastron and Lense-Thirring precession frequency of a neutral test particle moving around the slowly rotating EAE BHs in the limit of a small perturbation with respect to the equatorial plane π/2. In order to calculate the Periastron precession, we assume that the particle is slightly perturbed from its stable position which in turn results in oscillations about the stable position rs with a radial frequency r . The periastron P frequency is defined as a difference in orbital frequency φ and radial frequency r , while the Lense-Thirring precession frequency LT is Fig. 11 Periastron and Lense-Thirring precessions for the first kind of slowly rotating EAE BH defined as the difference in orbital frequency φ and latitudinal frequency θ , given by the relation where the orbital φ , radial r , and latitudinal frequencies θ are given by Eq. (34). In contrast to Newtonian gravity, when general relativistic effects are taken into consideration (close to BH), the radial and orbital frequencies are no longer equal, i.e., θ = φ . The graphical bahaviour of periastron and Lense-Thirring precession frequencies for both first and second kinds of EAE BHs have been depicted in Figs. 11 and 12, respectively. Both periastron and Lense-Thirring precession frequencies increase as the coupling parameters c 13 or c 14 increase. However, P decreases while P increases when the BH rotates rapidly. The radial profiles of P for the case of slowly rotating Kerr BH lie above the first and second kinds of slowly rotating EAE BHs, while LT lies below the EAE BHs.

Discussion and conclusions
We have explored the orbital and oscillatory motion of neutral test particles moving around slowly rotating first and second kinds of EAE BHs and examined the influence of the BH parameters. The positions of horizons have been investigated. The radii of the BH horizon increase with the increase of the coupling parameter c 13 , however, it decreases as the coupling parameter c 14 or electric charge of BH Q is increased.
Using the formalism of forces [31], we have found the analytical solutions of radial profiles of specific energy and specific angular momentum of equatorial stable circular orbits in dependence on the BH parameters and compared them with the case of slowly rotating Kerr BH. Initially, the specific energy and specific angular momentum of circular orbits grow fastly as the coupling parameter c 13 is increased, while it decreases with the increase of charge parameter Q of BH, and then becomes almost constant when the particles move away from the BH. The circular orbits close to the BH have less energy as well as angular momentum when the BH rotates rapidly or the coupling parameter c 14 is enhanced but increases as the radial distance is increased. The circular orbits with high energy and large angular momentum can be seen around slowly rotating Kerr BH as compared to the first and second kinds of slowly rotating EAE BHs.
We have determined the positions of ISCOs for the corotating and contra-rotating neutral test particles orbiting slowly rotating first and second kinds of EAE BHs. We observed that the co-rotating ISCOs shift towards the BHs, while the contra-rotating ISCOs move away from the BH with the increase of the coupling parameter c 14 or electric Fig. 12 Periastron and Lense-Thirring precessions for the second kind of slowly rotating EAE BH charge Q. However, the co-rotating ISCOs shift away from the BH as the coupling parameter c 13 is increased but move towards the BH when the rotation of BH is increased. It is noted that when the BHs rotate rapidly, smaller radii of ISCOs can be observed. The ISCOs with smaller radii can be found around the first and second kinds of slowly rotating EAE BHs as compared to slowly rotating Kerr BH.
We have studied the fundamental frequencies of small harmonic oscillations of test particles in the radial and vertical directions related to the equatorial circular orbits and studied the radial profiles of the radial, vertical, and orbital frequencies. The orbital and vertical frequencies coincide when the rotation of BH is switched off, however, different profiles can be observed when the BH rotation is switched on. We found that the radial profiles move toward the BH when the BH rotates rapidly or charge Q or parameter c 14 is increased. However, coupling parameter c 13 contributes to moving the radial profiles away from the BH. The particles moving around slowly rotating Kerr BH have high frequencies as compared to moving around the first and second kinds of slowly rotating EAE BHs.
We have explored the radial profiles of upper ν U (r ) and lower ν L (r ) frequencies for many different HF QPO models, i.e., ER model and it's variants (ER0, ER1, ER2, ER3, ER4, ER5), RP model and it's variant (RP1, RP2), TD and WD models and determined the position of resonance radii r 3:2 . It is worth mentioning that the positions of resonant radii r 3:2 for all considered HF QPO models for the case of slowly rotating first and second kinds of EAE BHs lie at smaller radii as compared to the slowly rotating Kerr BH. Furthermore, the resonant radius r 3:2 for the ER4 HF QPOs model is closer to the BH among all considered HF QPO models. In addition, we have also analysed the bahaviour of periastron and Lense-Thirring precession frequencies for both the first and second kinds of slowly rotating EAE BHs. Both periastron and Lense-Thirring precession frequencies can be seen high with the increase of the coupling parameters c 13 . However, these frequencies are small when BH rotation or coupling parameter c 14 is increased.