Accretion disc around black hole in Einstein-SU(N) non-linear sigma model

The accretion of matter onto celestial bodies like black holes and neutron stars is a natural phenomenon that releases up to 40%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$40\%$$\end{document} of the matter’s rest-mass energy, which is considered a source of radiation. In active galactic nuclei and X-ray binaries, huge luminosities are observed as a result of accretion. Using isothermal fluid, we examine the accretion and geodesic motion of particles in the vicinity of a spherically symmetric black hole spacetime in the Einstein-SU(N) non-linear sigma model. In the accretion process, the disk-like structure is produced by the geodesic motion of particles near the black hole. We determine the innermost stable circular orbit, energy flux, radiation temperature, and radioactive efficiency numerically. In the equatorial plane, we investigate the mobility of particles with stabilities that form circular orbits. We examine perturbations of a test particle by using restoring forces and particle oscillations in the vicinity of the black hole. We analyze the maximum accretion rate and critical flow of the fluid. Our findings demonstrate how parameter N influences the circular motion of a test particle as well as the maximum accretion rate of the black hole in the Einstein-SU(N) non-linear sigma model.


I. INTRODUCTION
The theory of general relativity (GR) predicts the existence of black holes (BHs), which are remarkably, extremely massive objects.These are the source of the strongest gravitational field in the universe and have magnetic intensity and high spin.Because of these properties, BHs are an ideal laboratory for astrophysical studies of matter and gravity.Some observational data has confirmed the existence of BH in recent years.The first was the discovery of gravitational waves from a binary BH merger by the collaboration of LIGO and Virgo [1].A further important contribution is the Event Horizon Telescope's use of very long Baseline Interferometry to get the first images of M87 BH shadow [2,3] and the recently released image of Sgr A * [4].These significant accomplishments help us to understand GR theory and the nature of accretion discs around supermassive celestial bodies in the strong gravity regimes around the event horizon.
It is speculated that celestial bodies, such as BHs, accumulate mass through the accretion process.These might be used to test modified theories of gravity.The presence of an accretion disc is also a necessary component to maintain the high accretion rate.By slowly spiraling into a centrally compact object, diffuse material forms an accretion disc and radiates gravitational energy.The phenomenon by which a dense object like a BH gathers particles from a nearby fluid is called accretion.These particles must pass through a critical point when they begin to accelerate, where the fluid velocity equals the sound speed.At supersonic speeds, the fluid flows onto the central mass.As a result of this, the BH mass should be increasing [5].Also, this phenomenon released energy, which could be the source of some astronomical events like the formation of jets, quasars, and radiation [6].As a result of analyzing the geodesic structure of particles in the vicinity of the BH, it is interesting to analyze several typical radii, such as marginally bound orbit (r mb ) and innermost stable circular orbit (ISCO).These radii are significant in the investigation of BH accretion discs.
In accretion discs, the interior edge of the disc corresponds to ISCO, and the efficiency of the energy emitted is a representation of how well rest-mass energy is converted into radiative energy [7], which may be calculated using these radii.The maximum or minimum of the effective potential is associated with the positions of unstable or stable circular orbits, respectively.According to Newtonian theory, the ISCO can have an arbitrary radius after the effective potential reaches a minimum value with any value of the angular momentum, indicating that there is no minimum radius for the ISCO [8].However, the situation differs once the effective potential demonstrates general relativistic effects or has a significant feature that depends on the angular momentum of particles and other parameters.For instance, the effective potential in GR and for particles revolving around the Schwarzschild BH has two extremal values for any value of the angular momentum (minimum or maximum).Consequently, the two points exclusively correspond to a particular value of angular momentum.This point examines ISCO at r = 3r g [8,9], whereas Schwarzschild radius is represented by r g .
The characteristics of spacetime affect the positions of these radii in distinct metrics, and several quantities such as angular velocity, angular momentum, and specific energy are significant in the location of these points.Because of this, the study of the particle dynamics and accretion process near black holes with different black hole spacetimes has attracted a lot of attention.Ruffini et al. [10] and Bardeen et al. [11] have investigated the impacts of ISCO near the Kerr BH.In their textbook on GR, Hobson et al. [12] have established these characteristics.The efficiency of Kerr and Schwarzschild BH accretion discs was calculated by Thorne and Novikov [13].Psaltis and Johannsen [14] proposed the Kerr-like metric, and Johannsen [15] developed the accretion discs near such BHs.Tursunov et al. [16] discovered the geodesic structure and circular orbits of charged particles around rotating BHs that are weakly magnetized.Because particles in an accretion disc rotate in stable orbits, whenever perturbations affect the particles, oscillations are obtained in vertical and radial directions with epicyclic frequencies.As a result, the analysis of orbital and epicyclic frequencies is essential to the physics of the accretion discs surrounding BHs.Isper [17,18], Wagoner [19], Kato [20], and Ortega-Rodriges et al. [21] have conducted investigations in this area.Tretyakova [22] has discussed geodesics for scalar-tensor empirical characteristics in Horndeski BH.Additionally, Salahshoor and Nozari [23] studied the circular orbits and accretion discs in a class of Horndeski/Galileon BHs using a similar approach.Moreover, in [24], Abbas and Ditta have investigated the circular accretion near a regular phantom BH.Recently, the particle dynamics and thin accretion disk around black holes in the Einstein-AEther theory have also been explored [25][26][27][28].
In this paper, we investigate the accretion disc of a black hole space-time in the Einstein-SU (N ) non-linear sigma model (NLSM).The non-linear sigma model is one of the most important theoretical frameworks with applications in many contexts, for example in the dynamics of Goldstone bosons [29], in condensed matter systems [30], in string theory [30], and in supergravity [31].One relevant example is the SU (N )-NLSM with N = 2, which provides an effective description of the low-energy dynamics of pions in 3 + 1 dimensions [32].Thus, it is interesting to explore the coupling of the SU (N )-NLSM to GR both from the theoretical and the phenomenological point of view.With this motivation, the Einstein-NLSM system has been studied in a lot of works [33][34][35][36].Recently, a black hole solution is obtained in Einstein-SU (N )-NLSM theory [37].The effects of the coupling constant of SU (N )-NLSM theory on the thermodynamics and stability of this black hole have also been explored [37].It is natural to ask whether the SU (N )-NLSM model can have any effects on particle dynamics and the accretion process near the black hole, and this is exactly the main purpose of the current paper.
To study the particle dynamics of the black hole in Einstein-SU (N )-NLSM theory, we limit our investigation to the polar coordinates system and equatorial plane for the conserved quantities of the test particles.Furthermore, we determine the event horizon and singularity of the BH.Moreover, we study circular orbits and effective potential to investigate the locations of circular orbits such as the photon sphere (r ph ), the r mb , and r isco .Ultimately, the critical accretion is determined by using the isothermal fluid parameters.This paper is arranged as follows.In section 2, we provide a brief introduction of the Einstein-SU (N )-NLSM theory and its black hole solution obtained in [37].In section 3, the general formulation for the mobility of particles is discussed, where circular motion, radiant energy flux, stable circular orbits, and oscillations are explored.The general formulas for several dynamical parameters, including accretion rate, critical flow speed, and accretion for an isothermal fluid, are determined in section 4. In section 5, we analyzed the solution of the black hole in the Einstein-SU (N )-NLSM and a circular geodesic in the equatorial plane.Finally, section 6 addresses conclusions and discussions.

II. BLACK HOLE SPACETIME IN EINSTEIN-SU (N ) NONLINEAR SIGMA MODEL
In this section, we present a brief introduction to the black hole in Einstein-SU (N )-NLSM.The action of the Einstein-SU (N )-NLSM is given by [37] where R is the Ricci scalar and L µ is Maurer-Cartan form components given by with U (x) ∈ SU (N ).Here, N denotes the flavor number, and t i represents the SU (N ) Lie group generators, where i = 1,... (N 2 − 1).Hence, the gravitational constant is represented by κ, and K is an experimentally determined positive coupling constant.Varying the action in Eq.
(1) with respect to the fields U and g µν yields the field equations, where ∇ µ is Levi-Civita covariant derivative, G µν = R µν− 1 2 gµν R is the Einstein tensor, and T µν is the energymomentum tensor of the SU (N )-NLSM which is given by Note that Eq. ( 3) contains (N 2 − 1) non-linear coupled second order differential equations.The spherically symmetric BH spacetime was derived in [37] and the metric is given as with where m denotes the constant of integration and a N is given by The given configuration clearly illustrates by pionic matter and generalization of SU (N ) BH obtained in Refs.[38] and [39].Furthermore, this solution is an asymptotically anti-de Sitter (AdS) form of the Barriola-Vilenkin metric [40], and one can recover Schwarzschild-AdS BH by setting K = 0.
It is interesting to note that this solution also reduces to the BH solution in the Einstein-SU (2)-Skyrme theory with N = 2 and without the Skyrme coupling.When the Skyrme coupling is present, it gives an explicit contribution to the metric function f (r) of order 1/r 2 .The effects of such Skyrme coupling on the deflection angle of light around the BH in the Einstein-SU (2)-Skyrme theory has been studied in [41].Thus it is natural to seek how the above solution can be modified by including a Skyrme coupling in the Einstein-SU(N)-NLSM model.As mentioned in [37] and using a similar formalism, the above solution can be generalized to solutions of the SU(N)-Skyrme model and an analytical form of it has not been reported yet.Since the main purpose of the current paper is not intended to construct a solution with the Skyrme coupling, we will only focus on the above solution without the contribution from the Skyme term.
The event horizon of the BH is located at . It is observed that, in order to get a real root in the earlier result, the constant of integration m must fulfill the relation m ≥ (Kκa N − 1) The above relationship establishes the event horizon's minimum radius, namely The horizon radius is depicted in Fig. 1 for distinct values of N .

III. GENERAL FORMALISM FOR TEST PARTICLE'S MOTION
In this portion, we derive the general formalism for test particles by considering the spherically symmetric and static space-time.We consider two killing vectors ξ t = ∂ t and ξ φ = ∂ φ along with two constants of motion L and E (angular momentum and conserved energy), which are given as follows and where the four-velocity for test particle is u µ = dx µ dτ = (u t , u r , u θ , u φ ) and satisfy the normalization condition u µ u µ = 1, which can be simplified to In equatorial plane (i.e.θ = π 2 ), by using Eqs.( 12), ( 13) and ( 14), we get From Eq.( 18), we have where V ef f is the effective potential for the motion of the particle, which can be expressed as It is obvious that the effective potential depends on spacetime metric function f (r), angular momentum of particles, and radial velocity.The evaluation of effective potential plays an important role in the geodesic motion of particles because the location of ISCO is obtained by local minima of the effective potential.

A. Circular Motion of test particles
In this subsection, we examine the circular motion of particles in the equatorial plane, so u r = ur = 0 must be fulfilled because the radial component r is necessarily constant.From Eq. ( 19), we have V ef f = E 2 and also d dr V ef f = 0. We figure out angular velocity Ω φ , specific energy E, specific angular momentum L and angular momentum l, which are given by From Eqs. ( 22) and ( 23), it is clear that the specific energy E and specific angular moment L must be real if The above expression is mandatory for the circular orbits because we can figure out the particular area of circular orbits to find the solution to inequality.For bound orbits and marginally bound orbits, Eq. ( 22) must fulfil the conditions E 2 < 1 and E 2 = 1, respectively.From Eq. ( 22), we get We can acquire marginally bound orbits by simplifying the above equation.From Eqs. ( 22) and ( 23), it can be observed that specific energy E and specific angular momentum L diverge if the following relation is satisfied By using Eq. ( 27), we can calculate the radius of the photon sphere, which plays a vital role to examine gravitational lensing.

B. Radiant Energy Flux and Circular Orbits
The existence of stable circular orbits relies on the local minima of the effective potential that is obtained from > 0. By using Eq. ( 20), we obtain By using the the conditions V ef f = 0, = 0, one can determine the radius of ISCO r isco .If r < r isco , then the accretion process is feasible.When particles fall from rest to infinity and fall into a black hole, these would release gravitational energy, which can be converted into radiations by the source of the powerful astronomical events.In [42], the energy flux radiant through the accretion disc can be represented in terms of Ω φ angular velocity, the specific energy E, and the specific angular momentum L as follows Here, K is the radiant flux, Ṁ is the mass accretion rate, Ω φ , r ≡ dΩ φ dr and g is the determinant of g µν given by where we put sin θ = sin π 2 = 1 as we are dealing in equatorial plane.By using Eqs.(21)(22)(23), we have with The associated expression K(r) = σT 4 (r) contains energy flux as well as temperature because the accretion disc assumes to be in thermal equilibrium, so the radiation emitted from the disk surface is assumed as black body radiation.By assuming thermal black body radiation, it is quite easy to attain the temperature distribution on the disc using this relation, also it helps to obtain luminosity L(ν) of the disc.The disc inclination angle γ along with d distance observed by the luminosity of the accretion disc [43] is given by In the above expression, thermal energy flux is expressed by I(ν).In the mass accretion process, the efficiency of accreting material is one of the most essential factors.
We can obtain the maximum efficiency, η * by the ratio of the specific binding energy of the ISCO to the specific rest mass energy illustrated given by where E represents the specific energy of a particle revolving in ISCO.This relationship is true if all photons emitted are capable of escaping to infinity.If a fluid element is perturbed, the motion of a test particle will be in a circular orbit in the equatorial plane.In the following section, we will investigate the problem of the particles oscillating in circular orbits with three components of motion.

C. Oscillations
By virtue of restoring forces, several forms of oscillatory motion are observed in the accretion process.Horizontal and vertical oscillations are produced when restoring forces affect perturbation in the accretion discs.A few of the restoring forces in accretion discs are caused by the rotation of the disc as well as a vertical gravitational field.Due to a restoring force carried by the fluid's rotation, when a fluid element is transferred in the radial direction, it will eventually restore its equilibrium state.
For the sake of central objects, the gravitational force in accretion discs counter balances the centrifugal force.The flow element would be dragged within or outside the original radius with epicyclic frequency, Ω r , once the former prevails compared to the latter or when the reverse happens.If the fluid element is perturbed vertically in the equatorial plane, the gravitational field squeezes the perturbed element into the equilibrium position.As a consequence of restoring force, the fluid element produces harmonic oscillations across the equatorial plane with vertical epicyclic frequency Ω θ .In an accretion disc, the movement of the particles in the fluid depends on circular motion with orbital frequency, harmonic radial motion with radial frequency, and harmonic vertical motion with vertical frequency.Thus, we discuss the analysis of vertical and radial motion on nearby circular orbits in the equatorial plane.Now we elaborate the radial motion ef f of the particles.From Eq. ( 14), we take u θ = 0 and u r = 0 for vertical and radial motion.By taking u r = dr dτ = dr dt u t and u θ = dθ dτ = dθ dt u t , we have 1 2 dθ dt In the equatorial plane, to examine the solution of vertical epicyclic and radial frequencies nearby circular orbits for small perturbations δr and δθ.We shall take time derivative of Eqs.( 35) and (37), which yields The equation of motion for a perturbed particle is where ef f , and δr = r − r 0 is perturbation radius that is displaced from its original radius r = r 0 and dot indicates the time derivative.Also δθ = θ − θ 0 is the deviation for perturbation in the vertical direction, we have where ef f .By Eqs.( 35) and (37) in equatorial plane, we have and where in Eq.( 40) prime indicates derivative for radial coordinate r.Basic dynamical equations for the black hole in the Einstein-SU (N )-NLSM theory are described in the following section.

IV. BASIC DYNAMICAL EQUATIONS
We analyze the basic formulation of accretion, which was determined by Babichev et al. [44,45].We consider ideal fluid, which is described by energy-momentum tensor where pressure and energy density of the fluid are represented by p and ρ also the four-velocity of the fluid elements is u µ .In the equatorial plane, the four-velocity is here τ is the proper time along the geodesic.By using normalization condition (u µ u µ = 1) and above expression, we acquire where the condition u t > 0 is for forward flow in time while for accretion (inward flow) u r < 0. For analyzing the accretion process, all the fundamental equations are determined by calculating the energy-momentum conservation equations as well as the particle-number conservation equations.Now the law of conservation for energymomentum tensor is where √ −g = r 2 sin θ, Γ expresses Christoffel symbol's of 2nd kind and (; ) is the representation for covariant differentiation.By using metric of BH, Eq.( 45) becomes after certain simplifications, we have integrating the above equation, we obtain where C 0 is the constant of integration.By applying the law of conservation adjacent to the four-velocity through u µ T µν ;ν = 0, we get Since g µν ;ν = 0 and by using normalization condition (u µ u µ = 1), we obtain as After a few simplifications and considering non-zero components of Eq. ( 51), we get by integrating the above equation, we have Here C 1 is a constant of integration.Since u r < 0, so in above equation C 1 > 0 .We compute finally where C 2 is a constant of integration and the expression for mass flux is given by Eq. ( 55) can be written as We neglect 1 √ −g ( √ −gρu θ ) ,θ from Eq. (56) because our concern is in equatorial plane.So, √ −gρu µ would be treated as constant and we have ρu r r 2 = C 3 . (57) In the above equation C 3 is the constant of integration.Now, we have to describe the characteristics of dynamical parameters, accretion rate, and critical accretion.

A. Dynamical Parameters
We consider isothermal fluids with p = kρ here k is the equation of state parameter.For isothermal fluid, the speed of sound must be constant and p ∝ ρ also by using Eqs.(53), ( 54) and (57), we can find where C 4 is a constant of integration.By inserting p = kρ in Eq. ( 58), we deduce Thus from Eq. ( 57), we get Also by using p = kρ, pressure p can be determined.

B. Mass Evolution
In astronomical observation, the mass of the BH would progressively change due to a variety of events including Hawking radiation and accreting mass from the accretion disc onto the BH.We can determine the rate of change of mass by integrating the flux of fluid over the BH surface which is Ṁ ≡ dM dt = − T r t ds where ds = √ −gdθdφ and T r t = (p + ρ)u t u r .By putting these expressions, we obtain Ṁ as follows By considering We can acquire the time evolution of the mass of the BH with initial mass M i and Eq. ( 62) can be revised as where F ≡ 4πC 1 (p + ρ) f (r ∞ ).Integration of Eq.( 63) yields where is the time of accretion.From Eq.(64), it is clear that at t = t cr the BH mass grows up in finite time up to infinity.

C. Critical Accretion
The fluid flow is static at a large distance from the BH but its movement is inward because of BH gravitational field.When the fluid moves inward it passes over the sonic point at this point fluid velocity must be equal to the speed of sound.By using Eqs.( 57) and (58), we can get and By solving the above equations, we have Also, we have and From Eqs.( 67)-(69), we obtain the following result For critical points, we take D 1 = D 2 = 0 and get also where index c stands for the value of the corresponding quantity at the critical point.The right side of Eq.( 70) is necessarily positive.The range of the critical radius is determined by the following inequality 4f (r) + rf ′ (r) > 0.
(73) By using Eq.( 59), we have which is the required equation for the speed of sound defined by c 2 s = dp dρ .

V. CIRCULAR EQUATORIAL GEODESICS
To analyze the circular motion of the test particle, we have to analyze the effective potential (Eq.( 20)) that is given by In Fig. 2a, the graph of the effective potential of the massive particles versus r is plotted.It expresses that the effective potential depends upon the angular momentum L. We observe the first extremum at L = 6, while can not see any extremum for L < 6.Also, the effective potential of the particles increases if the values of angular momentum L increase.The solid circle in Fig. 2a, denotes the position of ISCO located at r = 6.089.Moreover, effective potential V ef f comprises two extrema for larger values of the angular momentum L, where the unstable circular orbit is located at the maximum of V ef f and stable circular orbit located at the minimum of V ef f .Figure 2b, shows the profile of the effective potential V ef f of the test particles and the dimensionless radius r.In plot 2b, of effective potential V ef f , we take different values of parameter N as N = 3,N = 4, and N = 5 with fixed angular momentum L = 10.Hence the effective potential of the test particles' motion decreases by increasing the values of parameter N .The ISCO is the smallest marginally stable circular orbit in which a test particle can stably orbit around a massive object in GR.The location of the ISCO depends upon the angular spin of the central object.The ISCO plays an important role in BH accretion discs since it marks the inner edge of the disc.Now we are interested to calculate ISCO of SU (N ) of non-linear sigma model BH, we are unable to find ISCO analytically, so we can find ISCO numerically by using the general formula of ISCO which is given in Ref. [46].When N = 3, K = 1, κ = 1/40 and Λ = −0.001 the ISCO of SU (N )-NLSM BH is r isco = 6.089, and further detail is given in Table 1.
To examine the accretion process around the BH, the role of the ISCO is important but analysis of other radii is also compulsory.As it is studied earlier, the circular orbit exists whenever r > r ph .For small perturbations, the particle's motion will be unstable if r ph < r < r isco .This relation implies that particles escape to infinity or fall into the BH.The motion of the particle remains in stable circular orbits if r > r isco .We can find out marginally bound orbit r mb , photon sphere r ph , and circular orbit r circ and by using Eqs.( 25)- (27), respectively which are given by We are unable to find marginally bound orbit analytically, so we determine numerically which is given in Table 1.
In the equatorial plane, the following quantities can be calculated as In observe that as the value of N increases the specific enand specific angular momentum of the particle are also increased.Now we are interested to find the specific energy, the specific angular momentum, the angular velocity, and the angular momentum i.e (E isco , L isco , Ω isco , and l isco ) for innermost stable circular orbit but we can not determine all of these analytically so we calculate numerically which are given in the table.

A. Mass evolution and critical points
The radial velocity and the energy density of an isothermal fluid in the accretion disc of Einstein-SU (N )-NLSM BH are given by and From Eqs.(59), ( 71) and (72), we get and Figure 4 depicts the radial velocity vs radius.The fluid has zero radial velocity away from the BH and flows with subsonic speed prior to the critical points.After passing this point, in the neighborhood of the BH, the flow accelerates and enters the supersonic domain because of BH gravity.The radial velocity increases as the Einstein-SU (N )-NLSM BH parameter is N increased, as shown in Fig. 4. The Fig. 5, demonstrates the density curve of the fluid surrounding the BH for various values of N .The density increases when the value of parameter N is increased.Furthermore, for isothermal fluids, the mass of the BH varies over time according to the follow- ing relation which shows that the mass accretion rate for Einstein-SU (N )-NLSM BH is quite different from Schwarzschild BH.The accretion rate entirely relies on the metric parameter N and the characteristics of accreting fluid onto the Einstein-SU (N )-NLSM BH.If Ṁ is positive, so (p + ρ) must also be positive.Figure 6 illustrates the graph of mass accretion rate versus r.For altered values of the parameter, N the accretion rate decreases (increases) for smaller (larger) radii.It is observed that the accretion rate raised nearby the BH as a result of immense gravitational influence.
We obtained critical accretion time and mass by using Eq. ( 64), which are given by and We observe that for ordinary flow, the mass of BH rises with the increase of parameter N .In Ref. [47] Rodrigues et al. investigated the characteristics of a Schwarzschild BH in the context of a non-minimally coupled scalar field.They discovered that for BHs with initial masses less than a specific critical value, scalar field accretion can lead to mass decrease even in the absence of phantom energy and Hawking radiations.

VI. RADIANT ENERGY FLUX
In the equatorial plane, we analyze the radiant flux from the outer layer of the accretion disc by using the values of E, L, and Ω φ .Radiation energy flux for accretion disc is determined (by Eqs.(31) and ( 32)) as follows where To examine the behavior of radiation flux of the accretion disc around the SU (N )-NLSM BH for various values of BH parameter N is plotted in Fig 7 .We observe that by increasing N , the radiation flux of the accretion disc decreases.

A. Radiation temperature
Assuming that the accretion disc must be in thermal equilibrium, one could conclude that disc emission is black body radiation, which is described by the Stefan-Boltzmann law, which relates energy flux and temperature as Here σ is the Stefan-Boltzmann constant.We analyze the temperature of the disc by using the parameter N .In Fig. 8, the profile of disc temperature is drawn for various values of N with Λ = −0.001.We examine that the temperature of the disc decreases as the value of the parameter N increases.

B. Radiative efficiency
Radiation from gravitational energy is produced as the disc's material progressively spirals toward the center.Radiative efficiency, which is the ability of the central object to transform rest mass into radiation, is a quantity that can be determined by identifying the specific energy in the ISCO radius.The formula for radiative efficiency is given by In Table 1, we display the numerical result of ISCO, marginally bound orbit r mb , specific energy, specific angular momentum, specific angular velocity, angular momentum (i.e., E 2 isco , L 2 isco , Ω 2 isco ,l 2 isco ), maximum energy flux, and maximum temperature distribution.Also, numerical investigation of all these quantities is represented in Table 2.

C. Epicyclic frequencies
A particle will undergo small oscillations in both radial and vertical directions if the particle's motion on a circular orbit is affected by perturbations in the equatorial plane.We determine radial and vertical epicyclic frequencies by using Eqs.( 40) and (41), as follows and In Fig. 9, the graph of epicyclic frequencies ver-  sus the radial parameter r is plotted.In the profile of epicyclic frequencies (Ω θ , Ω r ) along dimensionless radial parameter r, the green curve represents the behavior of vertical epicyclic frequency Ω θ which shows that the vertical frequency is independent on the parameter N and completely rely on the radial distance r.It is observed from Fig. 9 that the vertical frequency of the particles decreases by increasing the radius r.While dotted curves indicate the behavior for radial frequency.From these curves, it can be explained that the radial frequency of the particles decreases by increasing the values of the parameter N .From Fig. 10, we observe that the ratio of epicyclic frequencies is a decreasing function of r.It is also obvious that in the vicinity of BH, this ratio is very greater than unity but far from the BH it turns to unity and decreases by increasing N .

VII. CONCLUSIONS
In this article, we examine the accretion and geodesic motion of particles around the Einstein-SU (N )-NLSM BH, in the equatorial plane.These orbits have been analyzed for their stability and circular geodesics, oscillations induced by small perturbations, unstable orbits, and ultimately the general formulation for the process of accretion flow in the vicinity of the BH has been presented.Also, the BH's effective potential, specific angular momentum, typical radius, specific energy, emission rate, epicyclic frequencies, dynamical parameters, and mass evolution are found.By assuming the equation of state p = kρ for an isothermal fluid, we can find some general solutions of fluid flow onto Einstein-SU (N )-NLSM BH.The effect of SU (N )-NLSM parameter N is investigated for each aspect.The loci of the stable and unstable circular orbits are transformed to other locations under the SU (N )-NLSM parameter N influence as per our observations.As the parameter N increases, the effective potential V ef f decreases, also we locate the position of the innermost stable circular orbit in Fig. 2a.The po-sitions of the characteristic radii such as r isco , r ph , r sin , and r mb in this space-time deviate significantly from the Schwarzschild solutions.In Table 1, we analyzed the numerical results for ISCO, marginally bound, E isco , L isco , Ω isco , l isco , and radioactive efficiency.In Fig. 3, the effect of parameter N on energy and angular momentum of the BH is represented.Furthermore, as N grows, the efficiency of accretion increases.As the BH parameter N becomes larger, the radiation flux, radiant temperature, and energy efficiency of massive particles decrease.In addition to investigating circular orbits and their characteristics, epicyclic frequencies are investigated in this article.The vertical epicyclic frequency decreases monotonically with r and has no extreme.
Finally, for the case of an isothermal fluid, the nature of the radial velocity, fluid particle density, and accretion processes are investigated using the equation of state parameter k = 0.6.Apart from BH, fluids have zero radial velocity; hence, the radial velocity will be a decreasing function of r.When accretion occurs, fluid flows through a critical point where the flow rate and sound speed are equivalent.The fluid flows with subsonic speed before the critical point.The flow speed increases and then becomes supersonic after reaching that point and near the black hole due to the huge gravitational field.Finally, when the accretion rate is examined, it is found that it is affected by the metric parameter N as well as the fluid's nature.In the case of normal fluid, the mass accretion rate increases because of the gravitational effect, and its value grows even more as one gets closer to the BH.In the case of the Schwarzschild black hole, the accretion rate also increases because of the positive deviation.It is essential to note that the particle being considered in this research is not spinning.An important factor in the path of a particle is its spin.Additionally, for the sake of clarity, the viscosity and the magnetic field of the accretion disc are disregarded and the fluid is considered to be ideal.However, these variables can affect the velocity of the test particle and extension, the accretion disc's structure, and the emission rate.To that end, we plan to look into how spinning particles and accreting viscous fluids react to the Einstein-SU (N )-NLSM BH under a magnetic field.In the future it would interesting to construct the BH solution in the Einstein-SU(N)-Skyrme model to see the effects of the Skyrme term explicitly on the accretion and accretion disk.

FIG. 1 .
FIG.1.The plot between f (r) and r for distinct values of N .

Fig 3 ,
we analyze the profile of specific angular momentum and the specific energy for parameter N along with radius r of the SU (N ).In Fig 3 (a)and (b), we

FIG. 4 .FIG. 5 .
FIG. 4. The graph between u and r for certain values of N and fixed values of k and C4 are plotted

FIG. 6 .
FIG.6.The profile shows the graph for accretion rate along dimensionless radius r for various quantities of N .

FIG. 7 .
FIG. 7. The plot displays the energy flux K verses r, for altered value of N .

FIG. 8 .
FIG.8.The plot displays the relation between radiation temperature verses r, for altered value of N .

FIG. 9 .
FIG.9.The graphical representation of epicyclic frequencies and the green curve represents the behavior of vertical frequency .

FIG. 10 .
FIG.10.The plot illustrates the ratio of vertical and radial frequency .

TABLE I .
The numerical values of ISCO, marginally bound orbit r mb , specific energy, specific angular momentum, specific angular velocity, angular momentum for an innermost stable circular orbit, maximum energy flux, and maximum temperature distribution for K = 1,κ = 1 40 and Λ = −0.001.

TABLE II .
The numerical values of ISCO, marginally bound orbit r mb , specific energy, specific angular momentum, specific angular velocity, angular momentum for an innermost stable circular orbit, maximum energy flux, and maximum temperature distribution for Schwarzschild case by considering K = 0, and Λ = 0.