Matter accretion onto the magnetically charged Euler–Heisenberg black hole with scalar hair

This paper deals with astrophysical accretion onto the magnetically charged Euler–Heisenberg black holes with scalar hair. We examine the accretion process of a variety of perfect fluids, including polytropic and isothermal fluids of the ultra-stiff, ultra-relativistic, and sub-relativistic forms, when fluid is accreting in the vicinity of the black hole. By using the Hamiltonian dynamical approach, we can find the sonic or critical points numerically for the various types of fluids that are accreting onto the black hole. Furthermore, for several types of fluids, the solution is provided in closed form, expressing phase diagram curves. We compute the mass accretion rate of a magnetically charged Euler–Heisenberg black hole with scalar hair. We observe that the maximum accretion rate is attained for small values of the black hole parameters. We may be able to understand the physical mechanism of accretion onto black holes using the outcomes of this investigation.


I. INTRODUCTION
We commence by emphasizing that one of the most unexpected events in our universe is the existence of black holes (BHs).We want to have a discussion about BH as an aspect of classical physical theory.Einstein's theory of gravity, which establishes time as well as space, is the classical theory.Black holes were first thought to exist exclusively in theory, and while their models were researched in considerable detail, many scientists, including Einstein, questioned whether or not they actually existed.How does Einstein's theory of gravity explain the concept of time and space around a huge object like a star? is a natural question to ask in this situation.Schwarzschild discovered the answer to this problem, which simply relies on the mass of all static round objects.Yet, when all of the mass is contained inside a specific radius known as the Schwarzschild radius, surprising things may happen.Then, as an event horizon occurs at the Schwarzschild radius, the term "BH" is used.
It is generally accepted that celestial objects, such as BHs, gather mass through a process known as accretion.The accretion phenomenon that surrounds enormous gravitating objects is a fundamental concept in astrophysics that plays a vital role in comprehending a number of astrophysical behaviors and speculations, such as the formation of super-massive (BHs), the expansion of stars, the emission of X-rays from compact star binaries, the luminosity of quasars, and other phenomena [1]- [3].The accretion of matter in a realistic astrophysical procedure is incredibly complex since it encompasses numerous significant problems of general relativistic mag-netohydrodynamics, such as nuclear burning, turbulence, radiation processes, etc.It is useful to summarize the challenge by establishing some assumptions or assuming some basic conditions in order to comprehend the general accretion processes.
The Bondi stationary, spherically symmetric solution [4] illustrates the fundamental accretion process by depicting an infinitely massive homogenous gas cloud steadily accreting onto a gravitational object in the center.Newtonian gravity serves as the foundation for the Bondi approach.Thereafter, Michel [5] examined the steady-state spherically symmetric flow of test fluids towards a Schwarzschild BH within the context of general relativity (GR).Subsequently, Shapiro and Teukolsky [6] also made contributions to the concept of relativistic accretion on compact objects.Furthermore, Babichev et al. [7] found that if phantom energy is permitted to accrete onto the BH throughout the accretion procedure, the BH mass may drop.Additionally, as demonstrated by Jamil et al. [8] the phantom accretion not only reduces the BH mass but also transforms it into a naked singularity.Debnath [9] elaborated the static accretion onto a general class of spherically symmetric BHs by examining the impact of the cosmological constant on the accretion rate, in accordance with the Babichev model.Bondi-type accretion onto the Reissner-Nordstrom antide-sitter spacetime was studied by Ficek [10].Using the methodology described in Ficek [10], Ahmed et al. [11] investigated the process of accretion onto the Reisner-Nordstrom anti-de-sitter BH with a global monopole.In the f (R) and f (T ) modified theories of gravity, they elaborated their prior research for accretion onto BH [12,13].Abdul Jawad and M. Umair Shahzad have investigated the accreting fluids onto regular BH using the Hamiltonian Technique [14].Matter Accretion onto a Conformal Gravity BH was determined by G. Abbas and A. Ditta [15].In [16], Astrophysical accretion near a regular Hayward BH was calculated by A. Ditta and G. Abbas.The accretion of matter onto a brane-world BH through the Hamiltonian approach has been described in Ref. [17].The research conducted by Sen Yang et al. focuses on the examination of the spherical accretion flow onto general parameterized spherically symmetric BH space-times [18].Moreover, the spherical accretion is discussed in [19,20].General relativistic dust accretion for stationary rotating BHs was established by Azreg-Ainou [21].
In literature [22]- [27] have discussed the accretion phenomena in numerous space-times.
The process of transonic accretion and the presence of the sonic point (or critical point) are significant aspects of spherical accretion onto the BH.During the accretion flow transitions occurs from a subsonic to a supersonic state at the sonic point.In a particular BH space-time, the sonic points are typically found close to the horizon.The narrow region near the sonic point is significant and fascinating, it is necessarily linked to current research on the gravitational and electromagnetic wave spectra.Consequently, investigation of the spherical accretion problem can not only help in our understanding of the accretion process in various BH but also, more significantly, convey a distinct perspective on how to investigate the nature of BH space-time under strong gravity.
This paper mainly focuses on the analysis of astrophysical accretion near the magnetically charged Euler-Heisenberg (EH) BH with scalar hair.The electrodynamics EH Lagrangian was first proposed in 1936 [28].In [29], a methodology for identifying the impact of the EH hypothesis was presented.It was only the natural way to relate the EH Lagrangian to the Ricci scalar via the volume element to investigate BH solutions because the EH theory possesses remarkable physical characteristics.In [30], analytical solutions were found for the magnetically charged situation while simultaneously discussing electric charges and dyons.This work provided the first BH solution to EH electrodynamics.In [31] and [32], electrically charged BH was addressed, whereas the geodesic structure was the focus of the research in [32].In [33] the authors have examined the charged particle motions surrounding the EH AdS BH.Whereas the quasi-normal modes were computed in [34], the thermodynamics of these BH were explored in [35,36].The BH together with the EH Lagrangian and modified gravity theories were investigated in [37]- [39].Ultimately, the shadow of EH BH was studied in [40].Recently, a BH solution is obtained in the EH theory [41].The effects of the coupling constant of EH theory on the thermodynamics and energy conditions of this BH have also been investigated [41].Also, the motion of particles in a magnetically charged EH BH with scalar hair is studied in [42].The fundamental objective of the present investigation is to address the obvious question of whether magnetically charged EH BH with scalar hair might affect astrophysical accretion processes using the Hamiltonian approach.Focusing on perfect fluid accretion onto the magnetically charged EH BH space-times, we investigate the transonic phenomena for various categories of fluid, including isothermal fluids (such as ultra-stiff, ultra-relativistic, radiation, and sub-relativistic) and polytropic fluid.
The following is the structure of our paper: A brief introduction to the magnetically charged EH BH with scalar hair is given in Sec.II.In Sec.III, we provide some helpful quantities and construct the fundamental equations for subsequent consideration of the spherical accretion of various fluids.We examine the accretion processes and determine the critical points of the system in Sec.IV.We apply the obtained formalism or findings to a number of well-known fluids and extensively investigate the transonic phenomenon for the accretion of these fluids in the magnetically charged EH BH with scalar hair in sec.V. Furthermore, we compute the BH mass accretion rate in Sec.VI by considering the impact of the accelerating parameters.Finally, we provide a summary of this article in Sec.VII.

II. THE SPACE-TIME OF MAGNETICALLY CHARGED EH BLACK HOLES WITH SCALAR HAIR
This section provides a concise overview of magnetically charged EH BH in the Einstein-Euler-Heisenberg theory.For this purpose, the Euler-Heisenberg action in the presence of a scalar field is given by [41], where L represents the Lagrangian of the Einstein-Euler-Heisenberg theory, R is the Ricci scalar, is the field strength, and ϵ µνρσ denotes the Levi-Civita tensor that fulfill The corresponding field equations can be obtained by varying the above action with respect to the spacetime metric g µν , the scalar field ϕ, and the electromagnetic field A µ , which are given respectively by where The mentioned spherically symmetric space-time metric ansatz is taken into consideration as follows by assuming the four-vector where Q m represents the magnetic charge for the BH.Also, we define the following quantities If dyons are not taken into account (both magnetic and electric charge), Q will become extinct.By using Eqs.
(3) to (7), we obtain and the metric function b 1 (r) and b(r) are given by where c 1 and c 2 represent constants of integration and ν is used to analyze the behavior of the scalar charge.At sufficiently large distances, we will impose that ν > 0. Hence, we have A new scale is introduced into the theory by the scalar charge, yielding an underlying cosmological constant.Since the generated mass term is determined by the scalar charge and an integration constant, the BH has a secondary scalar hair, so where m = ν+c2 6 and Λ eff = −(3c 1 + 12 ν 2 ).The condition for the small scalar case is achieved when ν → 0, thus Here m, Q m , α, Λ eff and ν represents the mass, magnetic charge, EH parameter, the cosmological constant, and the coupling small scalar hair charged parameter of the EH BH, respectively.By setting b = 0, one can determine the horizon of the BH.By using Eq. ( 17), we are unable to determine the horizon of BH analytically, so it is expressed graphically in Fig. 1.

III. FUNDAMENTAL EQUATIONS FOR SPHERICAL ACCRETION FLOW
In this section, we calculate the fundamental equations of accretion around the magnetically charged EH BH with small scalar hair.To analyze this, we use two fundamental laws: the conservation of the number of particles and the conservation of energy.We assume perfect fluid is flowing around the BH.The energy-momentum tensor for the perfect fluid is given by where e and p represent the energy density and pressure, respectively.If the proper number density is n, then the flux density is defined by J µ = nu µ , where u µ = dx µ dτ is the 4-velocity of the particles.In the accretion process, we presume that no particles are formed or destroyed, which means that the total number of particles is conserved, so particle conservation and energy conservation are given as follows By solving Eq. ( 21) in equatorial plane (θ = π 2 ), we have where C 3 is an integration constant.Since we consider the flow of the fluid to be in a radial direction in the equatorial plane, only two components u t and u r = u are different from zero.By utilizing the normalization condition, we attain with F being given by and u t takes the form Moreover, the first law of thermodynamics for a perfect fluid is stated as [43] dp = n(dh − T ds), de = hdn + nT ds, where s represents entropy, T indicates the temperature, and h denotes the specific enthalpy which is denied by There is a scalar hu µ ξ µ in relativistic hydrodynamics that is conserved along the fluid's trajectories [43], so where ξ µ stands for the Killing vector of space-time.If we assume ξ µ = (1, 0, 0, 0), we attain Integrating the above equation, we have where C 4 is integration constant.It is simple to demonstrate that the specific entropy of the fluid is conserved across the flow lines u µ ∇ µ s = 0.If we rewrite T µν as nhu µ u ν + (nh − e)g µν and then apply the conservation equation of T µν onto u µ , we get In the specific scenario, we assume that the motion of the fluid is radial, static (it does not change over time), and it conserves the BH spherical symmetry, so the above equation reduces to ∂ r s = 0, which means that s is constant everywhere.As a result, the fluid's motion is isentropic and Eq. ( 27) becomes we will study the flow by using Eqs.( 23), ( 31) and (33).Also, it transforms the equation of state (EOS) of a simple fluid, e = e(n, s), from its canonical form to its barotropic form because s is a constant, so by using second Eq. ( 33), we obtain h = de dn , which gives us here, ′ indicates the derivative with respect to n.Also, p ′ = nh ′ is produced by the first equation of Eq. ( 33), when h = F ′ (n), we get by integrating Eq. ( 36), we have We know that the EOS of the form p = G(n) cannot exist independently of the EOS of the form e = F (n).By solving the above differential equation, the connection between F and G can be determined.
We can determine the sound speed in a local inertial frame by using the following formula a 2 = ( ∂p ∂e ) s , which is given in [44].Because entropy s is constant, we can reduce this to a 2 = dp/de.From Eq. ( 33), we determine the useful formula for the subsequent parts as follows utilizing Eq. ( 33), in Eq. ( 39), we have Another useful expression is the three-dimensional fluid velocity ω, determined by a local stationary observer.As the motion in the equatorial plane is radial, so dθ = dϕ = 0, and then Eq. ( 8) becomes In the usual relativistic method [45,46] as seen by a local, stationary observer, the typical three-dimensional velocity ω can be defined as which yields By using u t = −F(r)u t , u r = u = dr dτ , u t = dt dτ , and Eq. ( 26), we have and Then employing Eq. ( 23), we have These equations will be utilized in the following Hamiltonian evaluation.

IV. DYNAMICAL SYSTEM AND SONIC POINTS
According to the fundamental Eqs. ( 23) and ( 31), there are integration constants C 3 and C 4 .Also, Hamiltonian H is read as the square of the L.H.S side of Eq. (31), which is given by Here, we establish the Hamiltonian dynamical system as a function of (r, ω) to study the Michal flow, which is found in Chaverra et al. [47,48], Ahmed et al. [11]- [13], and is written in the given scenario as Furthermore, the dynamical system corresponding to the Hamiltonian is expressed as follows where the t-derivative is represented by dot.According to Eq. ( 49), H ,ω is the partial derivative of H with respect to ω when r is assumed to be constant, and H ,r stands for the partial derivative of H with respect to r when ω is taken to be constant.Ultimately, the system (49), reduces to and We set Eq. ( 50) and ( 51) equal to zero and solve simultaneously to find the critical points, which are given below where r c , a 2 c and ω c represent the distance, speed of sound, and three-velocity of the fluid at the critical point also, we assume a 2 c = k.Furthermore, we can use Eq. ( 46), to find the constant C 2 3 in terms of the critical points, we have By using Eqs.( 46) and ( 52), the result would be If the solution of Eqs. ( 50) and ( 51), do not exist at the sonic point, one can define the reference points (r 0 , ω 0 ) through the phase depiction to arrive at [13] The previous equations will also be used to investigate the spherical accretion in several fluids.

V. APPLICATIONS TO TEST FLUIDS
In this section, we will analyze different types of fluids that flow around the magnetically charged EH BH with scalar hair using the above results.Conclusively, we assume the polytropic and isothermal test fluids given in the following subsections.

A. Isothermal test fluid
The motion of fluid at a constant temperature is regarded as an isothermal flow.In other words, the speed of sound is constant in the process of accretion.This makes sure that at any radius, the speed of sound at the critical point is identical to the speed of sound of the accretion flow.Hence, our system is adiabatically operating in this situation; it is more conceivable that our fluid is flowing in an isothermal fashion.Because of this, we drive the general solution of isothermal EOS in this subsection with the form p = ke.Additionally, by using Eqs.(34) and (38), in EOS, we have p = kF (n) and G(n) = kF , where k denotes the state parameter and 0 < k ≤ 1 [49].The definition of the adiabatic sound speed is typically given as a = dp/de.Hence, when we compare the adiabatic sound speed to the equation of state, we get a 2 = k.From Eq. ( 38), we attain which yields By using p = ke and Eq.(57), in Eq. ( 29), we get From Eq. (54), and Eq.(58) we have and where all constant components are integrated into the characterization of time t and Hamiltonian H. Now, we will study how the fluid acts by using various values for the state parameter k.
1. Solution for ultra-stiff fluid (k = 1) Now, we assume the ultra-stiff fluid, which is obtained by setting k = 1 and p = ke as the equation of state.This fluid has the characteristic that its energy density and isotropic pressure are both equal.The Hamiltonian Eq. ( 60), can be transformed into the form From Eq. ( 61), it can be seen that the flow will be physical if | ω |< 1.As a result, for the case of the ultrastiff fluid, the minimum value of Hamiltonian (61), is Differentiating Eq. ( 61), with respect to ω and r, we get the following system of equation From Eqs. ( 62) and (63), it is clear that the above dynamical system does not have a critical point.The minimum value of H is The curves in between the two magenta curves in Fig. 2, represent physical flows.Curves that lie on the upper half-plane with ω > 0 show the outer flow of the fluid or the emission of particles, while curves on the lower half-plane with ω < 0 show the fluid accretion.
2. Solution for ultra-relativistic fluid (k = 1 2 ) Now we study the ultra-relativistic fluid with the equation of state parameter k = 1/2.In this scenario, the energy density of the fluids is greater than their isotropic pressure.By assuming the value of state function k = 1 4 , Eq (52) gives us We can not find the roots of the above polynomial analytically; its roots will be computed numerically.For this, we consider the following values α Further, the system of equations in Eqs. ( 50) and ( 51), takes the following form Also, we can find critical points for ultra-relativistic fluid, when the RHS of Eqs.(66) and (67), vanish.By considering the BH parameters m = 1, Q m = 1, α = 1, Λ eff = −1 and ν = 0.2, we can obtain the physical critical points (r c , ±ω c ), which are (1.13718,−0.707107) for the outflow of the fluid and (1.13718, 0.707107) for fluid accretion, respectively.By substituting these critical points into Eq.(65), we obtain the critical Hamiltonian H c = 0.8849.The critical points r c , ω c , and H c for various values of the BH parameter are listed in tables I and II.In Fig. 3, we represent the physical behavior of an ultra-relativistic fluid by numerous curves with BH parameters m = 1, Q m = 1, α = 1, Λ eff = −1 and ν = 0.2.From Fig. 3, we observe that the critical points (r c , ω c ) and (r c , −ω c ) are the saddle points of the given dynamical system.In Fig. 3, the red curve (with H = H c + 0.2 branches) and the green curve (with H = H c + 0.4 branches) show purely supersonic outer flow (ω > ω c branches), supersonic accretion (ω < −ω c ), or purely subsonic accretion followed by subsonic outer flow (−ω c < ω < ω c ), respectively.The black (with H = H c -0.2 branches) and magenta (with H = H c − 0.3 branches) curves represent no physical behavior.
In Fig. 3, the blue curves represent the fascinating solution of the fluid and reveal the transonic behavior of the fluid outside the BH horizon.For ω < 0, the curves that traverse across the sonic point (r c , ω c ).One solution be- gins at spatial infinity with subsonic flow and proceeds to supersonic flow after crossing the sonic point.This solution refers to the standard non-relativistic accretion investigated by Bondi [4].In accordance with the investigation [13], the alternative solution that proceeds at spatial infinity with the supersonic flow but transforms to subsonic after crossing the sonic point is unstable, so it is extremely challenging to find these behaviors.If ω > 0, there are two possible solutions.One solution is detailed in [4] for non-relativistic accretion, which correlates to the transonic solution of the stellar wind, initiating flow at the horizon with supersonic flow and switching to subsonic flow after crossing the sonic point, while the remaining solution is identical to the case ω < 0, which is hard to attain and unstable [13].Generally, various Hamiltonian values correspond to various initial states of the dynamical system.If the ultra-relativistic fluid has a transonic solution, then the Hamiltonian can be examined at the sonic point.The Hamiltonian with dis- tinct values from the transonic one cannot represent any transonic flow solutions.For instance, the magenta curve displays the subcritical flow of the fluid since flows will not reach the critical point.In reality, the solutions have a turning point, which is the closest point that such fluids can reach before rebounding back or turning around infinity.Similarly, the black curves can be explained.Moreover, super-critical flows can be seen in the green and red curves.Although fluids do not reach critical points, their velocities are already higher than the permissible critical value.Such flows finally enter the BH horizon.It is also worth noting that the same evaluation applies to other fluids, such as radiation, sub-relativistic and polytropic fluids.The fluid that absorbs radiation emitted from the BH is said to be radiation fluid.By considering the value of state function k = 1 3 , Eq (52) gives us We can not find the roots of the above polynomial analytically; its roots will be computed numerically.For this, we consider the following values α = 1 m = 1 ν = 0.2 Λ eff = −1 Q m = 2 of the BH parameters and corresponding real roots are r c1 = 0.702594, r c2 = 1.27678.The Hamiltonian (60), becomes and the above system (50,51) reduce to We can determine the critical points of the above dynamical system by setting the right-hand side of Eqs.(r c , ω c ) seem to be saddle points of the above dynamical system.Furthermore, from Fig. 4, it is easy to observe that the motion of the radiation fluid (k = 1 3 ) is identical to that of the ultra-relativistic fluid, which is depicted in Fig. 3. Here, red and green contours represent supersonic flows when ω < −ω c or ω > ω c , and subsonic flows when −ω < ω c < ω, while the same behavior is observed in Fig. 3, for ultra-stiff fluid.The blue curves illustrate the transonic solutions.For ω < 0, one of the blue curves (which begins at spatial infinity with subsonic flow followed by supersonic after crossing the sonic point (r c , ω c ) shows the usual transonic accretion, while the other blue curve demonstrates the solution is unstable.For ω > 0, one blue curve represents the transonic outer flow of wind, while the other indicates that the flow is unstable, analogous to ultra-relativistic fluid.The magenta and black curves demonstrate the non-physical solution.We assume that the state equation for a sub-relativistic fluid is p = e/4, which means that the energy density is higher than its isotropic pressure.By considering the value of state function k = 1 4 , Eq (52) gives us To investigate heteroclinic flow, we compute the roots of the above polynomial.But we are unable to determine the roots of the above polynomial analytically; its roots will be computed numerically.For this, we consider the following values of the BH parameters and corresponding real roots are r c1 = 1.25842, r c2 = 2.6262.Â The Hamiltonian (60) becomes and dynamical two-dimensional system (50, 51) read as Figure 5 shows the phase space profiles for the subrelativistic fluid.From Fig. 5, it can be seen that the motion of the sub-relativistic fluid (k = 1 4 ) is identical to the motion of the radiation (k = 1 3 ) and ultra-relativistic fluids (k = 1  2 ).The green and red curves are purely supersonic outer flows for ω > ω c , while they represent supersonic accretions for ω < −ω c .Also these curves display subsonic flows when −ω c < ω < ω c .The blue curves in Fig. 5, are fascinating because they depict the transonic solution of outer flow for ω > 0 and the spherical accretion for ω < 0. The magenta and black curves are un-physical solutions, just like the radiation fluid and ultra-relativistic fluid.In accordance with the previous two cases, we provide the values of r c , ω c , and H c for such a dynamical system in tables V and VI with various values of the BH parameters.Figure 6 is the contour of H (60) in (r, ω) plane.From Fig. 6, we observe that (r c1 , ±ω c ) are the two saddle points.We assume (r rm , −ω c ), (r rm , −ω c ) are the two rightmost points in the lower and upper branches of the given plot, respectively.If we consider dω dr continuous, as the fluid passes through the saddle point, the processes of accretion will begin from (r rm , −ω c ) on the Black curve in the lower curve of plot 6.If the mobility is subsonic, it flows into the upper branch of the lower plot, passes through saddle point (r c1 , ω c ), and then crosses the horizon.Alternatively, if the motion is supersonic, it travels across the lower branch of the lower plot, eventually throughout the saddle point (r c1 , ω c ), when more till ω vanishes as the fluid nears the horizon.From there, the fluid travels along the upper branch of the upper plot while traveling at supersonic speeds until the point (r rm , ω c ) Â of the upper plot.First, it may be demonstrated that such motion is unstable using considerations similar to those provided in the case when k = 1/2.Secondly, the movement might become periodic, but this is extremely difficult to achieve.

B. Polytropic test fluid
The equation of state for polytropic test fluid is characterized by here, γ and K are taken to be constants.The constraint γ > 1 is typically used when working with ordinary matter.We acquire the subsequent equations for the specific  enthalpy using [12] in which the baryonic mass m is known integration constant.The speed of sound is determined by From Eqs. ( 54) and (78), we have where From Eq. (79), Z is positive constant.By using Eqs.( 79) and ( 48), we obtain here, m 2 is consumed in the re-definition of ( t, H).As the Hamiltonian is constant on the solution curve, so no global solutions exist.Using the process described in [12][13][14], the following relationship can be obtained In order to prevent the Hamiltonian (81), from diverging, the solution curve must not cross the r-axis at the point where ω = 0 and r = r h .The dynamical system of Eqs. ( 82) and (83), has numerical solutions, which are depicted in Fig 7 .One may observe that there is only one critical point, also referred to as a saddle point, in the accretion of a polytropic fluid.Furthermore, as shown in Fig. 7, the mobility of polytropic fluids is identical to that of isothermal fluids with k = 1/2 (see Fig. 3), k = 1/3 (see Fig. 4), and k = 1/4 (see Fig. 5).

VI. BLACK HOLE MASS ACCRETION RATE
In this section, we are interested in determining the mass accretion rate of the magnetically charged EH with scalar hair BH.We are aware that when matter accretes in the vicinity of compact objects, its mass varies with the passage of time.The given relationship Ṁ = − T 1 0 dS, where dS = √ −gdθdϕ, can be used to determine the rate of change of mass of the BH over time.To calculate the mass accretion rate, we can use the general formula, which can be found as [9] Ṁ = 4πLM 2 (e + p), By using Eqs.( 21) and ( 22), one can get The relativistic flux equation can be used in the above equation to calculate where K 1 denotes the integration constant.Inserting p = ke in Eq. ( 86), we obtain By using Eqs.(58), (87) and (85), we have By solving the above equation, we can analytically determine u r for specific values of k.For ultra-stiff fluid Eq. ( 88), reduce to From Eqs. (89),(87) and (84), we have where Ṁ is the mass accretion rate for ultra-stiff fluid of magnetically charged EH with scalar hair as well as K 0 and K 1 are constants.Similarly, if we set k = 1/2, we obtain the following expression The graph of the mass accretion rate is depicted in Fig. 8. From Fig. 8, it can be seen that as the values of the BH parameters increase, the mass accretion rate decreases.

VII. CONCLUSION
This paper studies the spherical accretion flow of a perfect fluid around the magnetically charged Euler-Heisenberg BH with scalar hair.We establish two basic formulas for the analysis of the accretion processes using energy and particle conservation equations.Then, utilizing these fundamental equations, we explore the various subcategories of perfect fluids, such as isothermal fluids of ultra-stiff, ultra-relativistic, radiation, sub-relativistic, and perfect fluids.Figures 2, 3, 4, 5,6 and 7 provide a graphical representation of the behavior of the fluids flow near the magnetically charged EH BH.It is essential to understand that for ultra-stiff fluids, the sonic point cannot be observed.The critical values for ultrarelativistic, radiation, and sub-relativistic fluids are given in tables.We employed the equation of state to analyze the physical properties of matter, including ultra-stiff, ultra-relativistic, radiation, and sub-relativistic fluids, in order to discover the nature of BH.In addition to the equation of state and model parameters, different types of accreting fluid exhibit varying accretion characteristics, including subsonic, supersonic, and transonic.We noticed precisely how supersonic accretion, subsequently followed by subsonic accretion, terminated inside the BH horizon.This relates to the fluid movement around the magnetically charged Euler-Heisenberg BH with scalar hair.The mass accretion rate of BH has been determined for ultra-stiff fluid (k = 1) and ultra-relativistic fluid (k = 1/2).It also graphically investigated how the BH accelerating parameter affects the spherical mass accretion rate of perfect fluid onto gravitational bodies.From Fig. 8, we observe that the mass accretion rate attains its maximum value for small radii and then decreases to its minimum for large radii.Furthermore, we can see that the mass accretion and BH parameters have an inverse relationship.The findings presented in the previous section show that the mass accretion rate de- pends on the accelerating parameter and that this influence only becomes noticeable when the parameter's value is small.We are unable to calculate the mass accretion rate for other values of the state parameter analytically because Eq. (88), becomes highly non-linear in terms of u, so it is not possible to find an explicit form of Ṁ .
Here we would like to mention that the study performed in this paper mainly focuses on the accretion properties of different fluids onto the Euler-Heisenberg BH with the scalar hair.It is interesting to extend the current work to other related astrophysical processes, such as the black hole images with different accretion processes or thin accretion disks, and explore how the scalar hair can affect the size and shape of the corresponding image.In particular, together with the recent observations of the black hole images by the Event Horizon Telescope and inspired by [50,51], we may run a systematic study in this direction by constraining the parameters of EH BH with scalar hair in our next step.We expect to return to the above issues soon in future work.

FIG. 2 .
FIG.2.The profile H (62) for an ultra-stiff fluid is as follows: for magnetically charged EEH BH with parameters m = 1, ν = 0.2, Qm = 0.2, Λ eff = −1, α = 1.The horizon can be located at r h ≃ 1.1723.The blue curve is an illustration of H = Hmin ≃ 0.38.The green and red curves refer to H > Hmin whereas the black and magenta curves represent the behavior for H < Hmin .
of the BH parameters and corresponding real roots are r c1 = 0.692258, r c2 = 1.16033.Â The Hamiltonian (60) has the form

4 ,
(70) and (71), equal to zero and then solving for r and ω.Moreover, in tables III and IV, the critical values r c , ω c , and H c are shown for the different values of the BH parameters.In Fig. 4, several curves illustrate the physical behavior of radiation fluid with BH parameters m = 1, Q m = 1, α = 1, Λ eff = −1 and ν = 0.2.From Fig. it is clear that the critical points (r c , −ω c ) and

TABLE I .
The values of ωc, rc and Hc at critical points with various values of BH parameters for sub-relativistic fluid are given

TABLE II .
The values of ωc, rc and Hc at critical point with several values of BH parameters for radiation fluid are given.

TABLE V .
The values of ωc, rc and Hc at critical point with various values of BH parameters for radiation fluid are given

TABLE VI .
The values of ωc, rc and Hc at critical point with several values of BH parameters for radiation fluid are given.