Slowly rotating black holes in the Einstein-Maxwell-scalar theory

We investigate a slowly rotating black hole solution in a novel Einstein-Maxwell-scalar theory, which is prompted by the classification of general Einstein-Maxwell-scalar theories. The gyromagnetic ratio of this black hole is calculated, and it increases as the second free parameter $\beta$ increases, but decreases with the increasing parameter $\gamma\equiv \frac{2 \alpha^{2}}{1+\alpha^2}$. In the Einstein-Maxwell-dilaton (EMD) theory, the parameter $\beta$ vanishes, but the free parameter $\alpha$ governing the strength of the coupling between the dilaton and the Maxwell field remains. The gyromagnetic ratio is always less than $2$, the well-known value for a Kerr-Newman (KN) black hole as well as for a Dirac electron. Scalar hairs reduce the magnetic dipole moment in dilaton theory, resulting in a drop in the gyromagnetic ratio. However, we find that the gyromagnetic ratio of two can be realized in this Einstein-Maxwell-scalar theory by increasing $\beta$ and the charge-to-mass ratio $Q/M$ simultaneously (recall that the gyromagnetic ratio of KN black holes is independent of $Q/M$). The same situation also applies to the angular velocity of a locally non-rotating observer. Moreover, we analyze the period correction for circular orbits in terms of charge-to-mass ratio, as well as the correction of the radius of the innermost stable circular orbits. It is found the correction increases with $\beta$ but decreases with $Q/M$. Finally, the total radiative efficiency is investigated, and it can vanish once the effect of rotation is considered.


I. INTRODUCTION
Despite the great success of Einstein's general gravity (GR) in continued consistency with observations, there are compelling reasons to investigate theories of gravity beyond GR, ranging from the attempts at quantum theories of gravity [1], to the explanation of phenomena such as inflation [2,3], dark matter [4] and dark energy [5,6].
The Brans-Dicke gravity theory is one of the earliest suggested theories in the literature. It is expressed in the Jordan frame, where the scalar-tensor theory was initially developed according t o [7]. Brans-Dicke's theory describes gravity by the metric tensor and a scalar field non-minimally coupled to gravity, which O'Hanlon [8], Acharia and Hogan [9] identified as the dilaton field. Based on the concept from particle physics, the dilaton field which manifests itself as the Nambu-Goldstone boson with broken scale -invariance, may mediate a limited range gravity [10]. Indeed, by using the conformal transformations, one may extract the action of the Einstein-dilaton gravity theory from scalartensor theory in Jordan's frame. For a thorough understanding of scalar-tensor gravity and its relation to the dilaton field, we refer the readers to the monograph by Fujii and Maeda [11].
In the low-energy limit of several underlying quantum theories of gravity, the dilaton field, in conjunction with the axion field, which is deduced from string theory's low-energy limit, has produced intriguing results in inflationary cosmology and, more recently, in the acceleration of the universe [12]. Exact static dilaton black hole solutions of Einstein-Maxwell-dilaton (EMD) gravity have been constructed by many authors [13][14][15][16]. However, exact rotating dilaton black hole solutions have been obtained only for some limited values of the dilaton coupling constant [17][18][19][20], the most notable of which is the charged rotating Kerr-Sen solution, which takes into account both the dilaton and axion fields. In turn, the Kerr-Sen solution may be used to conduct an indirect test of string theory. Although the Kerr-Sen metric bears a strong resemblance to the Kerr-Newman metric, the inherent geometry of the two black holes varies significantly. The distinguished properties of the two spacetimes have been extensively studied in [21][22][23][24]. For general dilaton coupling constants, the characteristics of charged rotating dilaton black holes have been explored exclusively in the situation of infinitesimal tiny charge [25] or angular momentum [26][27][28]. The slowly rotating black hole solutions in Einstein-Maxwell-Scalar theories in asymptotically flat or AdS spacetime have been widely researched (for references, see [29][30][31]).
Recently, the EMS models have attracted attentions once again owing to the studies on black hole spontaneous scalarization [32,33], which arises from the earlier established spontaneous scalarization of neutron stars in the setting of scalar tensor theories. This kind of EMS model admits both a Reissner-Nordström (RN) solution and a scalar solution. For a sufficiently large charge-to-mass ratio, however, the RN BH becomes unstable to scalar perturbations and dynamically evolves with a scalar field profile, making scalarization energetically advantageous. Astefanesei et al. [34] proposed a classification of the BH solutions in EMS models, based on the behavior of the coupling function. They considered the action EMS models are categorized into two categories based on whether or not the field equations admit the RN BH solution (φ(r) = 0, alternatively, K ,φ (0) ≡ dK(φ) dφ φ=0 = 0). When RN BHs solve the field equations, the type is termed the scalarised-type; otherwise, the dilatonic-type. For example, the case of K(φ) = e 2φ , i.e. the EMD theory, where RN BHs don't solve the field equations, presents a specific instance of the dilatonic-type. Scalarised-type is further divided into two subclasses according to whether or not the scalar field profile is continuously connected with RN black holes (scalarised-connected-type or scalarised-disconnected type) by examining the linearization of the field equations for small φ, If the condition of µ 2 eff < 0 holds, then the scalarized BHs bifurcate from the RN BHs due to tachyonic instability, a process known as spontaneous scalarization. The scalarized BH reduces to the RN BH when φ satisfies φ = 0, and then this type is termed as the scalarised-connected-type. A particular type, for instance, is K(φ) = e 2αφ 2 .
After examining the broad categories of EMS models, namely the dilatonic and scalarised type, one finds that dilaton black holes (allowed in dilatonic models) and RN black holes (allowed in scalarised models) both play critical roles. Then what will happen if the two mutually incompatible solutions are combined? To this end, we list the line elements of these two theories in the following. The line element of a dilaton black hole for K(φ) = e 2φ is while the line element of a RN BH, i.e. K(φ) = 1, is A natural combination of the two line elements is [35] where the angular component stems from that of dilaton BHs, while the non-angular component combines those of dilaton BHs and RN BHs. β is the second dimensionless free parameter that determines the coupling of the dilaton and Maxwell field. The corresponding action is then where F µν denotes the component of the Maxwell 2-form, φ a scalar field and R the Ricci scalar. Turimov et al. [36] have investigated the geodesic of this theory. We extended it by combining RN BHs and general dilaton BHs [35] and by considering higher dimensions [37]. The present article is devoted to the study of the slowly rotating black hole in four-dimensional and asymptotically flat spacetime. In this approach, we aim to get a better understanding of this theory. The article is organized in the following way.
In Section II, we briefly review this theory and derive the equations of motion for the slowly rotating black holes. Section III presents the numerical solution to black holes. In Section IV, we make a research on the properties of the slowly rotating black holes in two subsections, A and B. In subsection A, we study the angular momentum, the angular velocity of the event horizon, and the gyromagnetic ratio of the black holes. In subsection B, we compute the correction to geodesics due to the rotating effect. Finally, in Section V, we make a discussion and summarize our results.

II. ACTION AND THE EQUATIONS OF MOTION
The scalar field φ and the vector potential A µ from action S 1 are [35] and respectively. To understand the role of coupling constant β, we display the evolution of K(φ) against φ for various β in Fig. 1 [35]. It is shown that with the increase of β, the effect of the Maxwell invariant diminishes and gravity takes over the electromagnetic interaction, allowing the electromagnetic field to be safely ignored. Apart from the situation of β = 0 (i.e. EMD theory), the potential K(φ) exhibits an extreme at φ 0 = 1 4 ln 2+β β . Thus, if we make a transformation of φ → φ + 1 2 ln(φ 0 ) on the coupling function K(φ) such that K(φ) ∝ 1 cosh(2φ) , we get a maximum for K(φ) at φ = 0. At the maximum, we have dK/dφ| φ=0 = 0, implying that RN spacetime is the solution to action S 1 . However, as can be shown, solution (5) is not identical to the RN solution, and the scalar field (7) isn't trivial. We'll show they constitute the second set of black hole solutions-the scalarised counterparts of the RN black holes according to the two Bekenstein-type identities in [34].
The first identity is given by For a purely electric field, one has F 2 < 0, which suggests should be satisfied in some region of r outside the event horizon. Otherwise, the two terms of the integrand would always have the same sign, implying that the identity holds if and only if φ = 0. The second identity is given by This reveals that for a pure electric field, the potential should satisfy the condition φK ,φ > 0 (12) in some region of r outside the event horizon. The graphs of K ,φφ and φK ,φ with K = 2e 2φ β+2+βe 4φ and Q 2 = 0.5, M = FIG. 2. The plots of K ,φφ and φK ,φ with respect to r, respectively. In some region outside the event horizon, one has K ,φφ > 0 and φK ,φ > 0. Therefore, the theory has a scalarised black hole solution with a pure electric field.
1, β = 0.2 are shown in Fig. 2. It is apparent that in some range outside the event horizon (which is denoted by the zeros of metric component g 11 ), we always have K ,φφ > 0 and φK ,φ > 0. Then, we conclude that the theory has the scalarised black hole solution.
It is straightforward to demonstrate that µ 2 eff > 0 for the electrical case. Therefore, it does not suffer the tachyonic instability problem. In other words, the RN solution is stable to the scalar perturbations. According to [34], this model is classified into scalarised-disconnected-type. It seems that the asymptotic value of scalar field in (7) is set to zero, rather than arbitrary, making the black hole solutions unnecessarily specific. Indeed, if we rescale φ → φ + 1 2 ln(φ 0 ) and A µ → β 2 + 2β to make the new K(φ) achieve extreme at φ = 0, which is different from FIG. 1, then the coupling constant β disappears and the action becomes The corresponding solution is where the physical electric charge of the black hole is given by . Then, the asymptotic value of φ is not vanishing any more, and has an effect on the metric through the relationship between β and φ 0 .
The coupling between the scalar field and the Maxwell field in S 2 is 1 cosh2φ whereas in the dilaton theory is e 2φ . This connection can be explained by Herdeiro and Oliveira's [38] broader notion of electromagnetic duality: by some transformations, two different models are related by a non-trivial duality map. Let's begin from the Einstein-Maxwellscalar class of model, whose action is given by whereF µν denotes the Hodge dual of Maxwell 2-form and f (φ) , g(φ) denote two unspecified coupling functions. The solution is described by One can then establish an electromagnetic duality transformation defined by an angle θ. D θ maps any solution (16) of a certain EMS model S 0 to a different solution of a different (dual) model, within the same family as follows In the case of which is first discussed by Gibbon and Maeda in [13] and later by Garfinkle, Horowitz and Strominger [14], after the value of θ = π/4 is taken, the model along the duality orbit has f θ = 1/ cosh 2φ, g θ = tanh 2φ and its action is Then, using equations (13) and (19), we find that S 2 differs from S ′ π 4 only in the axion term.
The generalization of S 1 is represented by the action after considering the general Einstein-Maxwell-dilaton solution [26]. The corresponding static solutions are described by the line element and two other fields It should be emphasized that γ is introduced in order to simplify the notation and it is defined by Here a and b are related to the mass M and electric charge Q of the black hole by Observing the expressions of U and f in Eq. (21) and Eq. (22), we see the metric combines the dilaton part 1 − b r 1 − a r 1−γ and RN part βQ 2 f 2 together. When β = 0, the action S 3 reduces to the model of gravity coupled to a Maxwell field and a dilaton field by Horne and Horowitz [26]. When α = 1, the action S 3 reduces to S 1 . One can again use a similar transformation like that from action S 1 to action S 2 and deduce the following action from S 3 It's obvious that the analysis above on the case of α = 1 can be extended to any α (see the analysis concerning S 1 and S 2 ). It should be stressed once again, that in the solution of action S 2 and S 4 , β does not act as a free parameter but act as a third hair representing the scalar field's non-vanishing asymptotic value.
The above study fully illustrates the motivations for the novel Einstein Maxwell-scalar theory S 3 [35][36][37], as well as the connection between it and the dilaton model. Now, let us extend our grasp of this theory by considering rotation, starting with the associated static solution of S 3 , and concentrating on the search for a slowly rotating black hole solution in asymptotically flat spacetime. Varying the action S 3 with respect to the metric, Maxwell and the scalar field, respectively, yields We can solve equations (27)- (29) to first order of the angular parameter ǫ. We further assume, in accordance with [26], that the unique term in the metric that changes to the order of O(ǫ) is g tφ , that the scalar field does not change to the order of O(ǫ), and that A ϕ is the only component of the vector potential that changes to the order of O(ǫ). As a result, we suppose the metric has the following form and the vector potential with A t in equation (23). Inserting the metric, the Maxwell fields and the scalar field into the field equations leads to the perturbation equations and where the prime denotes the derivative with respect to r.
Combing the equations of R θθ and electromagnetic field we find equations (32) and (33) can be reduced to and respectively. Integrate equation (36), and then we achieve Inserting it into equation (37), we obtain We have now obtained two main equations (38) and (39). Before we solve the two equations, we must first consider the boundary conditions, as seen in the next section.

III. METHODOLOGY FOR SOLVING THE EQUATION
Observing equations (38) and (39), one finds the value of the const hasn't been determined, and we'll see later that the angular momentum ofspacetime is proportional to const. As usual, we require the solutions have the behavior of k(r) → O(1/r) and B(r) → 1/r when r → ∞. Solving equation (38), we obtain Given that k(r) → O(1/r) at infinity, the asymptotic behavior requires that c 1 must vanish. Now we focus on the equation of B(r), i.e. equation (39). First, we search for the Frobenius series solution of B(r) at infinity. Using the reciprocal substitution r ≡ 1/y, we transform equation (39) into For the sake of calculation of the coefficients in the series, we multiply equation (41) by f 4 /K(φ) 2 , and then the equation is transformed into By expanding the coefficients of d 2 B dy 2 , dB dy and B at y = 0, one finds they are in the order of O(1), O(1/y) and O(1/y 2 ), respectively. This demonstrates y = 0 is a regular singular point of the equation. Two roots of the indicial equation of the homogeneous equation are −2 and 1, respectively. Therefore, we will omit the exponent of −2 since B(y) = O(y) is required when y → 0 . As a result, we substitute into equation (42) and the coefficient of y 0 term gives As in customary, we need c 0 = 1(required by B(r) → 1/r), but there are still two variables to be determined, namely const and c 1 . The equations of higher orders will introduce c 2 , c 3 , and so on, indicating that the equations of coefficients are not closed. When β = 0 (i.e. EMD theory), Sheykhi et al. provided a particular solution, B(r) = 1/r in [29]. However, since finding an exact solution to B(r) for an arbitrary value of β is difficult, we must resort to numerical approaches, where one should set the boundary condition first. In [39], Barausse et al. made a study on slowly rotating black holes in Einstein-aether theory where they impose the condition that the solutions are regular everywhere, except for their center (singularity of the black hole). Equation (42) exhibits apparent singularities on the horizon of the black holes as defined by U = 0, i.e.
If the solution is to be regular there, B ′ (r) and B(r) should satisfy the following equation on the horizon y h (or equivalently, 1/r h ), Combining this condition and the asymptotic behavior of B(r) at infinity, we can specify the exact value of const. However, it should be noticed that there is not one but two horizons in the general case from the quadratic equation (45). Actually, there will be two horizons (r h > a) as long as 0 < βQ 2 < (b−a) 2

4
. In general, when two horizons exist, one would anticipate that two singularities would appear in the equation. As mentioned before, to preserve the regularity on the horizon, we must implement the local regularity condition. Thus, one has to impose two local conditions in the presence of two horizons. However, by applying one regular condition alone, the solution is already specified without any tuning of other parameters in order to impose further regularity conditions. Therefore, it is sufficient to impose the regular requirement on just the outermost horizon although the solutions with many horizons will display several singularities on the horizons. This is acceptable because the outermost horizon can be rendered regular by using the usual regularity condition while the remaining singularities are hidden in the outermost horizon. Therefore, the subscript y h in (46) represents the reciprocal radius of the outermost horizon.

IV. PROPERTIES OF THE SLOWLY ROTATING BLACK HOLES
In this section, we investigate the case of β ≥ 0, where the sign of K(φ) is always positive in the whole spacetime. Otherwise, there would exist regions where K(φ) < 0, in which the Maxwell field would become a phantom field and won't be considered in this article.
A. angular momentum, gyromagnetic ratio, and angular velocity Now, let us consider the influence of β on const. We fix a = 1, b = 3 and set β = 0, 0.2, 0.4, 0.6, 0.8 each time, and then conduct research on the relationship between const and γ, which is displayed in Fig. 3 (As for the physical meaning of a and b, keep in mind that they are connected with the mass M and electric charge Q through equation (25)). The figure shows that several curves cannot be extended into γ = 0 since we require 0 < βQ 2 < (b−a) 2

4
; otherwise, there will be naked singularity which violates the cosmic censorship conjecture. Considering that γ ≡ 2α 2 α 2 +1 , one will find the range of γ is [0, 2). Once the numerical value of const is obtained, we may get a truncated series solution for B(y). As a final check, our numerical solution is compared with the series solution to B(y) in Fig. 4. Next, we calculate the angular momentum of the slowly rotating black hole. It can be calculated using the method provided by Brown and York [40]. Conserved charges such as angular momentum are defined using the surface stress tensor and Killing vector fields on the boundary of spacetime. The boundary of Σ is B with the metric σ ij , and the product of B with segments of timelike world lines orthogonal to Σ at B is denoted as 3 B. The surface stress energy tensor is then defined by which is derived from the variation of the action with respect to the metric γ ab . We decompose the boundary metric γ ab into the standard ADM form, where we have chosen the two-surface as a two-sphere, and the coordinates ϕ i are the angular variables parameterizing the hypersurface of constant r around the origin. σ is the determinant of σ ij . N and V i are the lapse and shift functions, respectively. Suppose 3 B possesses an isometry concerned with Killing vector field ξ. then, the corresponding conserved charge is defined by where n a is the normal vector of B and is tangent to 3 B. The conserved charge associated with the rotational Killing vector field ∂ ∂ϕ is the angular momentum, Expanding τ ab n a ∂ ∂ϕ b and collecting the terms of order O(ǫ), we obtain Perform the integration and now the angular momentum equals to Taking account of B(r) = O(1/r) at infinity, the angular momentum equals to −ǫ const 6 . Angular momentum can also be calculated using the Komar angular momentum and the same result obtained (for details see [41]). When β = 0, we have const = (2γ − 3)a − 3b and the angular momentum is which is the same as the value in [29]. Once we get the numerical value of const, we can calculate the gyromagnetic ratio of the black hole. One of the remarkable facts about a Kerr-Newman black hole is that it has the same gyromagnetic ratio as an electron in the Dirac theory, g = 2. Scalar fields, such as the dilaton field, modify the value of gyromagnetic ratio of the black hole, as a result, it doesn't possess the gyromagnetic ratio of g = 2 in [26]. We will now examine the effect of β on the gyromagnetic ratio g. The magnetic dipole moment for this asymptotically flat, slowly rotating black hole can be defined as Substituting M = 1 2 (b + (1 − γ)a) (25) and J = − ǫ 6 const into above equation, we obtain g = −6 (b+(1−γ)a) const . In the case of β = 0, i.e. the four dimensional dilaton black hole, we have const = (2γ − 3)a − 3b, and the gyromagnetic ratio is For a dilaton black hole, r = a denotes a curvature singularity and r = b represents an event horizon. Thus, a < b is ensured by cosmic censorship. When the denominator of the fraction of (55) is examined, it's found to be always positive, indicating that the gyromagnetic ratio has an upper limit of 2, the well-known gyromagnetic ratio of the Kerr-Newman black hole. By including scalar hair, the gyromagnetic ratio is suppressed. This has been shown in [42], which extends the Kerr-Newman black hole to include scalar hair.
Turning to the behavior of the gyromagnetic ratio of the slowly rotating black hole corresponding to S 3 , we plot the gyromagnetic ratio g versus γ (recall that γ ≡ 2α 2 1+α 2 ) in Fig. 5. From this figure, we find that the gyromagnetic ratio decreases with γ or α when the value of β is fixed. In the case of β = 0 (black line) [29], the gyromagnetic ratio decreases starting from 2, the value for the Kerr-Newman black hole. However, as seen in the figure, when β = 0, the gyromagnetic ratio can surpass 2. Taking γ = 0 as an example, we have B(r) = 1/r and in this case const = − (3a+3b) β+1 , implying that gyromagnetic ratio is 2(β + 1). This is a specific example that the non-vanishing positive β increases the gyromagnetic ratio. Based on the analysis above, the Kerr-Newman black hole's gyromagnetic ratio may be obtained in the Einstein-Maxwell-scalar black hole by simultaneously raising the values of β and γ.
For the sake of simplicity, we shall restrict ourselves to the case of γ = 1, i.e, α = 1. Then the corresponding action is S 1 . Expanding the expression of K(φ) = 2e 2φ β+2+βe 4φ in terms of β K(φ) = 2e 2φ β + 2 + βe 4φ = e 2φ − e 4φ + 1 e 2φ β 2 + e 2φ e 4φ + 1 one finds that K(φ) is actually the combination of infinite distinct dilaton couplings. In Fig. 6, we show the gyromagnetic ratio of the slowly rotating black hole in relation to Q M for various values of β . The black line denotes the usual dilaton field K(φ) = e 2φ , whereas the other lines represent the situations where β = 0. We find the upper limit for Q M with fixed β is in order to ensure that the singularity r = a of the black hole is not naked. It indicates that the gyromagnetic ratio increases with increasing β but decreases with increasing Q M , which means we can achieve the same ratio as Kerr-Newman black holes by simultaneously increasing β and charge-to-mass ratio.
Next, we will investigate the angular velocity of the black hole horizon. The coordinate angular velocity of a locally non-rotating observer is defined by Ω = − gtr grr = ǫ k(r) f (r) 2 . It is known that one of the important quantities is the angular velocity on the horizon Ω h = Ω(r = r h ), which affects the region where the super-radiation occurs in the black hole background [43][44][45].
Because the angular velocity of the horizon is equal to ǫ 1 for a slowly rotating charged Kerr black hole, we can make the angular velocity dimensionless as followsΩ h = Ω h /( ǫ Fig. 7, we exhibit the dimensionless angular velocityΩ h as the function of the charge-to-mass ratio. As can be seen, the dimensionless angular velocity on the horizon increases with increasing charge-to-mass ratio for fixed β. However, since the angular velocity decreases with β when Q M is fixed, one may obtain the same angular velocityΩ h as the Kerr-Newman black holes by simultaneously increasing the value of β and Q M .

B. The innermost stable circular orbit, the radiation efficiency and their corrections
In this section, we will focus on the circular orbits in the equatorial plane in order to investigate the geometry of the spacetime above. In the stationary and axially symmetric spacetime, one can find the equations of motion for geodesics in the form [46] with the effective potential given by where the overhead dot stands for the derivative with respect to the affine parameter, and the constants E and L z correspond to the conserved energy and the (z-component of) orbital angular momentum of the particle, respectively. For simplicity, we put the orbits on the equatorial plane. With the constraint that θ = π 2 , one finds the effective potential V eff (r) must satisfy in order that the circular orbit in the equatorial plane is stable. Solving the above equations, one obtains Then the corrections to energy, angular momentum, and period up to the first order of ǫ are given by and We don't bother showing the detailed formula for energy, angular momentum, and the orbital period since they are rather lengthy. The relative correction of the period is which is shown in Fig. 8. As can be shown, for fixed values of radius for circular orbits and fixed β, the relative correction of the period increases with increasing Q M . The correction becomes smaller and smaller with the increase of β.
The innermost stable circular orbit (ISCO) of the particle around the black hole is given by the equation V eff,rr = 0, i.e.
Substituting the expressions of E, L z , g tt , g tφ and g φφ into Eq. (66), we obtain the equation denoted by P (R, ǫ) = 0, where R denotes the radius of ISCO of the rotating black hole. Assuming R 0 is the radius of ISCO in the corresponding static spacetime, i.e. P (R 0 , 0) = 0, the correction to the radius of ISCO up to the first order of ǫ is then R 0 + ǫR 1 .
Again, we don't bother giving the expression for R 0 since it is the root of a quartic equation. We only display the radius of the innermost stable circular orbits concerning Q M for various β in Fig. 9(a). The relative correction, denoted by R1 R0 is shown in Fig. 9(b). The graphic shows that when the perturbation parameter ǫ > 0, the relative corrections to the radius of the innermost stable circular orbits are always negative, and the absolute value increases as the charge-to-mass ratio increases. On the other hand, with the increase of β, the relative corrections become smaller and smaller.  The relative correction to the dimensionless radius of the innermost stable circular orbits R1/R0 with respect to Q M for various β = 0, 0.5, 1, respectively. Now we consider the influence of β on the radiative efficiency η in the thin accretion disk model, which is defined by This quantity indicates the maximal fraction of energy being radiated when the test particle is accreted by a central black hole. The radiative efficiencies of the Schwarzschild black holes and extreme Kerr black holes are 0.057 and 0.42, respectively. From equation(61), we know the that energy of ISCO is where E 1 and E 2 are the derivatives of E(r, ǫ) with respect to the first and second variables, respectively. Since we have we can denote 1 − E(R 0 ) and −(R 1 E 1 (R 0 , 0) + E 2 (R 0 , 0)) by η 0 and η 1 , respectively. The Fig. 10(a) shows that the radiative efficiency η 0 , starting from 1 − 2 √ 2 3 ≈ 0.057, increases with Q M for a fixed value of β. When Q M is fixed, the radiative efficiency η 0 increases as β increases. Concerning the first order correction η 1 , to the radiative efficiency, we see from the Fig. 10(b) that the correction increases with regard to Q M when β is fixed. For a fixed Q M , the radiative efficiency η 1 decreases with the rise of β in general, following the same trend as the correction to R 0 .
Taking into account the zeroth order efficiency, η 0 and the first order efficiency, η 1 , we obtain the total efficiency up to the first order in Fig. 11. It's found that for small values of β and ǫ < 0, the efficiency will vanish for some value of Q M . In this case, the black hole's capacity to capture particles becomes so weak that the accreted matter around it is greatly diluted. This leads to the failure of accretion disk to create radiation due to the insufficient amount of stresses and dynamical frictions.

V. SUMMARY
In this paper, inspired by the categorization of EMS models, we investigate the Einstein-Maxwell-scalar theories S 1 and S 3 that admit both RN and dilaton solutions. The theories can also be obtained from EMD theories by an electromagnetic duality after omitting the axion term. They are classified as the scalarised-disconnected-type according to [34]. We conduct study on the slowly rotating black holes and summarize the key findings as follows.
1. The requirement of asymptotic flatness is not enough in order to determine the solution of the slowly rotating black hole. Therefore, to specify the constant of integration for any value of β, the regular condition on the outermost event horizon must be imposed.
2. We constrain the range of β by demanding that the black hole singularity is not naked and the coupling between the scalar field and the Maxwell field K(φ) is normal (not phantom). Then, we find the range of γ cannot be [0, 2) for several values of β. What really intrigues us is the gyromagnetic ratio of the rotating black hole. We know the rotating charged black hole, i.e. the Kerr-Newman black hole, is noteworthy for having the same gyromagnetic ratio of 2 as the Dirac electron. Additionally, it is well established that the dilaton black hole always has a gyromagnetic ratio smaller than 2, implying that the gyromagnetic ratio is suppressed in the presence of scalar hair. The model we investigate combines these two theories naturally and we are interested in its gyromagnetic ratio too. It's found that, when b/a and β are fixed, the gyromagnetic ratio decrease with the increase of γ. For fixed b/a and γ, the gyromagnetic ratio also increases with increasing β. Thus, by simultaneously raising β and γ, the gyromagnetic ratio of the four-dimensional Kerr-Newman black hole can be obtained. Due to the presence of the second free parameter β, the suppressing of gyromagnetic ratio by the scalar hair may be overcome.
3. We have investigated the case of α = 1 as a specific example. The evolution of gyromagnetic ratio with respect to the charge-to-mass ratio is studied. We find that the gyromagnetic ratio increases with the increase of β, and decreases with increasing Q M . In view of this point, the gyromagnetic ratio may be restored by concurrently raising β and Q M . The dimensionless angular velocity of a locally non-rotating observer on the event horizon is also studied since it is related to the phenomenon of super-radiation. We find that the same value of the angular velocity as in the Kerr-Newman black hole can be obtained by increasing Q M and β. 4. Finally, the corrections to the period of circular orbits and the radius of the innermost stable circular orbits are studied. It is found the relative correction decreases with increasing β and increases with increasing Q M . We also studied the radiative efficiency and the related correction in the thin accretion disk model. It's found that the black hole's capacity to capture particles becomes so weak that the accreted matter becomes very dilute provided that β is very tiny and the perturbation parameter ǫ is negative. As a result, radiation cannot be created owing to the lack of sufficient stresses and dynamical frictions.
In all, the novel Einstein-Maxwell-scalar theory combines the well-known RN and dilaton black hole solutions. Additionally, the study on the slowly rotating black hole indicates it can mimic the properties of Kerr or Kerr-Newman black holes by adjusting the free parameter β. It is worth mentioning that a limitation of our work is that we don't analyze the stability of the slowly rotating black hole and this should be included in future work.