Circular orbits and accretion process in a class of Horndeski/Galileon black holes

In this paper the geodesics motion and accretion process around a subclass of Horndeski/Galileon black holes are investigated. Firstly, we present spherically symmetric geometries in a Horndeski/Galileon black hole spacetime by considering an isothermal fluid around the black hole. Then we focus on three main issues: in the first step circular orbits of test particles and their stability in equatorial plane are examined in details. Then, by treating perturbations via restoring forces, oscillations of particles around the central object are studied. Finally, the accretion process, the critical speed of the flow and accretion rate are investigated in this setup properly.

general class of scalar-tensor field models with second-order field equations and it may be considered as a proper generalization of General Relativity in high energy regime (Latosh 2016). We consider the following action which is a subclass of general Horndeski/Galileons class Here R is the Einstein-Hilbert term, G µν is the Einstein tensor, φ is the scalar field, Λ is a cosmological constant term, ζ > 0, η and β are model parameters. On astrophysical scales, which covers also static and spherically symmetric solutions, this scalar-tensor theory may play a crucial role. On the other hand, any modification of General Relativity must be consistent with astrophysical observations. So, it is important to see how astrophysical processes such as accretion onto black holes work in this scalar-tensor framework and can be used also as a probe to see viability of these theories from experimental viewpoint. The issue of black hole accretion disk is studied in some subclasses of scalar-tensor theories. However, there is a gap in literature since black hole accretion disk has not been studied in the mentioned scalar-tensor theory as a subclass of the general Horndeski/Galileon scenario. This is the motivation of the present study and we are going to fill this gap in this paper. Accretion disks are constructed by rotating gaseous materials that move in bounded orbits because of the gravitational force of central mass, such as Young Stellar Objects (YSO), main-sequence stars (MSs), neutron stars (NSs), and supermassive black holes in Active Galactic Nuclei (AGN). In such systems, particles orbits are stable, but when the orbits of these materials become unstable, following it, accretion will be happened. Accretion is the process by which a massive central object such as a black hole captures particles from a fluid in its vicinity. The particles which accelerate from rest must be passed through a critical point, the point where the velocity of the gas matches its local sound speed. Then the gas falls onto the central mass at supersonic velocities. This process leads to increase in mass of the black hole (Martnez 2014). In addition, extra energy would be released in this process where this energy can be source of some astrophysical phenomena, such as the production of powerful jets, high-energy radiation, and quasars (Kato et al. 2008). Therefore, the study of the geodesic structure of particles in the vicinity of black holes and specially investigation about some characteristic radii such as marginally bound orbits (r mb ) and innermost stable circular orbits (r isco ) are interesting issues for a careful study of the subject matter. These radii are very important in the study of black hole accretion disks. For example, in thin accretion disks, the inner edge of the disk coincides with the innermost stable circular orbit (ISCO) and the efficiency of the energy released, which describes the significance of converting rest-mass energy into radiative energy (Xie et al. 2012), can be determined from this radii.
The location of unstable or stable circular orbits is consistent with the maximum or minimum of the effective potential respectively. In Newtonian theory, for any value of the angular momentum, the effective potential has a minimum and then stable circular orbit is free to have arbitrary radius, that is, there is no minimum radius of stable circular orbit, (ISCO) (Kaplane 1949). But this situation is different when the effective potential has a complicated form depending on the particle angular momentum and other parameters or when one incorporates general relativistic effects. For example, in General Relativity and for particles moving around the Schwarzschild black hole, for any value of the angular momentum, the effective potential has two extrema (minimum or maximum). But, only for a specific value of the angular momentum the two points coincide. This point introduces ISCO where is located at r = 3r g (Landau 1993, Kaplane 1949 where r g is the Schwarzschild radius. In different metrics, the properties of spacetime affects the locations of these radii and some parameters such as: specific energy, angular momentum and angular velocity are important in the position of these points. A lot of research programs are devoted to study these radii and their physical significance. Ruffini et al. (1971) and Bardeen et al. (1972) studied the properties of innermost stable circular orbits around the Kerr black hole. Even Hobson et al. (2006) described these features in details in their textbook on General Relativity. The radiation efficiency of accretion disks, η, for Schwarzschild and Kerr black holes was obtained by Novikov and Thorne (1973) which its value lies in the range of 0.057 − 0.43 depending on the black hole spin. The Kerr-like metric was constructed by Johannsen and Psaltis (2011) and then Johannsen (2013) has studied the accretion disks around such black holes. The study of the geodesic motion and the circular orbits of charged particles around weakly magnetized rotating black holes are carried out by Tursunov et al. (2016).
In an accretion disk particles move in stable orbits but when a perturbations, as a result of restoring forces, act on the particles, oscillations around the circular orbit can take place in vertical and radial directions with epicyclic frequencies. Happening the oscillations (in response to perturbations) in the inner region of an accretion disk is another important characteristic of these regions. Oscillations can be source of strong and chaotic time variations in spectrum of such systems. Therefore, study about orbital and epicyclic frequencies (radial and vertical) play an important role in the physics of relativistic accretion disks around the black holes. Isper (1994Isper ( , 1996, Wagoner (1999), Kato (2001) and Ortega-Rodriges et al. (2006) have studied in this field. Resonance between such a frequency modes which proposed by Kluzniak and Abromowicz (2000) can be a physical mechanism for existing Quasi-periodic oscillations (QPOs). QPOs in the X-ray fluxes of some astrophysical objects such as a neutron star and black hole sources have been reviewed by many researches including van der Klis (2000) and McClintock et al. (2003). Johannsen (2013) has examined the radial and vertical epicyclic frequencies in the Kerr-like metric.
With these preliminaries, in this paper we study nonrotating black hole solutions with accretion disk in a subclass of Horndeski/Galileons spacetime general class. For simplicity we restrict our study to equatorial plane in a polar coordinates system. Firstly, the singularity and event horizon in this spacetime geometry are presented. Then, in order to investigate the circular orbits, effective potential is obtained in this setup. We study the locations of several characteristic radii, such as: marginally stable circular orbits r isco , marginally bounded circular orbits r mb and photon orbits r ph in equatorial plane. Also, the ISCO binding energy, the maximum radiation efficiency, the emission and temperature in equatorial epicyclic frequencies are computed. Finally, some dynamical parameters and critical accretion of isothermal fluid are investigated in details.
This paper is organized as follows: In section 2 we introduce Horndeski/Galileon spacetime. The general formalism of a test particle's motion is discussed in section 3 where circular motion, stable circular orbits and oscillations are examined in subsections 3.1, 3.2 and 3.3. In section 4 and it's subsections, the general form of some dynamical parameters such as critical speed of the flow, accretion rate and the time of accretion for an isothermal fluid are obtained. In section 5 we have explained physically all of these results for a subclass of general solution of the Horndeski/Galileon black hole. Finally, section 6 is devoted to summary and discussion.

Horndeski/Galileon Spacetime
We study static and spherically symmetric limit of black hole solutions in a subclass of general Horndeski/Galileon theories (see Maselli (2015) for the case of slowly rotating black holes in Horndeski theory and also Babichev et. al. (2016) for black hole and star solutions for Horndeski theory). The general form of the line-element for such systems with the metric signature (+, −, −, −) is described by The metric functions, f (r) and h(r), depend only on the radial coordinate r and are given as follows (Tretyakova 2016) In these relations, µ plays the role of the mass term and k is obtained from the following constraint equation where C 0 is a constant of integration. It is important to note that the static metric (2) has the time rescaling symmetry. So, if h(r) is a solution then ch(r) should be also a solution where c is a constant. This means that c should be set for t to describe the proper time in the relevant region, namely h → 1. We note also that equation (5) is not a solution, but rather an ansatz on the scalar field. This type of scalar configuration was firstly considered by Babichev and Charmousis (2014) in which the solution for ψ(r) is also presented. About the stability of solutions, Ogawa et al. (2016) and Takahashi and Suyama (2017) showed that solutions with nonzero q are generically plagued by ghost or gradient instability. However, for solutions with q = 0, there exist some stable solutions (see Kobayashi et al. (2012), Takahashi and Suyama (2017), and Tretyakova and Takahashi (2017)). We note that recently it has been pointed out by (Babichev et. al. (2017) and (2018)) that the conclusion of Ogawa et al. (2016) and Takahashi and Suyama (2017) is incorrect. Now from the recent relations, various solutions for different values of C 0 and q can be obtained. We are going to discuss a common expression in this paper.
Now we study the properties of a subclass of black hole solutions in Horndeski/Galileon gravity as has been introduced above. The mentioned Horndeski spacetime contains a singularity at the location where the following condition holds f (r) = 0 .
The event horizon is a null surface. A surface that is defined as f (x µ ) = 0 will be null if where n µ is the normal 4-vector to the surface and it is defined as n µ = ∇ µ f . Since we are interested in to study the problem in the equatorial plane, then the relation (8) can be written as g rr (∂ r f ) 2 = 0. Therefore, in the radial distance that g rr = 0 or equivalently f (r) = 0, we would have an event horizon. The location of the event horizon in the geometry is a radial distance from the center of the core where the metric is singular, except the intrinsic singularity which cannot be removed via coordinate transformation.

Test Particle's Motion: General Formalism
The motion of a test particle is governed by the geodesic structure of the underlying spacetime manifold. In this section, we study general form of timelike geodesics around a subclass of Horndeski/Galileon black hole. Spacetime around this black hole is static and symmetric with two Killing vectors ξ t = ∂ t and ξ ϕ = ∂ ϕ which imply two constants of motion E and L (conserved energy and angular momentum per unit mass) along the trajectory as follows where u µ = (u t , u r , u θ , u ϕ ) is the four-velocity of the test particle. Using the normalization condition for four-velocity, that is u µ u µ = 1, we have From equations (9) and (10) and in equatorial plane with θ = π 2 , four-velocity will be given by the following components Also the following equation can be derived easily In this equation V ef f is the effective potential for the test particle motion that is given by It is clear that effective potential depends on the particle's specific angular momentum radial distribution and the spacetime parameter via h(r). The study of effective potential is very useful in geodesic motion. For example, the local exterma of the effective potential determine the location of the circular orbits.

Circular Motion
For circular motion in the equatorial plane, radial component r must be constant and so u r =u r = 0 must be satisfied. Therefore, from equation (12) we would have V ef f = E 2 and d dr V ef f = 0. From these relations the specific energy E, the specific angular momentum L, the angular velocity Ω ϕ and angular momentum l can be obtained by the following relations respectively In order the energy and angular momentum to be real, the following condition must be satisfied By solving this inequality, the limited area of circular orbit can be obtained. Therefore, this is the condition for existence of the circular orbits. For bound orbit the relation E 2 < 1 must be hold and in marginally bound orbits we have E 2 = 1. Then from equation (14) we fined By solving this equation, marginally bound orbits can be obtained easily. From equations (14) and (15) it is seen that the energy and angular momentum diverge at the radius where the following relation holds Photon sphere can be obtained by solving this relation. In a photon sphere, photon moves on circular orbits. This region plays a crucial role in the study of gravitational lensing, since lensing effect cannot be observed below this region.

Stable Circular Orbits and Radiant Energy Flux
The local minima of the effective potential correspond to the stable circular orbits. Thus a stable circular orbit exists if d 2 dr 2 V ef f > 0 and in addition to this condition, in marginally stable circular orbits, r isco , the condition d 2 dr 2 V ef f = 0 must be satisfied. From equation (13) we have Accretion process is possible in r < r isco . When falling particles from rest at infinity accrete onto the central mass, the released gravitational energy of falling particles can convert into radiation where this energy is the source of the most energetic phenomena in astrophysics. The flux of the radiant energy over the disk can be expressed in terms of the specific angular momentum L, the specific energy E and the angular velocity Ω ϕ by the following relation (see for instance Kato et al., (2008)) whereṀ is the accretion rate, Ω ϕ,r ≡ dΩϕ dr and the parameter g is determinant of g µν given by We set sin θ = 1, since we restrict our studies in equatorial plan. From relations (14)- (16) we would have where by definition The steady-state accretion disk model is supposed to be in thermodynamical equilibrium. Then the radiation emitted from the surface of the disk can be as a black body radiation. So, the relation K(r) = σT 4 (r) can be hold between energy flux emitted at the surface of the disk and effective temperature of the disk (σ is the Stefan-Boltzman constant). Using this relation, temperature distribution on the disk by assuming thermal black body radiation can be obtained easily and then we can compute the luminosity L(ν) of the disk. The observed luminosity at the distance d to the source with the disk inclination angle γ has the following form (Torres 2002) where I(ν) is the thermal energy flux. In this relation r i indicates the position of the inner edge and we take r i = r ms . Also r f indicates the outer edge of the disk. Since for any kind of general relativistic compact object the flux over the disk surface could be vanishing at r → ∞, we take r f → ∞. The emitted frequency is given by ν e = ν(1 + z) where the redshift factor z, by neglecting the light bending, can be written as follows The efficiency of the accreting flow is another important characteristic of the mass accretion process. The maximum efficiency of transforming gravitational energy into radiative flux of such particles between innermost circular orbit and infinity, η * , is defined as the ratio of the specific binding energy of the innermost circular orbit to the specific rest mass energy which is given by the following relation where E isco is the specific energy of a particle rotating in an innermost stable circular orbit. This relation is valid for the case where all the emitted photons can escape to infinity. Now we focus on perturbations. If a perturbation acts on the fluid element, the motion of a test particle will be nearly circular orbit in the equatorial plane and the particle will oscillate around the circular orbit with three components of motion, the issue which is discussed in the next section.

Oscillations
In an accretion disk, various types of oscillatory motions as a result of restoring forces are expected. Restoring forces act on perturbations in the accretion disks resulting Horizontal and Vertical oscillations. Some of these restoring forces in accretion disks are resulting from rotation of the disk and from a vertical gravitational filed. When a fluid element is displaced in the radial direction, it will return to its equilibrium position due to a restoring force resulting from rotation of the fluid. In accretion disks, because of existence of central object, centrifugal force is balanced by the gravitational force. When the former dominates over the latter or the reverse happens, the element of flow will be pushed inward or outward to return to the original radius with epicyclic frequency Ω r . On the other hand, when a fluid element is perturbed in the vertical direction, the vertical component of the gravitational field returns the perturbed element toward equilibrium position, that is, the equatorial plane. As a result of this restoring force, the element of the fluid makes harmonic oscillation around the equatorial plane with vertical epicyclic oscillations Ω θ (Kato at al. 2008).
In a general relativistic discussion about the motion of the fluid in an accretion disk, three frequencies around the central object are important. Circular motion at the orbital frequency Ω ϕ , harmonic radial motion at the radial frequency Ω r and the harmonic vertical motion at the vertical frequency Ω θ . As we have stated, resonance between such frequencies can be source of quasi-periodic oscillations which leads to chaotic and quasi-periodic variability in X-ray fluxes from many galactic black holes. Study in this field is important in some sense. For this purpose, radial and vertical motions around a circular equatorial plane are discussed in this section. (10), to describe radial motion u θ = 0, and also for describing the vertical motion we have u r = 0. By setting u r = dr dτ = dr dt u t and u θ = dθ dτ = dθ dt u t , we find

Radial and vertical motions can be explained by
In order to investigate the radial and vertical epicyclic frequencies, small perturbations δr and δθ around the circular orbit in equatorial plane are considered. By taking the timederivative of the first equation in (29), equation describing the radial oscillations can be obtained as follows For a particle which is perturbed from its original radius at r = r 0 by a deviation δr = r −r 0 , the perturbed equation of motion is given by where a dote denotes differential with respect to time coordinate t and Ω 2 r ≡ − dr 2 . By the same procedure, for a perturbation in the vertical direction by a deviation given as where ef f . Then from equations (29) in equatorial plane we would have respectively and In these equations, a prime denotes differential with respect to the radial coordinate, r.
To proceed further, now we present basic dynamical equations in this subclass of general Horndeski/Galileons black hole spacetime.

Basic dynamical equations
In this section we provide the basic dynamical equations for our forthcoming arguments (we refer to Babichev et al., (2005) and (2013) for more details). Here we consider a perfect fluid which is specified by the following energy-momentum tensor where p and ρ are pressure and energy density of the fluid respectively. In this relation, u µ is the fluid elements four-velocity. Because of background symmetries, in relation (35) all of the components are functions of only the radial coordinate, r. Since we are assuming the fluid is flowing radially in the equatorial plane (θ = π 2 ), the general form of the four-velocity will be as follows where τ is the proper time along the geodesic. From this relation and under the normalization condition u µ u µ = 1, we obtain where for forward flow in time, u t must be positive and for accretion (inward flow), u r < 0. By deriving the energy-momentum and also particle-number conservation equations, all of the required equations for studying the accretion are obtained. Conservation of the energymomentum tensor is given by where in this relation (; ) shows the covariant differentiation, √ −g = r 2 sin θ h(r) f (r) and Γ is the second kind Christoffel symbol (affine connection) where its non-zero components are as follows From these relations, equation (38) yields where after some manipulations we obtain with A 0 as an integration constant. Projecting the energy-momentum conservation law onto the four-velocity via u µ T µν ;ν = 0, yields By considering the normalization conditions as u µ u µ = 1 and since g µν ;ν = 0, this relation reduces to Since A b ;a = ∂ a A b + Γ b ac A c , we find u r ρ ,r + (p + ρ) Γ 0 0c u c + (u r ,r + Γ 1 1c u c ) + Γ 2 2c u c + Γ 3 3c u c = 0 .
By using the non-zero components of the connection, this relation after some simplification yields which after integration, we would have where A 1 is an integration constant. Since in the left hand side u r < 0, the right hand side takes a minus sign too. So we find finally where A 2 is an integration constant. The equation of mass flux in this setup which is given by can be rewritten as Since we are interested in to study just in equatorial plane, the second term in equation (49) vanishes. Therefore, √ −gρu µ would be as a constant, that is where A 3 is an integration constant. We note that while we have restricted our study to equatorial plane by symmetry considerations, the general case is not so complicated in essence. Because of symmetry all important characteristics of the model can be obtained in θ = π 2 plane as well. Now we are going to determine dynamical parameter, critical accretion and accretion rate in this setup.

Dynamical parameters
To proceed further, we assume isothermal fluids. These fluids flow at a constant temperature. Therefore, p ∝ ρ and then the speed of sound throughout the accretion process remains constant. For such fluids, the equation of state is of the form p = kρ where k is the equation of state parameter. Then equations (46), (47) and (50) yield where A 4 is an integration constant. By substituting p = kρ, we can obtain u as follows Then from eq. (50) the density of the fluid can be obtained as Finally, from the relation p = kρ, the pressure can be obtained easily.

Mass evolution
In realistic astrophysical cases, the mass of the black hole is not fixed in essence. By some processes such as accreting of mass from accretion disk onto black hole and also Hawking radiation, its mass will be changed gradually. The rate of change of mass can be obtained by integrating the flux of the fluid over the surface of the black hole, that isṀ ≡ dM dt = − T r t ds where a dot denotes the time derivative, ds = √ −g dθ dϕ and T r t = (p + ρ)u t u r . By substituting these relations,Ṁ can be obtained as followṡ where A 0 = −A 1 A 2 , and A 2 = (p ∞ + ρ ∞ ) h(r ∞ ). Therefore, we obtaiṅ Now, time evolution of mass of the black hole with initial mass M i can be obtained by integration of equation (55) that can be rewritten as where F ≡ 4πA 1 (p + ρ) h(r ∞ ). By integrating from equation (56) we obtain is the critical accretion time. In the case t = t cr the denominator of equation (57) vanishes and the black hole mass grows up to infinity in a finite time. After determining the time evolution of the disk and black hole mass, now we study critical accretion in this setup.

Critical Accretion
Very far from the black hole, the flow is at rest but gravitational field of black hole tends to accelerate it inwards. When flow moves inwards, it must pass through critical point (sonic point) where in this point r = r c , the four-velocity of the fluid matches the local speed of sound, u = c s . In order to obtain sonic point, an expression for the radial velocity gradient with no other derivatives is required. From equations (50) and (51), the following two equations are obtained and From these equations we obtain where by definition and With these two equations, the following relation can be obtained The condition for critical points is D 1 = D 2 = 0. This condition yields and where index c refers to the critical point. Since the right hand side of the relation (64) must be positive, then if dependence of h(r) is known, by solving the following inequality the range of critical radius can be obtained Finally the sound speed c 2 s = dp dρ can be obtained from (52) as

A Subclass of Horndeski/Galileon Black Hole Solutions
As we have stated in section 2, a huge variety of solutions can be obtained from equations (3)-(6). Here we present a solution of this subclass of the Horndeski/Galileon setup that are characterized by the parameters A, B, C and γ (see for instance Tretyakova 2016). The solutions are as follows h(r) = C − µ r + Ar 2 + ∆ , and where we have set ∆ ≡ B tan −1 (rγ)

rγ
. The coefficients C, A, B and γ are defined as follows Here ǫ is a small parameter which marks the deviation from the inherent de Sitter solution. Also the scalar field in this case is given by where ψ = dφ dr . As we have said, singular points can be obtained by solving the relation f (r) = 0. From this relation singular point is located at r sing = µ B+C .
An important issue should be stressed here: about the behavior of the scalar field at the black hole horizon (say, r h ), the scalar field seems to be divergent at the horizon from Eq. (71). But since f (r h ) = 0 and df dr (r h ) = 0, one can expand f (r) in Taylor series to find in near horizon which gives for the scalar field φ(r) = φ 0 + φ 1 (r − r h ) 1/2 + φ 2 (r − r h ) 3/2 + ... . Therefore, the scalar field remains finite in the near horizon region (a similar analysis can be found in the paper by Miao and Xu (2016)). We note that although the scalar field itself does not diverge at the horizon, its derivative does. However, there is no physical divergence since all invariants, such as g µν ∂ µ φ∂ ν φ remain finite. For more discussion on this issue we refer to Feng et al. (2015) (see also Hadar and Reall (2017) and Caceres et al. (2017)).
We note that for nonzero η the model given by Eqs. (68)-(70) admits solutions in which the Λ-term in action (1) is totally screened. The metric then is not asymptotically flat but rather it is de Sitter with the effective cosmological constant proportional to η/β since the scalar kinetic term tends to a constant around the present time (Gubitosi and Linder, 2011). It is important to note that as has been shown by Tretyakova (2016), for 1−(B+C) < 3×10 −4 this metric must be equivalent to the Schwarzschild metric in the sense that it matches with the observations of the gravitational light deflection and perihelion precession. This feature guarantees the existence of bounds orbits for (B + C) > 0 (for more discussion on this issue see Tretyakova (2016)).
It is necessary to mention that for calculations, the assumption A ≃ 0 is considered. Then the components of metric (2) will be as h(r) = C − µ r + ∆ and f (r) = h(r) B+C . On the other hand, since γ ≃ √ A, due to smallness of γ, it is more suitable to substitute tan −1 (rγ) ≃ rγ and therefore ∆ ≃ B. All calculations are done with these approximations and also with assumption C = 1 and 0 < B + C < 1.125 (Tretyakova 2016). By these assumptions, we focus mainly on the role of parameter B as Horndeski/Galileon correction factor in our forthcoming treatment.

Circular Equatorial Geodesics
In order to investigate circular geodesics in equatorial plane, we need the explicit form of the effective potential which is governed by equation (12) as where ∆ ≡ B tan −1 (rγ) rγ . From the condition d 2 dr 2 V ef f > 0 for existence of the stable circular orbits, we see that for r < 3 equation (20), the location of the stable circular orbits would be at r ≥ 3 µ B+C . Then is introduced as the radius of the innermost stable circular orbit. The left panel of figure 1 represents the effective potential versus r for several values of the angular momentum L in the case with B = −0.1. We see that for L < 2.2 √ 3 no extremum can be observed and the first extremum is observed at L = 2.2 √ 3 (solid circle in the figure). This point represents the location of the innermost stable circular orbit located at r = 6.6. For larger values of the angular momentum, V ef f has two extremum where the maximum one denotes the location of the unstable circular orbit and the minimum one denotes the stable circular orbit. By increasing the angular momentum, V ef f will be larger and the maximum point turns to the smaller radii whereas the minimum point goes to larger radii. From the right panel of figure  1 the effect of Horndeski/Galileon correction factor B on the effective potential can be seen. The effective potential achieves larger values for larger values of B and by increasing this parameter, the loci of unstable circular orbit becomes closer to the central mass and stable orbits will be located farther from the central mass. Also the enhancement of the distance between these points by increasing B is obvious.
In addition to innermost stable circular orbit which is very important in studying the accretion around the black hole, there are other special radii where considering them is necessary. As we have stated previously, the circular orbits exist only for radii larger than the photon radius r ph . For r ph < r < r ms , the motion of the particle will be unstable against the small perturbations. This means that particle falls into the black hole or flee away to infinity. In the region r > r ms the particle moves on stable circular orbits.
From equations (18)-(20), other characteristic radii including the photon sphere r ph , circular orbit r circ and marginally bound orbit r mb can be obtained respectively as and We have plotted the characteristic radii versus the Horndeski/Galileon correction factor B in figure 2. The value of this parameter affects the location of the characteristic radii in the vicinity of the black hole. For larger values of B, the location of r ph , r isco and r sing will be closer to the black hole, whereas the behavior of r mb for negative and positive values of B = 0 is different. For B < 0, marginally bound radius decreases by increasing B, but for B > 0 this behavior becomes reverse as it finally matches the innermost stable circular orbit. Since r isco represents the inner edge of the accretion disk, we see that in larger values of B, the disk will be extended close to the central mass.
For a particle which moves in a circular orbit, the specific energy, specific angular momentum, angular velocity and angular momentum in equatorial plane can be derived as respectively. For the innermost stable circular orbit, these relations reduce to    The maximum efficiency at the innermost stable circular orbit is represented by a circle in each case in the right panel.    respectively.
In figure 3, the behavior of specific energy and angular momentum versus the radius are shown and the effect of the Horndeski/Galileon correction factor B is studied. The upper left panel represents variation of the specific energy versus r. In the upper right panel the location of this radius in innermost stable circular orbits is shown by the solid circle. The empty circle denotes the limit of the bound orbit. Increasing the Horndeski/Galileon correction factor enhances the energy and decreases the range of bound orbit radius. The lower left panel represents the angular momentum where it gets smaller values by increasing B. The loci of these parameters in innermost stable circular orbits are represented by circles on the lower right panel of the figure.
Now by knowing E in ISCO, we are able to determine the radiation energy efficiency of accretion. As we have said previously, the efficiency of accretion is 1 − E, where the maximum efficiency of accretion given by 1 − E isco in innermost stable orbit in this setup is 1 − 8 9 (B + C). In figure 4, the efficiency of accretion is plotted versus r for different values of the Horndeski/Galileon correction factor B where we have denoted the maximum with a circle in this figure. In the case B = 0 where coincides with the Schwarzschild black hole (dashed line), the efficiency equals to 0.057 at r = 6r g as usual. For negative values of the Horndeski/Galileon parameter, efficiency goes up and becomes greater than the Schwarzschild black holes ones. For example, in the case with B = −0.1, the efficiency equals to 0.105 at r = 6.66r g where is tangible. On the other hand, positive values of B is accompanied with decreasing of the efficiency. Now we can study radiation flux from the surface of the accretion disk by knowing E, L and Ω ϕ in equatorial plane and r ms . The Flux of the radiation energy of the accretion disk is obtained from equations (24) and (25) as where by definition Then the temperature can be obtained by using the equation K = σT 4 . In figure 5, the relation between the radiation flux and temperature is shown. The radiation flux has a maximum in the vicinity of the black hole and decreases at the smaller radii. By increasing the Horndeski/Galileon correction factor B, the energy flux raises and the maximum of emission flux tends to the smaller radii but the reverse happens after this point. The dependence on this parameter is very considerable in the vicinity of the black hole, but it is relatively weak far from the black hole. These behaviors are the same for the temperature. As we have said in subsection 3.2, the luminosity can be compute from equation (26), but because of complexity of the required equations, solving this equation is not possible analytically.

Epicyclic Frequencies
If a perturbation acts on a particle moving on a circular orbit in the equatorial plane, the particle experiences small oscillations in the vertical and radial directions. Using the epicyclic frequencies given by equations (33) and (34), we derive the radial and vertical epicyclic frequencies as follows and Ω 2 r = − 3 2 (1 + η β r 2 ) − 8 9 Aγ 4 r 8 ( 15 8 ∆ + C) + 5 3 Aγ 4 r 7 µ + γ 2 r 6 1 9 ∆γ 2 (3∆ + C) + A(− 16 9 C − 10 3 ∆ + B) + 10 3 µγ 2 r 5 (− 1 30 C − 1 5 ∆)γ 2 + A 1 3 µ 2 γ 2 + 1 9 (2∆ + B)C − 1 3 ∆(−2∆ + B) γ 2 + 7 9 A(B − 15 7 ∆ − 8 7 C) r 4 respectively. From equations (78), (86) and (87) it is clear that In the left hand side of figure 6 we have plotted epicyclic frequencies versus the radius for several values of the Horndeski/Galileon correction factor B in order to have a comparison. The angular and vertical epicyclic frequencies are shown in the left panel of this figure by the dashed lines which are coincide and will decrease by increasing r. We see explicitly that they are not dependent on parameter B. The radial epicyclic frequency is shown by the solid and dotted curves (Schwarzschild black hole) with a maximum where increasing the parameter B shifts this maximum to the smaller radii. The radial oscillation frequency increases by increasing the values of the parameter B. Dependence on this parameter is significant close to the black hole, but far from the central mass this effect is weak. We see that Ω r < Ω θ . The ratio of Ω θ Ωr is plotted in the right hand side panel of figure 6 which is a decreasing function of r. It is clear that in the vicinity of the black hole, this ratio is very greater than unity but far from the black hole it turns to unity and decreases by increasing B. In the left panel of figure 7, the locations of three particular resonances such as the parametric resonance with condition Ω θ Ωr = 3 2 and the forced resonance with ratio 3 : 1 and 2 : 1 are shown. The dependence of these characteristic radii to the metric parameter B is shown in the right panel of figure 7 where such radii are monotonically decreasing functions of this parameter. On the other hand, resonance will be happened in smaller distance from the central mass for larger values of the metric parameter, B.

Mass Evolution and Critical Points
In this Horndeski/Galileon accretion disk, the energy density and radial velocity for an isothermal fluid are given by and respectively. For critical points we find and It is clear that r c depends on the equation of state parameter k. This is means that the location of a critical point is not the same for all fluids. The profiles of radial velocity, density and accretion rate are presented in figures 8, 9 and 10 versus the dimensionless parameter r for equation of state parameter k = 1 2 and constants integration A 4 = 1.4 and A 3 = 1 respectively. The dotted, dashed and solid curves are corresponding to cases with B = 0, B = −0.05 and B = −0.1 respectively.
The left panel of figure 8 represents the radial velocity versus the radius. The locations of the critical points are marked by solid circles and their coordinates are stated in the right panel. The fluid has zero radial velocity far from the black hole and flows at the sub-sonic speed before the critical points. In a critical point, the speed of flow matches the speed of sound. After passing this point, in the vicinity of the black hole the speed of flow increases and turns to the super-sonic domain because of strong gravity. It can be observed that velocity decreases by increasing the Horndeski/Galileon correction factor B and the loci of the critical points get shifted to the black hole. Therefore, the speed of infalling particle reaches the speed of sound closer to the central mass.
The density profile of the fluid around the black hole for different values of B is shown in figure 9. Increasing the value of the parameter B increases the density. In addition, for such isothermal fluids, the mass of the black hole changes with time by the following relatioṅ We see that accretion rate for a general spherically symmetric static black hole in Horndeski/Galileon gravity is different from the case of a Schwarzschild black hole. In Horndeski/Galileon case, the accretion rate completely depends on the nature of the accreting fluid and also the metric parameter (here, parameter B). SoṀ > 0 for a normal fluid which satisfies (p + ρ) > 0. The change of accretion rate for different values of the parameter B with the same equation of state parameter are shown in figure 10. It is seen that accretion rate is higher in the vicinity of the black hole because of strong gravitational effect. Also, increment of the parameter B enhances the accretion rate. By using equation (57), the critical accretion time and the mass of the black hole are given by respectively. We see that for normal fluid, the black hole mass increases by accretion matter subject to Horndeski/Galileon gravity and increasing correction factor B will increase the black hole mass further. This behavior is studied by Rodrigues et al. (2009) for a Schwarzschild black hole in the presence of a non-minimally coupled scalar filed. They found that for black hole with initial masses smaller than a certain critical value, the accretion of the scalar filed can led to mass decreasing even in the absence of Hawking radiation and phantom energy. Also the black holes with initial masses greater than critical value grow by accreting the scalar filed similar to the minimally coupled scalar case.

Summary and Conclusion
In this paper, the geodesic motion and accretion process of a test particle in a subclass of the general Horndeski/Galileon gravity theories in the equatorial plane of a non-rotating black hole are investigated. In this framework, the circular geodesics, the stability of such orbits, oscillations under the action of small perturbations, unstable orbits and finally accretion process of the fluid flowing around the black hole have been investigated in a general form. Expressions for the effective potential, energy, momentum, characteristic radii, emission rate, epicyclic frequencies and dynamical parameters of the system and also the mass evolution of the black hole are derived in details. Then isothermal fluid with equation of state p = kρ is considered and some discussions are done for this subclass of solutions of the Horndeski/Galileon black holes. In this manner, the metric parameters with some approximations are obtained as h(r) ≈ C − µ r + B and f (r) ≈ h(r) B+C , where by assumption we have set C = 1. The effect of the Horndeski/Galileon correction factor B is considered for each case and our solutions are compared with the Schwarzschild black hole solutions. Our analysis has revealed that these Horndeski/Galileon solutions have deviations from the Schwarzschild solutions (which is recovered where B = 0) substantially. Our results show that Horndeski/Galileon correction factor affects the effective potential and as a result changes the loci of the stable and unstable circular orbits. For larger values of this parameter, V ef f achieves larger value and unstable circular orbits will be located at smaller radii, whereas stable orbits will be located at farther distances from the central mass. When the metric parameter B decreases, two points joint together at ISCO.
In this spacetime the location of the characteristic radii such as r isco , r ph , r sing and r mb have considerable deviation from the Schwarzschild solutions. These radii, except r mb , are decreasing functions with respect to the Horndeski/Galileon correction factor B, that is, they will be closer to the central mass for larger values of B. As r isco represents the inner edge of the accretion disk, our results show that for larger deviations, the disk will be extended close to the central mass. The behavior of r mb is different and has a minima at the Schwarzschild case. This radii is a decreasing function for B ⊂ (−1, 0) and it is a growing function for B ⊂ (0, 1.125), and finally it coincides with the innermost stable circular orbit.
As Horndeski/Galileon correction factor B grows, the energy raises whereas angular momentum decreases and one can see from the energy diagram that the range of bound orbit will be smaller for larger deviations. Increasing the parameter B enhances the E isco and then the efficiency of accretion will be decreased accordingly. In the case of Schwarzschild black hole the efficiency equals to 0.057 at r = 6r g . For negative values of the correction factor B, efficiency grows up and becomes greater in the Schwarzschild black hole limit (B = 0), where this behavior is reverse for positive values of the correction factor. The flux of the radiation energy has a maximum in the vicinity of the black hole and decreases at smaller radii. By increasing the parameter B, the flux raises and the maximum of the flux turns to the smaller radii but the reverse happens after this point. The dependence on this parameter is very considerable in the vicinity of the black hole, but it is weak far from the black hole. These behaviors are the same for temperature.
In this paper, in addition to investigation of the circular orbits and their properties, epicyclic frequencies are studied. Vertical epicyclic frequency is monotonically decreasing function of r and has no extrema. Deviation from the Schwarzschild case in this Horndeski/Galileon setup has no effect on the vertical epicyclic frequency, while the radial epicyclic frequency always has a maximum and the effect of Horndeski/Galileon correction factor on it is considerable. Dependence on this parameter is significant close to the black hole, but far from the central mass this effect is weak. Increasing the parameter B shifts the maximum to the smaller radii. We found that Ω r < Ω θ and the ratio of Ω θ Ωr is decreasing function with respect to r. In the vicinity of the black hole, this ratio is very larger than unity but far from the black hole, it turns to the unity and it decreases by increasing B. The dependence of some important resonances such as parametric and forced resonances with respect to Horndeski/Galileon parameter B is obtained, which this dependence is a monotonically decreasing function, that is, larger values of the Horndeski/Galileon correction factor lead to happening the resonance at the smaller radii.
Finally, the accretion process of the isothermal fluid is discussed for equation of state parameter k = 1 2 and the behavior of the radial velocity and density are studied. Radial velocity is a decreasing function of r as fluids have zero radial velocity in far from the black hole. When accretion happens, fluid passes at a critical point where in this point, the speed of flow matches the speed of sound. The fluid flows at sub-sonic speed before the critical point. After passing this point and in the vicinity of the black hole, because of strong gravitational field, the speed of flow increases and will be super-sonic then after. Our results show that velocity decreases by increasing B and the loci of the critical point shifts towards the black hole. Therefore, the speed of infalling particles reaches the speed of sound closer to the central mass. Finally, the rate of accretion is discussed where we found that this rate depends on the nature of the fluid and also the metric parameter. For normal fluid,Ṁ > 0 and its value is larger in the vicinity of the black hole because of strong gravitational effect and positive deviation from the Schwarzschild case increases the accretion rate.
We note that in this paper a non-spinning particle is considered. When a particle has spin, this spin has considerable influence on the particle's orbit. Also, a perfect fluid is considered and viscosity as well as magnetic field of the accretion disk are ignored for simplicity. However, these effects can affect the motion of the test particle and therefore they can affect the structure and the emission rate of the accretion disk. So, we are going to study the behavior of spinning particles and accretion of viscose fluids subject to the Horndeski/Galileon gravity in presence of a magnetic field in our future work.